郭商勇 胡 雄 闫召爱 程永强 涂 翠

(中国科学院国家空间科学中心 空间天气学国家重点实验室 北京 100190)

(中国科学院空间环境态势感知技术重点实验室 北京 100190)

0 引言

作为中间层与热层的边界层，中间层顶区域存在多种能量交换机制，例如太阳辐射、大气分子的红外辐射、光化学反应、重力波和潮汐波等大气动力学过程，这些因素都影响着该区域的大气温度分布特征[1]。钠原子层正好分布于中间层顶区域，因此开展钠层大气温度分布特性研究有助于更准确地掌握中高层大气规律。

20 世纪初，研究发现了存在于地球上空约80～110 km 高度范围的钠金属层，便开始把钠金属原子作为示踪物，利用其共振荧光特性获取该层大气参数。对钠原子的早期观测主要是通过分析其对太阳光的共振散射实现的。1969 年Bowman 等[2]首次利用激光雷达实现了钠层观测。随后，钠荧光激光雷达的研究不断深入，并逐步改进观测技术。早期的钠荧光激光雷达受激光器技术的限制，只能探测钠原子密度。随着激光器技术的发展，在20 世纪80 年代，采用窄线宽激光器，成功研制了可同时探测钠层大气温度和三维风场的钠荧光多普勒激光雷达[3-5]。钠荧光激光雷达具有高探测精度和高时空分辨率的特点，是可实现中间层顶大气常规观测的重要设备，其探测数据已被广泛用于地球大气特性的研究[6，7]。

本文利用中国科学院国家空间科学中心部署在河北廊坊临近空间大气探测站（39°N，116°E）的钠荧光多普勒激光雷达的探测数据，研究了2013 年廊坊上空钠层大气温度的年度和季节分布特性，并分析了影响温度分布的多种因素。

1 激光雷达系统的组成及探测数据

中国科学院国家空间科学中心于2011 年成功研制了中国首台车载钠荧光多普勒激光雷达，部署于中国科学院廊坊临近空间大气探测站开展中高层大气观测[8]。该激光雷达采用三频探测技术，发射位于钠原子D2谱线上三个频点的窄线宽脉冲激光束，激发高空钠原子产生共振荧光散射信号，通过分析后向散射信号的多普勒展宽和多普勒频移，反演出钠层大气的温度、风场及钠原子密度三个参量。廊坊钠荧光多普勒激光雷达的系统参数列于表1。

表1 廊坊钠荧光多普勒激光雷达系统参数Table 1 Langfang sodium fluorescence Doppler lidar system parameters

本文利用廊坊钠激光雷达在2013 年间单次观测时长大于3 h，共24 个夜晚、115 h 的有效探测数据，对廊坊上空的大气温度分布特性进行分析研究。

2013 年廊坊钠激光雷达探测数据在每个月份的分布情况如图1 所示。激光雷达原始数据的积分时间为1 min，积分高度76.5 m。为减小反演误差，数据处理时将每个廓线对积分时间累加为15 min、积分高度累加为306 m。日平均值由每日所有温度廓线值相加除以廓线数得出，月平均值由一个月中所有日平均值廓线累加除以天数得出，年平均值则由月平均值计算得出。

图1 2013 年廊坊钠激光雷达探测数据时间分布Fig. 1 Monthly observation hours and nights of sodium lidar for Langfang in 2013

图2 为2013 年廊坊钠原子数密度的年平均值分布廓线，廓线呈高斯状分布，半高宽约12 km，钠原子数密度最大值约为3.15×109m–3，中心高度约为91.5 km。数据结果与其他站点钠激光雷达观测结果一致[9]，表明探测数据可靠。

图2 2013 年廊坊钠原子数密度年平均值分布廓线Fig. 2 Annual mean number density profile of sodium atoms for Langfang station in 2013

2 中间层顶年均温度分布特性

