2. 新疆大学, 新疆 乌鲁木齐 830046;
3. 中国科学院射电天文重点实验室, 江苏 南京 210008
2. Xinjiang University, Urumqi 830046, China;
3. Key Laboratory of Radio Astronomy, Chinese Academy of Sciences, Nanjing 210008, China
为获得高质量的光学星像,不仅要求望远镜在大口径的前提下保证设计完善和加工精准,而且要安装在光学观测环境好的台址。衡量台站光学观测环境质量的因素有很多,主要可分为两方面:(1) 大气因素,主要包括:大气视宁度[1]和夜天光亮度[2];(2) 气象因素,主要包括:大气水汽含量、晴夜数、温度、空气相对湿度、风速、风向。
大气视宁度是衡量台站光学观测环境质量的重要参数,是在大气湍流的影响下对星光产生的综合作用结果。目标点源星像经过圆形口径的光学系统所成的像并不是理想的像,因为衍射极限的存在使得它变成一个由中心亮斑和围绕它的衍射环组成的夫琅禾费衍射斑。星光在到达地面望远镜前需经过处于湍流状态的地球大气,大气的密度和温度时刻发生变化,导致大气折射率变化进而引起光波波前改变,此时的星像在密度、形状、位置上发生变化,不再是稳定的夫琅禾费衍射斑,其衍射环被破坏,中心亮斑变得大而模糊,整个星像不规则地迅速移动和闪烁。大气视宁度是评价望远镜台址质量的重要指标, 对台址的大气视宁度进行监测,分析视宁度随时间、季节的变化,可确保望远镜的工作效率和观测数据的质量。目前,国际上已勘选出许多观测条件优秀的天文台址,如夏威夷的Mauna kea、智利北部Paranal、西班牙La Palma,这些台址的视宁度都低于1″,其中Mauna kea甚至能达到0.5″。一般而言,如果台址一年四季的平均视宁度为1″左右[3],就可以看作是较优秀的台址。
夜天光亮度是评价光学台址观测质量的重要参数之一,直接影响观测的极限星等和观测精度,其来源一部分源于地面以上的自然光,包括大气外来的未分辨的星光、弥漫银河光、黄道光等,另一部分则是人为光源,比如城市灯光等。观测表明,太阳活动对夜天光长时标的变化也有明显的影响,即太阳活动剧烈时,夜天光亮度增加,反之减小。自然天光无法避免而光污染是人为造成的,选址中的夜天光亮度测量有助于确定天光背景,而后期对台址夜天光亮度监测,则能定量分析人为原因造成的光污染,保护天文台址。对于夜天光亮度的测量,如果V波段的夜天光亮度达到21.5~22.0 mag arcsec-2,就可以看作是优良的台址。
另外,客观存在的气象因素对观测的影响不可忽视,温度的变化能够影响近地气流变化从而影响观测质量,适宜的温度是光学观测系统正常工作的保证,相对空气湿度不大于60%以及平稳风速与风向更是开展光学观测的重要前提。
南山观测站位于乌鲁木齐县甘沟乡,海拔2080 m,当地气候属于典型温带大陆型气候,干燥少雨,夏季和冬季都无极端天气出现,常年无极端大风。自2004年开展光学观测以来,目前已有多台光学望远镜运行,如2012年安装完成的1 m大视场光学望远镜系统,用于开展光学时域天文和空间目标与碎片的观测研究。为配合光学观测工作的进行,2010年在南山站安装光学环境监测系统,包括SBIG视宁度仪①、夜天光亮度测量计(Sky quality meter-LE, SQM)②、鱼眼云量相机。另外,配合监测的设备还包括2006年安装的PC-4型自动气象站,用于测量温度、相对湿度以及风速和风向数据。本文着重对以上长期监测数据进行统计分析处理,为评估南山站的光学观测环境提供量化数据。
① https://www.optcorp.com/pdf/SBIG/SeeingMonitorManual.pdf
② http://unihedron.com/projects/sqm-le/
1 观测原理与方法 1.1 大气视宁度表征视宁度的物理量有很多。早期采用照相拖影的方法测量大气视宁度,具体是以星像拖影轨迹的抛物线拟合方差σ表征[4],σ越小表明视宁度越好。现代天文台址使用SBIG、DIMM等设备测量大气视宁度,测量结果主要以星像亮度半峰值全宽(Full Width at Half Maximum, FWHM)和大气光学湍流参数——大气相干长度r0[5]表征。测得视宁度FWHM值越小,表明视宁度越好,而r0值越小,表明视宁度越差。星像亮度半峰值全宽若忽略机械振动且望远镜达到衍射极限,与r0有如下关系:
$ FWHM = 0.98\frac{{\bar \lambda }}{{{r_0}}}, $ | (1) |
大气相干长度r0定义为:
$ {r_0} = 0.185{\lambda ^{\frac{6}{5}}}{\left( {\sec \psi } \right)^{- \frac{3}{5}}}{\left[{\int {C_n^2\left( h \right){\rm{d}}h} } \right]^{ -\frac{3}{5}}}, $ | (2) |
其中,λ为波长;ψ为天顶角;Cn2(h)为折射率结构常数的高斯分布;h为测站高度。折射率结构常数Cn2主要由温度起伏决定,即
$ {C_n} = 7.9 \times {10^{-7}}\frac{P}{{{T^2}}}{C_{\rm{T}}}, $ | (3) |
其中,CT为温度结构常数;P为气压,单位为Pa;T为温度,单位为K。
在众多影响光学天文观测的参数中,大气视宁度的测量相对复杂,测量方法也较多。下面介绍几种通过测量星像破坏和星像运动得到大气视宁度的方法。早期观测视宁度是通过Danjon[6]方法目视直接观测星像的破坏程度,并由破坏程度的分级推算大气湍流角作为大气视宁度好坏的衡量,这一方法后来被证明存在缺陷,即忽略了星像运动而只关注星像的破坏。文[7]提出的拖影方法估计大气视宁度,即通过测量星像随周日运动在底片上留下“拖影”的抖动程度测量视宁度。后来拓展为使用望远镜拖影测量,北半球则广泛使用SBIG视宁度仪以北极星为目标星进行拖影,并通过计算机程序能够直接计算得半峰值全宽,即大气视宁度。到了20世纪80年代,开发出较差Hartmann视宁度监视仪,即按照Stock原理,通过统计单星像经瞳孔平面上两个小孔径所成像的相对运动实时测量整个大气湍流强度,即大气视宁度。1986年,文[8]根据该原理研制出视宁度监测仪DIMM[9],是目前测量视宁度较为精确的方法之一。DIMM可观测(2) 式中的r0,并且不存在跟踪误差,测量结果有效性高,但目前还不能实现使用DIMM长期在无人值守状态下的监测。
南山观测站对大气视宁度的监测选用SBIG视宁度仪, 配备焦距为150 mm,焦比为F/5.3,非制冷、无快门的ST-402ME天文相机,整个装置安装在一个带有透明窗口的防雨金属盒中,窗口指向北极星。北半球以北极星为目标星,观测时目标星轨迹一直在仪器视场中,无需跟踪,由计算机程序直接实时获得天顶的视宁度半峰值全宽,长期监测能够获取大量视宁度资料进行统计分析。本文大气视宁度的观测记录时间为2010年7月至2012年1月。
1.2 夜天光亮度2011年,南山观测站开始使用Unihedron公司的夜天光亮度测量计测量在可见光波段的夜天光亮度。该设备被广泛应用于区域光污染监测工作,体积仅为3.6×2.6×1.1立方英寸,采样时间最小为1 s,最大为80 s,经校准后观测绝对精度为±10%,夜天光亮度A(magarcsec-2)其与绝对亮度F(cd/m2)存在如下转换关系:
$ A = 12.41-2.5\log F. $ | (4) |
天文台址使用夜天光亮度测量计能够获得常年的夜天光亮度监测数据,监测中不区分晴和阴以及有月夜和无月夜,为了能够通过这些数据统计获得更客观的夜天光亮度,首先需要选择全晴夜晚的观测数据,因为云、雨、雪等天气对观测结果的影响很大。其次是只统计全晴夜晚的无月夜数据。根据每天的观测记录数据可以发现,月光是可见光波段的天光亮度增加的主要原因。本文通过观测气象日报并结合当日的数据质量选测出全晴夜的观测数据,利用贝塞尔公式剔除粗大误差,即剔除天光亮度数据中过高的数据,最后只统计每日观测天光亮度的最大值均值做为台址在一段时间内的夜天光亮度水平,这样基本避免了有月夜的数据影响。本文夜天光亮度数据记录时间为2011年1月至2013年12月。
1.3 气象因素对于南山站的气象数据监测,选用技术符合气象观测规范的PC-4型自动气象站,该设备广泛应用于气象、海洋环境监测、水利等领域,可采集温度、相对湿度、风速、风向、太阳辐射等数据,工作稳定且精度高,各要素测量分辨率分别为:温度0.1℃、相对湿度0.1%、风速0.1 m/s、风向3°。2006年,在南山观测站安装了此自动气象站,该设备同时配备温度、湿度、风速、风向传感器,其中温、湿度传感器距离地面高度1.5 m,风速和风向传感器距离地面高度10 m,各项数据采集时间间隔均设置为1 min。
2 观测结果通过对南山站采集的大量视宁度资料进行统计,获得南山站大气视宁度数据的主要分布区间。使用SBIG视宁度仪每分钟能够采集数据9次,电脑记录为每分钟采集数据的中值。图 1给出南山站大气视宁度数据统计分布,共采集观测数据156061组,视宁度优于1″占5%,优于1.5″占38%,统计均值1.74″,中值1.67″。
图 2给出视宁度每月统计均值的变化情况,8月以后视宁度月均值呈变好趋势,到11月视宁度条件则突然变差,从观测当日的气象观测记录发现,2010年11月阴天居多,这使得仪器记录数据相对偏高,而2011年11月视宁度均值则只有1.37″,可以看出SBIG视宁度仪对天气因素的影响敏感,所以获得较为精确的统计数据需要选择天气情况佳的观测数据。通过曲线的整体趋势发现,秋季、冬季视宁度优于春季、夏季。视宁度月均值的最差值出现在2月为2.02″,最好值出现在10月为1.51″。
图 3给出南山观测站各月大气视宁度夜间时段的均值变化情况。具体统计各月每小时时间区间内的视宁度均值,各时段均值统计不少于25天的数据资料量。南山站的视宁度夜间变化整体较平稳,其中5月、7月、8月、9月出现黑夜与白昼转换时视宁度变差的现象。变差的原因可能是SBIG视宁度仪本身对于风的影响比较敏感,而在夜昼转换过程中温度变化比黑夜时大,造成近地空气流动,使得仪器记录的视宁度值偏大,视宁度条件变差。
根据2011年至2013年的气象观测记录以及每天夜天光亮度数据的记录情况,共选出375天全晴夜观测数据, 图 4给出这期间全晴夜夜天光亮度数据的统计分布结果,统计显示,夜天光亮度主要布在21.4~22 mag arcsec-2区间,在这375天的观测数据中,可见光波段夜天光亮度均值大于21.5 mag arcsec-2达157天。图 5给出根据气象日志筛选的全晴夜夜天光亮度监测数据的记录情况,图中显示,记录的可见光波段的夜天光亮度甚至能够达到23 mag arcsec-2,出现这样的情况是因为仪器在记录数据的过程中遇到观测天区有时有云或者冬天积雪覆盖仪器的情况,使得记录的数据偏大,对于此类数据在数据处理过程中将其作为粗大误差,通过贝塞尔公式,将大于3倍标准差的数据剔除,最后统计每天夜天光亮度的最大值均值为21.7 mag arcsec-2,另外,图中偏低的数据主要是月光对可见光波段的夜天光亮度影响造成的。
图 6给出南山站2006年5月~2016年8月的温度按月统计结果,年平均温度为3.4℃,夏季平均温度为14.1℃,记录最高温度为30.7℃,冬季平均温度为-7.8℃,记录最低温度为-31.0℃,每月平均温度与中值很接近。春季和冬季的最高温与最低温差较大为40℃左右,夏季和秋季为30℃左右。图 7给出相对湿度按月统计结果,年平均相对湿度为57.3%,夏季平均相对湿度为60.8%,冬季平均相对湿度为53.0%,每月平均相对湿度与中值很接近。南山夏季雨水较多,7月、8月相对湿度较高,分别为62.7%、60.1%,其他月份相对湿度均小于60%。
2013年7月至2016年7月通过自动气象站获得这期间的风速和风向数据,风速为每分钟所测风速的平均值。由图 8看出,这期间南山站24小时主要风向为西南风和西南风偏西,风速主要分布区间为2~4 m/s。图 9给出风速频次统计,分析表明,白天平均风速1.9 m/s,夜间平均风速1.8 m/s,出现0.2~1.5 m/s风频率占28.7%,1.5~3.3 m/s风频率占53.2%,3.3~5.4 m/s风频率占9.8%,5.4~10.7 m/s风频率出现的概率不足0.6%,无6级以上大风天气(表 1)。
风速(风级)/(m/s) | 累计百分比/% |
≤0.2 (0级) | 7.80 |
≤1.5 m/s (1级) | 36.50 |
≤3.3 m/s (2级) | 89.70 |
≤5.4 m/s (3级) | 99.40 |
≤7.9 m/s (4级) | 99.90 |
>10.7 m/s (5级) | 0.04 |
文章统计了使用SBIG视宁度仪得到南山站长期监测的视宁度结果,通过分析得到南山站视宁度均值为1.74″,中值1.67″,各季节均值有一定的变化,秋、冬季要优于春、夏季。通过夜天光亮度测量计获得南山站连续3年的夜天光亮度监测数据,统计最大值均值为21.7 mag arcsec-2,比较国内外优秀台址的夜天光亮度,南山站属于优良水平。通过自动气象站获得的历年南山站的气象资料,分析得南山站干燥、少雨,年均温度为3.4℃,相对湿度低于60%,无6级以上大风的记录。虽然在南山站使用的SBIG和夜天光亮度测量计不是目前天文界对台址观测环境监测的主流设备,但是通过其测量得到的长期连续的监测数据依然能够在一定程度上反映南山站的光学观测环境良好,适于小型光学望远镜的建设与运行。
[1] |
仇朴章, 刘忠. 现代天文选址中的视宁度[J]. 天文学进展, 1993, 11(1): 55–61 Qiu Puzhang, Liu Zhong. Seeing in modern astronomical site testing[J]. Progress in Astronomy, 1993, 11(1): 55–61. |
[2] |
宗伟凯, 付建宁, 牛家树, 等. 2009年南极冰穹A的双色天光背景统计和分析[J]. 天文研究与技术——国家天文台台刊, 2014, 11(1): 89–94 Zong Weikai, Fu Jianning, Niu Jiashu, et al. Sky brightness values in the SDSS g and r Bands at the dome A of the Antarctica in 2009[J]. Astronomical Research & Technology——Publications of National Astronomical Observatories of China, 2014, 11(1): 89–94. |
[3] | 谭徽松, 李银柱, 金振宇. 实测天体物理[M]. 北京: 国防出版社, 2014: 1-11. |
[4] | 恒星物理研究室. 滇西宾川选址总结报告(Ⅱ)——大气宁静度、大气透明度及天空背景亮度[J]. 云南天文台台刊, 1981(1): 18–34 |
[5] | Fried D L. Statistics of a geometric representation of wavefront distortion[J]. Journal of the Optical Society of America, 1965, 55(11): 1427–1431. DOI: 10.1364/JOSA.55.001427 |
[6] | Danjon A. Vistas in astronomy[M]. New York: Dover Publication, 1955: 377. |
[7] | King E S. A manual of celestial photography; principles and practice for those interested in photographing the heavens[M]. Boston: Eastern Science Supply Company, 1931: 15-17. |
[8] | Sarazin M, Roddier F. The ESO differential image motion monitor[J]. Astronomy & Astrophysics, 1990, 227(1): 294–300. |
[9] |
辛玉新, 范玉峰, 伦宝利, 等. 长期视宁度监测DIMM系统设计方案[J]. 天文研究与技术——国家天文台台刊, 2012, 9(4): 397–404 Xin Yuxin, Fan Yufeng, Lun Baoli, et al. A design of a DIMM system for long-term seeing measurement[J]. Astronomical Research & Technology——Publications of National Astronomical Observatories of China, 2012, 9(4): 397–404. |