2. 中国科学院地质与地球物理研究所 中国科学院电离层空间环境重点实验室, 北京 100029
2. Key Laboratory of Ionospheric Environment, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
1 引言
太阴潮源于月球的引力变化,因为引起太阴潮的源很清楚,尽管相对于太阳潮而言其幅度比较小,广大学者仍对太阴潮有浓厚的研究兴趣.1850年,Kreil首次从地磁日变化中发现太阴潮汐对上层大气的影响.近年来,随着雷达技术的发展,有很多利用MF雷达或者流星雷达的观测来分析中层和低热 层的大气太阴潮汐风(Stening,1989; Stening et al., 1987,1990,1994,1997,2003,Stening and Jacobi, 2000; Forbes et al., 2013).主要分析了周期为12.42太阳时的大气太阴半日潮汐(M2潮),研究结果表明,大气太阴潮汐不仅在气压场上有反应,在大气风场、地球磁场等方面也有反映,且在高层和高纬度地区比低层和低纬度地区更加明显.太阴潮中除主要的M2潮之外,与M2潮周期接近的N2潮(周期为12.66时)也引起关注.N2潮和M2潮的引力势幅度之比为0.191(Doodson,1921). N2潮最早被Bartels和Johnston在Huancayo的地磁场水平分量中确定(1940).之后的学者主要分析了磁场中的变化(Leaton et al., 1962,Winch et al., 1972)及地面压力潮汐的变化(Malin and Chapman, 1970).对于中低热层大气太阴N2潮的分析非常少,Schlapp等(1996)利用数个台站的风场数据及地磁数据分析得到N2潮与M2潮幅度之比大于0.191. 本文基于武汉流星雷达和阿德莱德中频雷达数年观测数据给出了中低纬地区大气太阴N2潮的季节、高度、年度变化及与M2潮的幅度对比分析.
2 雷达及数据分析方法 2.1 武汉流星雷达及阿德莱德中频雷达武汉流星雷达是一种全天空无线电雷达,可以用于探测约70~110 km高度上200~300 km范围内出现的流星余迹的回波.该雷达系统由发射机、多通道接收机、控制器、采集和数据分析计算机以及发射和接收天线阵构成.天线系统包含1根发射天线和5根接收天线,其布阵方式与加拿大REsoluteBay 的流星雷达基本一致(Hocking et al., 2001).通过测量不同接收天线的回波信号的相位差可以确定流星的位置,再利用回波的多普勒频移可获得流星位置处大气运动的径向速度,并研究大气运动的各种特性(熊建刚等, 2003,Niu et al., 2005,2007,Zhao et al., 2005a,2005b,Xue et al., 2007,2008).武汉 流星雷达的峰值功率为7.5 kW,工作频率为38.7 MHz,2002年9月18日以前,雷达发射的脉冲重复频率为500 Hz,其间探测的流星数目比较少,9月18日以后,脉冲重复频率升级到1980 Hz,探测的流星数目有明显的提高.
阿德莱德中频雷达是由阿德莱德大学大气物理试验组于1983年建立并开始工作.其工作频率是1.98 MHz,峰值功率最初是25 kW,从1995年开始峰值功率升级为50 kW.雷达通过一个空间天线系统接收中层大气离化的回波信号(Vincent,1986).白天和晚上的脉冲重复周期分别是12.5和25 ms.对应的高度覆盖范围分别是60~98 km和80~98 km.脉冲宽度为30μs.
2.2 数据分析方法分析太阴潮汐的方法有两种,Chapman-Miller方法(Chapman et al., 1940)和Malin和Schlapp(1980)提出的最小二乘法.最小二乘法相对Chapman-Miller方法具有一定的优点.最小二乘法不需要海量数据,对数据连续性方面的要求不高,允许数据有断点.另外,最小二乘法很容易实现误差估计.近年来,在分析太阴潮汐时被大多学者所采用.本文采用Malin和Schlapp的最小二乘法.假定风由平均风,太阳周日、半日潮汐及太阴M2(sin2τ)和N2(sin(2τ-s+p))潮汐组成.大气太阴潮汐通常幅度很小,而且很难从噪声中确定出来.因此,如果要考虑所研究的大气太阴潮汐的意义及可靠性,可信的误差估计就显得尤为重要.最小二乘法首先把数据随机地分成10个子集,每个子集分别用最小二乘法求得一组未知系数.然后对这10组系数分别求得每个系数的误差估计.如果幅度小于标准偏差,则对应的相位不可信;如果幅度大于良配的标准差,则结果可信;如果幅度在一个标准差和两个标准差之间,那么结果适度可信(Stening et al., 1994).本文采用最小二乘法对武汉(2002.2—2005.12,其中2003年3、4月份因故未曾工作)、阿德莱德(Adelaide(2002.1—2003.10))雷达观测数据进行了分析.分 别讨论了N2潮的变化及其与M2潮幅度的对比分析.
3 分析结果基于武汉流星雷达2002—2005年80~100 km的高度上将近4年的观测数据和阿得莱德2002—2003年84~98 km高度上的近两年的观测数据,我们分析了周期为12.66h的大气太阴N2潮的幅度和相位,并将其幅度和大气M2潮的幅度做了对比.表 1分别给出了上述台站的N2和M2潮的幅度,标准误差及相位,并给出了N2潮和M2潮的幅度之比.由表 1武汉风场的数据可以看到,南北向和东西向N2潮的幅度很接近,其与M2潮幅度的比值为0.383和0.368,均大于引力势之比0.191.这与 Schlapp(1996)给出的低纬Christmas Isl and (2°N,157°W)站南北向比值非常接近.阿得莱德N2潮南北方向的幅度不可信,东西向幅度强于武汉的N2潮幅度,但阿德莱德M2潮南北和东西向幅度均大于武汉的M2潮,故阿德莱德东西向N2潮和M2潮的幅度之比略小于武汉风场的结果,其值为0.295. 与阿得莱德1985—1990年的结果几乎一致(0.251). 近乎关于赤道和武汉对称的阿得莱德的M2潮幅度明显强于武汉的风场,但N2潮的幅度和武汉的非常接近.武汉和阿德莱德的风场数据均表明,N2潮与M2潮的幅度之比均大于其引力势之比的0.191,接近于引力势之比的2倍.
我们将数据分为D、E及J三个月份分析N2及M2潮的季节变化,其中D月份包括11、12、1和2月,E月份包括3、4、9和10月,J月份包括5、6、7和8月.表 2给出了武汉和阿德莱德84~98 km高度上的数据分析结果.其中,武汉J月份N2潮的南北向幅度不可信.由表 2可得:武汉的E月份,M2潮较弱,D月份(冬季)最强,J月份强度介于E、D月份之间.阿德莱德J月份较弱,D、E月份幅度很接近.N2潮武汉同样是D月份最强,J月份最弱.阿德莱德的N2潮在E月份最强,D月份最弱.M2潮和N2潮均表现出明显的季节变化特征.武汉和阿德莱德的M2和N2潮都在秋冬季较强,夏季较弱.N2与M2潮的幅度之比,武汉在J月份东西向比值为0.193,阿德莱德在D月份南北向比值为0.191,与引力势之比相符.在其余月份,N2与M2潮的幅度之比均大于0.191.
图 1给出了武汉北向太阴M2潮和N2潮在12 个月份上幅度随高度的变化.1、8月份,96 km以下 的高度上,M2潮幅度明显大于N2潮幅度.5、6、7和11月份,M2潮和N2潮幅度相近.3、10和12月份的部分高度上,N2潮的幅度甚至略大于M2潮的幅度.2月份,在84~94 km高度上,M2潮强于N2潮.武汉南北向太阴M2潮和N2潮的幅度均表现出高度和季节变化.图 2展现了武汉东向M2潮和N2潮在12个月份上幅度随高度的变化.东西向在1、7月份M2潮幅度明显大于N2潮的幅度,4、6月份在86 km高度以上,M2潮幅度明显强于N2潮,在其余月份,M2潮幅度和N2潮幅度非常接近,幅度随高度变化的趋势也十分相似,尤其是3、5、12月份的较高高度上.武汉的中低热层大气太阴M2潮和N2潮均有明显的高度和季节变化.其中在一些月份的某些高度上,N2潮的幅度甚至大于M2潮.
阿德莱德从较大尺度上来看,几乎和武汉关于赤道对称.我们采用阿德莱德2002.1—2003.10将近两年观测数据分析,阿德莱德中低热层大气太阴M2潮和N2潮幅度随季节和高度的变化特征.图 3给出了阿德莱德84~98 km高度上北向M2潮和N2潮在12个月份的幅度随高度的变化.由图 3可得:南北向分量在每个月份都有明显的幅度随高度的变化.在1、2、3、10及11月份,M2潮幅度大于N2潮的幅度,4月份的92 km高度以下、9月份的90 km高度以下,M2潮也比N2潮强.6、7、8月份,M2潮和N2潮的幅度非常接近. 6、11月份,M2潮和N2潮幅度随高度的变化趋势一致.
同时,图 4给出东西向M2潮和N2潮幅度的高度变化.在1、2、3、4、8、10月份,M2潮幅度强于N2潮.5月份96 km高度以下,6月份的88~94 km高度上,8、9月份的90~94 km高度上、11月份及12月份的86~96 km高度上,N2潮幅度略强于M2潮.6、7、8月份,M2潮和N2潮幅度大小非常接近.近乎和武汉关于赤道对称的阿得莱德中低热层大气太阴M2潮和N2潮幅度同样显示出明显的高度和季节变化特征.且在部分月份和高度上,N2潮的幅度接近甚至超过M2潮的幅度.中低热层的大气太阴N2潮值得我们关注.
武汉D月份的N2潮相对E、J月份幅度较强,我们选取D月份分析M2和N2潮的年度变化及幅度之比.表 3给出武汉D月份80~100 km高度上M2和N2潮在不同年份的幅度、标准差和相位以及N2和M2潮幅度之比.在2002—2005年的D月份中,M2潮的东西向分量均强于南北向.而N2潮除 2004年外,其余三年南北向幅度均略强于东西向.
在2002和2003两年中,N2潮的南北向幅度均强于M2潮,在其余年份,N2潮幅度均小于M2潮,但N2潮和M2潮的幅度之比均大于引力势之比.最小的比值为0.334,大于0.191.M2潮和N2潮均表现出明显的年度变化.
4 结论本文利用武汉全天空流星雷达和阿德莱德中频雷达观测的风场数据,分析了武汉和阿德莱德中低热层大气太阴M2潮和N2潮,并将N2潮和M2潮的幅度进行了对比分析.研究了大气太阴M2潮和N2潮的季节、高度及年度变化,主要结论如下:
(1)关于赤道和武汉对称的阿得莱德的M2潮幅度明显强于武汉的M2潮,但N2潮的幅度和武汉非常接近.N2潮与M2潮的幅度之比均大于其引力势之比0.191,接近于引力势之比的2倍.
(2)武汉和阿德莱德的M2潮和N2潮幅度均表现出明显的季节变化.武汉的D月份,M2潮和N2潮幅度均最强.阿得莱德在E月份最强.均在秋冬季比较强.N2与M2潮的幅度之比,武汉在J月份东西向比值为0.193,阿德莱德在D月份南北向比值为0.191,与引力势之比相符.在其余月份,N2与M2潮的幅度之比均大于0.191.阿德莱德D月份的东西向其幅度之比大于50%,比值为0.653.
(3)武汉和阿德莱德的M2潮和N2潮在大多数月份表现出明显的高度变化.在少数月份,幅度随高度的变化不明显.武汉的流星雷达有将近4年的连续观测值,我们借此分析了太阴M2和N2潮的年度变化.在不同的年份,M2潮和N2潮的幅度不同.东西向和南北向分量的强弱也表现出年度变化.在某些年份,N2潮的幅度甚至强于M2潮.且N2潮和M2潮的幅度之比均大于引力势之比.最小比值为0.334,大于0.191.
本文分析了分别位于南北半球中低纬度的阿德莱德和武汉两个雷达的观测数据,均表明,N2和M2潮的幅度之比大于引力势之比,在太阴潮汐中值得引起注意并进一步研究.
致谢 感谢由中国科学院电离层空间环境重点实验室提供武汉流星雷达的观测数据.[1] | Bartels J, Johnston H F. 1940. Geomagnetic tides in horizontal intensity at Huancayo. Terr. Magn. Atmos. Electr., 45(3): 269-308. |
[2] | Chapman S, Miller J C P. 1940. The stasistical determination of lunar daily variations in geomagnetic and meteorological elements. Mon. Not. R. A. S. Geoph. Supp., 4(9): 649-669. |
[3] | Doodson A T. 1921. The harmonic development of the tide-generating potential. Proc. R. Soc., A100: 305-329. |
[4] | Forbes J M, Zhang X, Bruinsma S, et al. 2013. Lunar semidiurnal tide in the thermosphere under solar minimum conditions. J. Geophys. Res., 118(4): 1788-1801. |
[5] | Hocking W K, Fuller B, Vandepeer B. 2001. Real-time determination of meteor-related parameters utilizing modern digital technology. J. Atmos. Terr. Solar. Phys., 63(2-3): 155-169. |
[6] | Kreil K. 1850. Einfluss des Mondes auf die magnetische Deklination. Denkschr Akad. Wiss Wein Math Naturw Kl, 3: 1-47. |
[7] | Leaton B R, Malin S R C, Finch H F. 1962. The solar and luni-solar daily variation of the geomagnetic field at Greenwich and Abinger. 1916-1957. R Obs Bull, 63. |
[8] | Malin S R C, Chapman S C. 1970. Lunar tidal components O1 and N2 in the atmospheric pressure. Pure. Appl. Geophys., 80(1): 309-318. |
[9] | Malin S R C, Schlapp D M. 1980. Geomagnetic lunar analysis by least-squares. Geophys. J. R. Astr. Soc., 60(3): 409-418. |
[10] | Niu X J, Xiong J G, Wan W X, et al. 2005. Lunar tidal winds in the mesosphere over Wuhan and Adelaide. Adv. Space. Res., 36(11): 2218-2222. |
[11] | Niu X J, Xiong J G, Wan W X, et al. 2007. A measurement of the lunar semidiurnal tide at Wuhan (30°40/N, 114°30/E). Earth Planets Space, 59(9): 991-997. |
[12] | Schlapp D M, Stening R J, Forbes J M, et al. 1996. N2 and M2 lunar tides: Atmospheric resonance revisited. Ann. Geophys., 14(8): 826-836. |
[13] | Stening R J, Meek C E, Manson A H. 1987. Lunar tidal winds measured in the upper atmosphere (78-105 km) at Saskatoon, Canada. J. Atmos. Sci., 44(8): 1143-1151. |
[14] | Stening R J. 1989. A diurnal modulation of the lunar tide in the upper atmosphere. Geophys. Res. Lerr., 16(4): 307-310. |
[15] | Stening R J, Avery S K, Tetenbaum D. 1990. Observations of lunar tides in upper atmosphere winds at Poker Flat, Alaska. J. Atmos. Terr.Phys., 52(9): 715-721. |
[16] | Stening R J, Manson A H, Meek C E, et al. 1994. Lunar tidal winds at Adelaide and Saskatoon at 80 to 100 km heights: 1985—1990. J. Geophys. Res., 99(A7): 13273-13280. |
[17] | Stening R J, Schlapp D M, Vincent R A. 1997. Lunar tides in the mesosphere over Christmas Island (2°N, 203°E). J. Geophys. Res., 102(D22): 26239-26245. |
[18] | Stening R J, Jacobi C. 2000. Lunar tidal winds in the upper atmosphere over Collm. Ann. Geophysicae, 18: 1645-1650. |
[19] | Stening R J, Tsuda T, Nakamura T. 2003. Lunar tidal winds in the upper atmosphere over Jakarta. J. Geophys. Res., 108(A5): 1192-1199. |
[20] | Vincent R A. 1986. Hardware requirements: A new generation partial reflection radar for studies of the equatorial mesosphere. MAP Handbook, 20: 85. |
[21] | Winch D E, Cunningham R A. 1972. Lunar magnetic tides at Watheroo: Seasonal, elliptic, evectional, variational and nodal components. J. Geomagn. Geoelectr., 24(4): 381-414. |
[22] | Xiong J G,Wan W X,Ning B Q,et al,2003. Meteor radar observation of circulation near mesopause over Wuhan. Chinese Science Bulletin, 48(15): 1634-1638. |
[23] | Xue X H, Wan W X, Xiong J G, et al. 2007. Diurnal tides in mesosphere/low-thermosphere during 2002 at Wuhan (30.6°N,114.4°E) using canonical correlation analysis. J. Geophys. Res., 112(D6): D06104, doi: 10. 1029/2006JD007490. |
[24] | Xue X H, Wan W X, Xiong J G, et al. 2008. The characteristics of the semi-diurnal tides in Mesos-phere/low-thermosphere (MLT) during 2002 at Wuhan (30.6°N, 114.4°E)-using canonical co-rrelation analysis technique. Adv. Space Res., 41(9): 1415-1421. |
[25] | Yu Y, Wan W X, Ning B Q, et al. 2013. Tidal Wind Mapping from observations of a meteor radar chain in December 2011. J. Geophys. Res., 118(5): 2321-2332. doi: 10. 1029/2012 JA01976. |
[26] | Zhao G X, Liu L B, Wan W X, et al. 2005. Seasonal behavior of meteor radar winds over Wuhan. Earth Planets Space, 57: 61-70. |
[27] | Zhao G X, Liu L, Wan W X, et al. 2005. The turdiurnal tide in the mesosphere and lower thermosphere over Wuhan. Earth Planets Space, 57(5): 393-398. |
[27] | 熊建刚, 万卫星, 宁百齐等. 2003. 武汉上空中层顶附近大气环流的流星雷达观测. 科学通报, 48(10): 1102-1106. |