Abstract
A model for estimating astronomical seeing at Kunlun Station (Dome A, Antarctica) is proposed. This model is based on the Tatarskii equation, using the wind shear and temperature gradient as inputs, and a seeing model depending directly on the weather data is provided. The seeing and near-ground weather data to build and validate the proposed seeing model were measured at Dome A during the summer of 2019. Two calculation methods were tested from the measured weather data relating the wind shear and temperature gradient to a combination of the two levels for the boundary layer. Both methods performed well, with correlation coefficients higher than 0.77. The model can capture the main seeing trends in which the seeing becomes small when weak wind speed and strong temperature inversion occur inside the boundary layer.
© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
“Seeing” is a critical parameter representing the optical turbulence intensity in the atmosphere, which affects the quality of astronomical images. The Antarctic Plateau is almost completely covered by snow and ice, but it has attracted interest from the astronomical community owing to its excellent seeing [1,2]. Estimating seeing in Antarctica would be significant, as many experiments of observational astronomy have been planned there [3].
Dome A (80.37° S, 77.53° E), located in East Antarctica, has the highest elevation at 4089 m on the Antarctic Plateau. On Earth, the best seeing is possible from Dome A [1,2]. Excellent seeing was directly measured by the Kunlun Differential Image Motion Monitor (KL-DIMM [4]) at Dome A during summer and winter in 2019 [2]. The optical turbulence above Dome A was mainly concentrated in the first tens of meters above the ground [5]. Similarly, at the South Pole [6] and Dome C [7] in Antarctica, only a very small fraction of the optical turbulence was found to be above the boundary layer, based on the analysis of sounding data. In addition, turbulence occurs sufficiently infrequently at high altitudes based on the computation of the Richardson number [8]; despite the occurrence of turbulence, its implications for seeing would be minor because of the relatively low air density at high altitudes [9]. Moreover, estimating the optical turbulence at high altitudes (or “free-atmosphere” seeing) has an expensive equipment cost and low data accuracy. Thus, this study focused on the boundary layer, which is a major contributor to seeing in Antarctica, and the free-atmosphere seeing was assumed to be a constant (similar to a previous report [10]).
The Monin–Obukhov similarity (MOS) theory [11] has been used to estimate the refractive index structure constant ($C_n^2$) over snow and ice from weather data [12–15]. However, the MOS theory cannot be used to estimate the optical turbulence above Dome A owing to its limited applicability under very stable conditions [16–18] because an extremely stable atmosphere was found above the snow surface of Dome A [19,20]. Applying the MOS theory was attempted at Dome A, but its estimated $C_n^2$ was unacceptably high. Fortunately, the Tatarskii equation can be used in stably stratified conditions [21].
In this study, a model for estimating the astronomical seeing at Dome A from weather data (wind speed, temperature, etc.) is proposed. An equation for seeing directly depending on weather data is given based on the Tatarskii equation. The Tatarskii equation, which was built in the context of boundary-layer turbulence [22,23], has been used to estimate the profile of $C_n^2$ in the atmosphere, and the seeing can be calculated by integrating $C_n^2$ [24,25]. In addition, the Tatarskii equation was tested at the Antarctic Taishan Station [26]. Therefore, the proposed seeing model has a rigorous physical basis.
The remainder of this paper is organized as follows. In Section 2, the experiment is described briefly, and data are provided. In Section 3, a model for estimating seeing is proposed. In Section 4, the methods for applying and evaluating the proposed model are introduced. The estimation results of the proposed model are analyzed in Section 5. In Section 6, a discussion is provided. Finally, Section 7 concludes the paper.
2. Experiment and data
The measured seeing and near-ground weather data at Dome A were acquired from previous work [2] and are available for download. The seeing was observed every minute by the KL-DIMM between January 25, 2019 and March 11, 2019. KL-DIMM can endure low temperatures (−80°C). Weather data, including temperature and wind speed, were simultaneously measured using a new second-generation Kunlun Automated Weather Station (KLAWS-2G [20]), which has a mast height of 15 m. The main technical indicators of the observational instruments are shown in Table 1, and details can be found on the website of the Kunlun Observatory, Dome A, Antarctica.
3. Theory
3.1 Optical turbulence
In this study, the boundary layer $C_n^2$ was estimated by the Tatarskii equation, which is expressed as follows [22].
where ${L_0}$ is the outer scale of turbulence, and M is the vertical gradient of the potential refractive index. The value of ${M^2}$ can be calculated from the pressure P (unit: hPa) and temperature $T$ (unit: K) [27].3.2 Boundary-layer thickness
The boundary-layer thickness (${H_{BL}}$, in m), which determines the affected scope of boundary-layer optical turbulence, generally increases with an increase in the wind speeds within the boundary layer. At Dome C, the values of ${H_{BL}}$ appear to increase linearly with the wind speed, and a linear relationship between ${H_{BL}}$ and wind speed at 3.6-m ${U_{3.6}}$ (m·s−1) can be established [30]:
The orography of the Dome A region is similar to that of Dome C (both located on a summit) [31]. Hence, it is assumed that the ${H_{BL}}$ values at Dome A could also be calculated using Eq. (6) using the measured wind speed (${U_{3.6\textrm{m}}}$ was linearly interpolated from 2 and 4 m) by KLAWS-2G. The range of the calculated ${H_{BL}}$ values coincided with the measurements in 2009 [5], with a value between 0 and 100 m.The entrance pupil of the KL-DIMM was approximately 9 m above the ground [2], which would be affected by only a portion of the boundary layer. Based on Eq. (6), the thickness of the boundary layer affecting the observation of KL-DIMM can be expressed as
3.3 Seeing model
The KL-DIMM observed seeing (${\varepsilon _{TOT}}$), which is contributed by the free atmosphere (${\varepsilon _{FA}}$) and the boundary layer (${\varepsilon _{BLA}}$), can be calculated by [32]
Here, the boundary-layer seeing can be calculated as follows (above 9 m, the height of the entrance pupil of the KL-DIMM):
The observational data of sonic radar at Dome A indicated that the optical turbulence within the boundary layer varies by only a factor of a few and then drops by several orders of magnitude within a few meters [5]. In addition, an experiment over an asphalted area showed that the vertical atmospheric turbulence was almost in a homogeneous state and was hardly affected by the temperature and wind speed changes [33]. Then, the optical turbulence intensity inside the boundary layer can be considered to be uniform in the vertical direction (similar to previous research [10]). Thus, Eq. (9) can be rewritten as
Finally, substituting Eq. (5) into Eq. (1) and substituting the resulting equation into Eq. (10) gives an equation for seeing that depends directly on the weather data.4. Method for applying and evaluating the proposed seeing model
4.1 Determination of seeing model
According to Eq. (11), the model for estimating seeing depends on S and $dT/dz$, and their values can be calculated by multiple-level measurements of KLAWS-2G. To test the robustness of the proposed model, two calculation methods were developed from the measured weather data relating the wind shear and temperature gradient to a combination of the two levels for the boundary layer:
The three undetermined coefficients ($A$, B, and $C$) in Eqs. (12a) and (12b) can be determined using a linear fitting, and the fitting data were obtained from the KL-DIMM (seeing) and KLAWS-2G (temperature and wind speed) between January 25, 2019, and February 25, 2019 (UTC). Such fitting only makes sense when the KL-DIMM is within the boundary layer; then, only the measurements under the condition of ${U_{3.6}}$ > 19/12 m·s−1 (when ${H_{BLA}}$ > 0 m) and ${\varepsilon _{TOT}}$ > 0.31 arcsec (the value of free-atmosphere seeing is assumed to be 0.31 arcsec) were used. Inspired by previous research [26], the fitting function was divided into two cases ($\partial \theta /\partial z$ < 0 and $\partial \theta /\partial z$ > 0). In addition, the Tatarskii equation shows weakness in the ability to estimate $C_n^2$ under neutral atmospheric conditions ($\partial \theta /\partial z$ ∼ 0) [26]. Therefore, in this study, the data meeting the condition $|{\partial \theta /\partial z} |$ < 0.05 K·m−1 also was not used. Finally, the coefficients ($A$, B, and $C$) in Eqs. (12a) and (12b) were fitted, as listed in Table 2.
Finally, ${\varepsilon _{BLA}}$ can be estimated by substituting the weather data into Eq. (12) using the values of A, B, and C in Table 2. Then, one can obtain ${\varepsilon _{TOT}}$ using Eq. (8), assuming ${\varepsilon _{FA}}$ = 0.31 arcsec.
4.2 Method of statistical evaluation
The statistical operators average bias ($Bias$), root mean square error ($RMSE$), and Pearson correlation coefficient (${R_{xy}}$), which have been used to assess the accuracy of the estimated astronomical seeing in previous studies [34,35], are expressed as
5. Results
5.1 Statistical analysis
To verify the prediction ability of the proposed seeing model, weather data for further eight days, measured by KLAWS-2G from February 26 through March 5, 2019 (UTC), were used. The missing measurements during these eight days were fewer. The total-atmosphere seeing values (${\varepsilon _{TOT}}$) were estimated using the model and Eqs. (8) and (12) with ${\varepsilon _{FA}}$ = 0.31 arcsec. Hence, the model estimated ${\varepsilon _{TOT}}$ cannot be smaller than 0.31 arcsec, which is one of the limitations of the model. Figures 1(a) and (b) show the scattering plots related to Eqs. (12a) and (12b), respectively. The estimated seeing using both models shows some negative bias when the measured seeing is large. Such deviations may be caused by the inaccurate estimation of free-atmosphere seeing.
Table 3 lists the $Bias$, $RMSE$, and ${R_{xy}}$ for seeing. Equation (12b) shows slightly better performance than Eq. (12a), and its ${R_{xy}}$ reaches 0.81. This could be caused by the larger height range (2 to 12 m), which provides a better representation of the boundary-layer atmosphere.
5.2 Temporal evolution
Table 3 shows Eq. (12b) using the data at 2 and 12 m high has better prediction performance; thus, the temporal evolutions of ${\varepsilon _{TOT}}$ from Eqs. (8) and (12b) are shown in Figs. 2 and 3 (four days per figure), and the corresponding temporal evolutions of wind speed and temperature are also displayed to help analyze the correlation between the weather data and the seeing. The empty area in Fig. 2(a), (c), and (e) indicate that the measured data are missing.
Figures 2 and 3 show the periods when seeing is relatively large and small, respectively. Comparing Fig. 2 with Fig. 3 reveals that a better seeing can be obtained when weak wind speed and strong temperature inversion occur. The seeing model can capture the main characteristics of the seeing trends, as shown in Figs. 2(e) and 3(e).
6. Discussion
Figure 1 and Table 3 show that the seeing models using Eqs. (12a) and (12b) perform similarly. This may be because of the stable boundary-layer atmosphere at Dome A, and the estimated optical turbulence intensity inside the boundary layer could be less likely to be fragmented, thus tending to be uniform with height [5]. The estimations show little effect from the selected levels of weather data. In this study, it seemed feasible to use the weather data (2 and 6 m) below the entrance pupil (9 m) of the KL-DIMM to estimate the seeing (${R_{xy}}$ = 0.77). The intensity of turbulence at Dome C has also been observed to be distributed almost uniformly within the boundary layer see a previoius article [30]. In summary, these results suggest that the assumption that the optical turbulence intensity inside the boundary layer is uniform is reasonable.
As discussed elsewhere [36], the temperature profile is mostly linear with height in the weakly stable regime, whereas the temperature profile in the very stable regime is exponential in shape. Thus, Figs. 2(d) and 3(d) show the mean vertical profiles of the weakly and very stable boundary layer, respectively. Comparing Figs. 2(f) and 3(f) reveals that a better seeing can be obtained in a very stable state. Infrared radiative cooling behaves like a diffusive process, giving the temperature profile an exponential shape [37]. Therefore, the low wind conditions and exponential shape of the temperature profile shown in Fig. 3 suggest very weak turbulence activity and a radiative-dominated stable boundary layer [38]. The proposed seeing model, Eq. (11), can capture such features. This is because the boundary-layer thickness (${H_{BLA}}$) is proportional to the wind speed (${U_{3.6}}$), as defined in Eq. (7), and the boundary-layer seeing value (${\varepsilon _{BLA}}$) is inversely related to the temperature gradient ($dT/dz$), as the fitted model coefficients presented in Table 2 show ($C$ is negative at $\partial \theta /\partial z$ > 0.05 K·m−1).
7. Conclusion
A model including the wind shear and temperature gradient for estimating astronomical seeing was proposed based on the Tatarskii equation. The seeing and near-ground weather data at Dome A, which were measured by the KL-DIMM and KLAWS-2G, were selected as inputs for the seeing model. Two methods of calculating the wind shear and temperature gradient were utilized to test the prediction ability of the proposed model, both of which provided good performance, with ${R_{xy}}$ values higher than 0.77.
The results show that the seeing became small when weak wind speed and strong temperature inversion occurred inside the boundary layer. The proposed seeing model can quantitatively express their relationship and achieve good estimation compared with the measurements.
To the authors’ knowledge, this study is the first to build a direct relationship between seeing and weather data over the Antarctic Plateau. Such a simplified calculation of seeing should be attempted at other sites, such as Domes C and F. In addition, the MOS theory becomes less credible when the boundary layer is under a stable atmospheric condition [16–18]. However, the air above the Antarctic Plateau is known to be a stable atmosphere [3]. In this study, the seeing model based on the Tatarskii equation was found to be a reliable alternative to the MOS theory because the Tatarskii equation has been recommended for use in stably stratified conditions [21].
Funding
Foundation of Key Laboratory of Science and Technology Innovation of Chinese Academy of Sciences (CXJJ-19S028); National Natural Science Foundation of China (41576185, 91752103).
Acknowledgments
We wish to thank Bin Ma, Zhaohui Shang, Yi Hu, et al. (National Astronomical Observatories, Chinese Academy of Sciences) for their shared data [2].
Disclosures
The authors declare no conflicts of interest.
Data availability
Data underlying the results presented in this paper are available in Ref. [2].
References
1. J. S. Lawrence, M. C. B. Ashley, A. Tokovinin, and T. Travouillon, “Exceptional astronomical seeing conditions above Dome C in Antarctica,” Nature 431(7006), 278–281 (2004). [CrossRef]
2. B. Ma, Z. Shang, Y. Hu, K. Hu, Y. Wang, X. Yang, M. C. B. Ashley, P. Hickson, and P. Jiang, “Night-time measurements of astronomical seeing at Dome A in Antarctica,” Nature 583(7818), 771–774 (2020). [CrossRef]
3. M. G. Burton, “Astronomy in Antarctica,” Astron. Astrophys. Rev. 18(4), 417–469 (2010). [CrossRef]
4. B. Ma, K. Hu, Y. Hu, W. Wang, Z. Shang, Y. Wang, and X. Yang, “An automatic DIMM for Dome A, Antarctica,” Proc. SPIE 10700, 1070052 (2018). [CrossRef]
5. C. S. Bonner, M. C. B. Ashley, X. Cui, L. Feng, X. Gong, J. S. Lawrence, D. M. Luong-Van, Z. Shang, J. W. V. Storey, L. Wang, H. Yang, J. Yang, X. Zhou, and Z. Zhu, “Thickness of the Atmospheric Boundary Layer Above Dome A, Antarctica, during 2009,” Publ. Astron. Soc. Pac. 122(895), 1122–1131 (2010). [CrossRef]
6. R. D. Marks, “Astronomical seeing from the summits of the Antarctic plateau,” Astron. Astrophys. 385(1), 328–336 (2002). [CrossRef]
7. H. Trinquet, A. Agabi, J. Vernin, M. Azouit, E. Aristidi, and E. Fossat, “Nighttime Optical Turbulence Vertical Structure above Dome C in Antarctica,” Publ. Astron. Soc. Pac. 120(864), 203–211 (2008). [CrossRef]
8. T. Travouillon, M. C. B. Ashley, M. G. Burton, J. W. V. Storey, and R. F. Loewenstein, “Atmospheric turbulence at the South Pole and its implications for astronomy,” Astron. Astrophys. 400(3), 1163–1172 (2003). [CrossRef]
9. E. Aristidi, A. Agabi, M. Azouit, E. Fossat, J. Vernin, T. Sadibekova, T. Travouillon, J. S. Lawrence, B. Halter, W. L. Roth, and V. P. Walden, “Site testing study based on weather balloons measurements,” in Dome C Astronomy and Astrophysics Meeting, M. Giard, F. Casoli, and F. Paletou, eds. (EAS Publications Series, 2004), pp. 227–232.
10. Q. Tian, P. Jiang, X. Jin, J. Li, C. Pei, F. Du, Z. Li, X. Li, H. Chen, T. Ji, X. Shi, S. Zhang, C. Yang, and H. Zhou, “Discovery of high-quality daytime seeing windows at the Antarctic Taishan station,” Mon. Not. R. Astron. Soc. 493(4), 5648–5652 (2020). [CrossRef]
11. A. S. Monin and A. M. Obukhov, “Basic Laws of Turbulent Mixing in the Atmosphere Near the Ground,” Tr. Akad. Nauk SSSR Geoph. Inst. 24(151), 163–187 (1954).
12. E. L. Andreas, “Estimating Cn2 Over Snow And Sea Ice From Meteorological data,” J. Opt. Soc. Am. A 5(4), 481–495 (1988). [CrossRef]
13. C. Qing, X. Wu, H. Huang, Q. Tian, W. Zhu, R. Rao, and X. Li, “Estimating the surface layer refractive index structure constant over snow and sea ice using Monin-Obukhov similarity theory with a mesoscale atmospheric model,” Opt. Express 24(18), 20424–20436 (2016). [CrossRef]
14. X. Wu, Q. Tian, X. Jin, P. Jiang, C. Qing, J. Cai, and H. Zhou, “Estimating optical turbulence of atmospheric surface layer at Antarctic Taishan station from meteorological data,” Acta Phys. Sin. 66(3), 039201 (2017). [CrossRef]
15. Q. Yang, X. Wu, Y. Han, and C. Qing, “Estimation of behavior of optical turbulence during summer in the surface layer above the Antarctic Plateau using the Polar WRF model,” Appl. Opt. 60(14), 4084–4094 (2021). [CrossRef]
16. A. A. Grachev, E. L. Andreas, C. W. Fairall, P. S. Guest, and P. O. G. Persson, “The Critical Richardson Number and Limits of Applicability of Local Similarity Theory in the Stable Boundary Layer,” Boundary-Layer Meteorol. 147(1), 51–82 (2013). [CrossRef]
17. D. Li, “Turbulent Prandtl number in the atmospheric boundary layer - where are we now?” Atmos. Res. 216, 86–105 (2019). [CrossRef]
18. D. L. Hutt, “Modeling and measurements of atmospheric optical turbulence over land,” Opt. Eng. 38(8), 1288–1295 (1999). [CrossRef]
19. Y. Hu, Z. Shang, M. C. B. Ashley, C. S. Bonner, K. Hu, Q. Liu, Y. Li, B. Ma, L. Wang, and H. Wen, “Meteorological Data for the Astronomical Site at Dome A, Antarctica,” Publ. Astron. Soc. Pac. 126(943), 868–881 (2014). [CrossRef]
20. Y. Hu, K. Hu, Z. Shang, M. C. B. Ashley, B. Ma, F. Du, Z. Li, Q. Liu, W. Wang, S. Yang, C. Yu, and Z. Zeng, “Meteorological Data from KLAWS-2G for an Astronomical Site Survey of Dome A, Antarctica,” Publ. Astron. Soc. Pac. 131(995), 015001 (2019). [CrossRef]
21. S. Basu, Faculty of Civil Engineering and Geosciences, Delft University of Technology, Delft, The Netherlands, and A. A. M. Holtslag are preparing a manuscript to be called “Revisiting and revising Tatarskii's formulation for the temperature structure parameter (${\$}$C_T^2${\$}$) in atmospheric flows”.
22. V. I. Tatarskii, Wave Propagation in a Turbulent Medium (McGraw-Hill Book Company, Inc., 1961).
23. E. M. Dewan and N. Grossbard, “The Inertial Range “Outer Scale” and Optical Turbulence,” Environ. Fluid Mech. 7(5), 383–396 (2007). [CrossRef]
24. S. Wu, X. Wu, C. Su, Q. Yang, J. Xu, T. Luo, C. Huang, and C. Qing, “Reliable model to estimate the profile of the refractive index structure parameter (Cn2) and integrated astroclimatic parameters in the atmosphere,” Opt. Express 29(8), 12454–12470 (2021). [CrossRef]
25. S. Wu, Q. Yang, J. Xu, T. Luo, C. Qing, C. Su, C. Huang, X. Wu, and X. Li, “A reliable model for estimating the turbulence intensity and integrated astroclimatic parameters from sounding data,” Mon. Not. R. Astron. Soc. 503(4), 5692–5703 (2021). [CrossRef]
26. Q. Yang, X. Wu, S. Wu, Y. Han, C. Su, S. Zhang, and C. Qing, “Simple method to estimate the optical turbulence over snow and ice,” J. Opt. Soc. Am. A 38(10), 1483–1488 (2021). [CrossRef]
27. C. E. Coulman, J. Vernin, Y. Coqueugniot, and J. L. Caccia, “Outer scale of turbulence appropriate to modeling refractive-indexstructure profiles,” Appl. Opt. 27(1), 155–160 (1988). [CrossRef]
28. H. Ruggiero and D. A. DeBenedictis, “Forecasting optical turbulence from mesoscale numerical weather prediction models,” in DoD High Performance Modernization Program Users Group Conference (2002), pp. 10–14.
29. H. Yang, Key Laboratory of Atmospheric Optics, Anhui Institute of Optics and Fine Mechanics, Chinese Academy of Sciences, Hefei 230031, China, and Z. Fang, C. Li, X. Deng, K. Xing, C. Xie are preparing a manuscript to be called “Atmospheric Optical Turbulence Profile Measurement and Model Improvement over Arid and Semi-arid regions”.
30. I. Petenko, S. Argentini, G. Casasanta, C. Genthon, and M. Kallistratova, “Stable Surface-Based Turbulent Layer During the Polar Winter at Dome C, Antarctica: Sodar and In Situ Observations,” Boundary-Layer Meteorol. 171(1), 101–128 (2019). [CrossRef]
31. F. Lascaux, E. Masciadri, and S. Hagelin, “Mesoscale optical turbulence simulations above Dome C, Dome A and South Pole,” Mon. Not. R. Astron. Soc. 411(1), 693–704 (2011). [CrossRef]
32. R. D. Marks, J. Vernin, M. Azouit, J. W. Briggs, M. G. Burton, M. C. B. Ashley, and J. F. Manigault, “Antarctic site testing – microthermal measurements of surface-layer seeing at the South Pole,” Astronomy and Astrophysics Supplement Series 118(2), 385–390 (1996). [CrossRef]
33. S. Rasouli and Y. Rajabi, “Investigation of the inhomogeneity of atmospheric turbulence at day and night times,” Opt. Laser Technol. 77, 40–50 (2016). [CrossRef]
34. E. Masciadri, F. Lascaux, A. Turchi, and L. Fini, “Optical turbulence forecast: ready for an operational application,” Mon. Not. R. Astron. Soc. 466(1), 520–539 (2017). [CrossRef]
35. C. Giordano, A. Rafalimanana, A. Ziad, E. Aristidi, J. Chabé, Y. Fanteï-Caujole, and C. Renaud, “Contribution of statistical site learning to improve optical turbulence forecasting,” Mon. Not. R. Astron. Soc. 504(2), 1927–1938 (2021). [CrossRef]
36. S. J. A. van der Linden, J. M. Edwards, C. C. van Heerwaarden, E. Vignon, C. Genthon, I. Petenko, P. Baas, H. J. J. Jonker, and B. J. H. van de Wiel, “Large-Eddy Simulations of the Steady Wintertime Antarctic Boundary Layer,” Boundary-Layer Meteorol. 173(2), 165–192 (2019). [CrossRef]
37. T. A. Cerni and T. R. Parish, “A Radiative Model of the Stable Nocturnal Boundary Layer with Application to the Polar Night,” J. Appl. Meteorol. Climatol. 23(11), 1563–1572 (1984). [CrossRef]
38. E. Vignon, B. J. H. van de Wiel, I. G. S. van Hooijdonk, C. Genthon, S. J. A. van der Linden, J. A. van Hooft, P. Baas, W. Maurel, O. Traullé, and G. Casasanta, “Stable boundary-layer regimes at Dome C, Antarctica: observation and analysis,” Q. J. R. Meteorol. Soc. 143(704), 1241–1253 (2017). [CrossRef]