图3 给出了廊坊钠激光雷达获取的83～105 km钠层大气温度年平均值廓线， 误差棒代表对应高度大气温度的标准偏差。钠层大气温度年平均分布特性可在以下三个高度范围内进行研究。（1） 83～87.5 km。大气温度随高度增加而减小，此分布特征由平流层顶一直延续至此。这是由于该高度范围内大气的热量主要来源于O3对太阳紫外辐射的吸收，而O3的密度随着高度升高而减小，则吸收的太阳辐射也随之减少，因此形成了大气温度随高度增加逐渐降低的现象[10]。

图3 2013 年钠激光雷达和SABER 探测的廊坊中间层顶区域大气温度年平均值廓线Fig. 3 Annual mean temperature profiles of Langfang mesopause region in 2013, detected by Na lidar and SABER respectively

（2） 87.5～97.5 km。大气温度先随高度增加而升高，在91 km 处达到最大值198 K，然后再次随高度增加而降低，直至中间层顶。中间层顶的年平均高度约为97.5 km，温度值约为191.2 K。在此高度范围内，存在多种能量收支形式，其中O2和O3对太阳辐射的吸收是最主要的能量来源，与CO2的红外辐射形成大气辐射平衡，在研究中常被用来计算大气背景温度。此外，放热化学反应是导致年平均温度廓线在91 km 高度处出现极大值的直接原因。为研究光化学反应，Mlynczak 等[11]建立了中间层顶区域的化学反应模型，模型包含了太阳辐射对O2和O3的光解反应、发生在该区域的Ox与Hx等粒子间的多种化学反应。结果表明，中纬度地区上空的化学反应加热率在91 km 附近达到最大值，约为10 K· d–1，表明光化学反应会导致该高度大气温度升高，这与廊坊激光雷达观测数据表现一致。

（3） 97.5 km 以上。大气温度随高度而急剧升高，在钠层顶部，达到224 K。中间层顶以上大气密度极为稀薄，能量主要来源于氧原子对太阳极紫外辐射的吸收，且受太阳活动影响较大。

数据对比表明，廊坊激光雷达的年平均温度结果与位于同纬度的美国Fort Collins （40.6°N，105°W）的钠激光雷达温度数据结果一致[12]。图3 同时给出了TIMED/SABER 卫星探测的廊坊上空温度的年平均值廓线，SABER 数据选取原则为：地理位置以廊坊钠激光雷达观测点为中心，纬度39°±2°N、经度116°±4°E 范围内；数据探测时间与钠激光雷达同日；数据平均值处理方法与钠激光雷达相同。由图3 可以看出，两条温度廓线的变化趋势保持一致，尤其在87～98 km 高度范围内，温度差维持在±3 K 范围内。87 km 以下和98 km 以上区域，两条温度廓线的差值逐渐增大，这是因为激光雷达的测量精度与钠原子密度有关，而在钠层上下边缘钠原子密度减小，导致激光雷达数据误差增大；在此高度，尤其钠层上边缘区域，稀薄的大气密度也会导致用于SABER 反演温度的波段偏离局部热力学平衡，使温度测量偏差增大[13]。

3 钠层大气温度季节分布特性

利用钠层大气温度的月平均值和温度拟合公式分析钠层大气温度的季节分布特性，拟合公式为

式中，x为 月份，y(x)为对应月份的大气温度数值，A0为 温度的年平均值,A1和A2分别为温度年变化和半年变化的幅值,φ1和φ2分别为温度年变化和半年变化的相位。

上述拟合公式考虑了温度的年变化和半年变化[14]，以月份为变量拟合出每个高度上温度变化的幅值和相位，再由拟合参数计算温度的月平均值，结果如图4 所示。由图4 可以看出中间层顶的季节变化特点，中间层顶在冬季位于99 km，最低温度181 K。冬季过后，中间层顶迅速降至夏季高度88 km，温度177 K，并在此高度持续整个夏季，然后在秋季又逐渐上升至冬季的高度。中间层顶这一季节性变化是太阳辐射和大气波动共同作用的结果。地球围绕太阳公转影响着地球对太阳辐射的接收，夏季半球接收的太阳辐射多，则地表温度较高，气流以绝热膨胀的方式从大气层底部向上传输，形成重力波，所经过区域大气温度降低。随着大气密度减小，重力波在90 km 附近不能继续向上传输，并在此区域耗散。而在100 km 以上，由于氧原子吸收了大量的太阳紫外辐射，大气温度较高。在两种因素共同作用下，使得中间层顶基本维持在90 km 以下高度。

在冬季，自对流层向上传输的重力波波谱较宽，由此引起的降温效应微弱，中间层顶区域大气温度主要受O3吸收大气辐射影响。而由于O3浓度随高度降低，大气温度也逐渐降低[15]。至100 km 以上，大气温度由于氧原子大量吸收太阳辐射而升高，因此中间层顶位于约100 km 高度以下位置。

由图4 可以看出，钠层大气气温在秋季较高。有研究认为，这是由于重力波与周日迁移潮汐发生非线性相互作用，造成周日迁移潮汐在此处破碎释放能量，从而引起温度升高的结果[16]。

图4 2013 年廊坊上空钠层大气温度月平均拟合值分布Fig. 4 Distribution of monthly average fitted values for atmospheric temperature in the sodium layer over Langfang in 2013

图5 给出了2013 年廊坊4 个季节的平均温度廓线。由图5 可明显看出，冬夏两季的中间层顶分别位于两个不同的高度，而春秋两个季节均存在两个温度低点，表现为冬夏两季温度分布的过渡。

图5 2013 年廊坊上空钠层大气温度季节平均值廓线Fig. 5 Seasonal mean temperature profiles in the sodium layer over Langfang in 2013

4 钠层大气温度的扰动特性

为研究钠层大气温度的年变化特性，由拟合公式提取出钠层大气月均温度的变化幅值和相位。2013 年廊坊钠层温度的年变化及半年变化幅值随高度分布廓线如图6 所示。由图6 可见，温度年变化的最大幅值位于87 km 高度处，约为16.5 K；最小值位于95.5 km 高度处，约为2.4 K。其中温度的半年变化幅值较小，仅在95～99 km 高度范围内与年变化幅值相当，这一特性在图4 中也显现出来，大气温度在这一高度范围表现出较为明显的半年分布特征。

图6 2013 年廊坊钠层温度年变化和半年变化的幅值廓线Fig. 6 Amplitude profiles of annual and semi-annual temperature changes in sodium layer 2013

温度拟合公式中的相位分布可用来分析影响温度变化的因素，相位对应着温度变化最大的月份。如前所述，通常将中间层顶区域大气对太阳辐射的吸收和大气分子红外辐射的综合效应作为大气背景温度，其他因素引起的温度变化作为温度扰动。分析图7中相位的年度分布可以看出，中间层顶区域上半部分的相位位于夏季，而夏季是太阳辐射最强烈的季节，表明温度变化与太阳辐射同相，太阳辐射是中间层顶上部的主要影响因素。在中间层顶区域下部，温度变化不再与太阳辐射同相，根据本文上述分析，大气动力学因素是影响此高度范围内大气温度分布的主导因素。

图7 2013 年钠层温度年变化和半年变化的相位廓线Fig. 7 Phase profiles of annual and semi-annual temperature changes in sodium layer 2013

5 结论

通过分析中国科学院廊坊临近空间大气探测站钠激光雷达2013 年间的观测数据，给出了廊坊上空钠层大气温度的年度及季节分布特征。分析了太阳辐射、大气动力学及光化学反应对中间层顶区域温度分布的影响，并从波动理论的角度分析了钠层区域大气温度的年变化和半年变化特征。由于中间层顶区域大气温度分布受多种因素影响，且探测手段不多，因此掌握的数据较少。地基钠荧光多普勒激光雷达仍然是目前、甚至未来一段时间内可连续开展中间层顶区域大气观测的重要设备。积累更多的探测数据有助于开展中间层顶区域大气温度和气候研究，也可为建立更准确的大气模型提供数据支持。