网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

金星上层霾消光和微物理特性的时空分布  PDF

  • 李宜奇 1,2,3
  • 孙晓兵 2,3
  • 黄红莲 2,3
  • 刘晓 2,3
  • 提汝芳 2,3
  • 郑小兵 1,2,3
  • 樊依哲 4
  • 余海啸 2,3
  • 韦祎晨 1,2,3
  • 王宇轩 2,3
  • 王宇瑶 2,3
1. 中国科学技术大学 环境科学与光电技术学院,安徽 合肥 230026; 2. 中国科学院合肥物质科学研究院 安徽光学精密机械研究所,安徽 合肥 230031; 3. 中国科学院通用光学定标与表征技术重点实验室,安徽 合肥 230031; 4. 河南理工大学,河南 焦作 454000

中图分类号: P185

最近更新:2024-11-20

DOI:10.11972/j.issn.1001-9014.2024.05.010

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

金星上层霾消光和微物理特性的变化会影响金星大气的化学和辐射平衡,为研究其时空分布,选取2006~2013年间金星快车SPICAV SOIR仪器的太阳掩星数据进行分析。首先利用MODTRAN建模去除金星中上层大气的吸收效应,后利用剥洋葱法计算67~92 km高度处金星上层霾的消光廓线。研究表明:1)金星上层霾的消光系数总体趋势是随着高度增加而减小。不同区域间的消光系数变化较大,低纬度霾在任务开始时大幅增加,昼夜间霾的平均消光系数略有变化。霾层的垂直光学厚度在10-2数量级。2)上层霾的数密度随高度增加而减小,从南极到北极,上层霾的数密度呈现先增加后减小的趋势。3)低纬度地区的云顶高度较高,为82.7 ± 5.8 km;极区的云顶高度较低,北极地区云顶高度为73.3 ± 2.4 km,南极地区为79 ± 3.5 km。北极地区上霾层平均标高为4.0 ± 0.9 km。

引言

金星有着太阳系四颗类地行星中最浓厚的大气层,其中超过96%都是CO2。其表面被厚云所覆盖。金星主云层分布在距地表47~70 km高度处,在其上方和下方有较薄的霾。早期的地面旋光观测表明,云顶的颗粒呈球形,半径约为1 µm。Venera 9号和Venera 10号着陆器以及四个先锋金星(Pioneer Venus)轨道器上的仪器可同时对金星云的垂直结构和霾粒的尺寸分布进行原位测量。Esposito等

1发现主云层由三个不同的层组成。即上层云(57~70 km)、中层云(49~57 km)和下层云(47~49 km),云层上下有弥散的霾层。硫酸气溶胶是各云层的主要成分,但各层中气溶胶的粒径分布不同。先锋金星轨道器上的光度偏振计(OCPP)测量表明,上层霾主要由非常小的“模式1”颗粒组成,其平均半径在0.25 ~0.4 µm之间。先锋金星轨道器上用于测量颗粒尺寸和形状的云颗粒尺寸光谱仪(LCPS)显示,低层的霾也主要由这些粒子组成。上层云主要由第二种粒子类型组成,称为“模式2”,其平均半径约为1 µm,但也包含大量模式1粒子。第三种粒子类型“模式3”是中层云和下层云的主要组成部分。它由相对较大的颗粒组成,平均半径约为3.85 µm,而一些颗粒的半径高达35 µm。模式3可能具有成分不确定的结晶颗粒。

了解上层霾的变化对于金星中间层的化学和辐射平衡具有重要意

2。金星的上层霾位于上层云之上,范围从云顶(70 km)到90 km1。1974年Hansen和Hovenier3发现云层内存在有效半径约为1.05 µm的颗粒。就折射率而言,这些颗粒很可能由浓硫酸组1。这种浓硫酸形成于紫外线引发的光化学反应,反应物包括SO2、CO2和H2O。紫外线分解CO2产生CO和氧原子,与大气中的SO2反应形成SO3,再与水蒸气结合生成硫酸。霾产生机制尚不清楚,Takagi4等人提出SO和SO2向上输送,经光解反应形成高海拔区域的霾。霾颗粒源以大于云层沉降速度的速度向上输送,与其他化合物反应产生SO和SO2,后通过一系列化学反应形成新的霾。先锋金星的测量表明,云层内部和云层上方存在第二种较小颗粒,其有效半径约为0.25 µm。Kawabata等5发现霾颗粒在550 nm处的折射率为1.45 ± 0.04,其有效半径reff为0.23 ± 0.04 µm,有效方差veff为0.18 ± 0.10。通过先锋金星的测量发现,在金星极地区域,波长365 nm处主云层上方霾的光学厚度约为0.83。相比之下,低纬度主云层上方霾的光学厚度为0.06。Sato等6对先锋金星2820天的数据进行了更详细的分析,发现亚微米霾呈现出较大的时间和空间变化。在北极地区,霾颗粒的有效半径为0.25 ± 0.05 µm,波长550 nm处的折射率为1.435 ± 0.02。在南极地区,霾颗粒的有效半径为0.29 ± 0.02 µm,波长550 nm处的折射率为1.45 ± 0.02。即南极地区的霾粒径比北极稍大。

2005年欧空局发射的金星快车(VEX)航天器上的几台仪器已经在广泛的波长范围内积累了金星大气层的新数据,其上的金星大气特征调查光谱仪(SPICAV/SOIR)能够研究主云层上方的霾。它具有三个工作模式:紫外(UV)、红外(IR)、太阳掩星红外(SOIR),可提供三个光谱范围(200~300 nm,0.65~1.7 µm和2.2~4.3 µm)的掩星测量。Wilquet

7等人基于2006~2010年SOIR仪器的数据,分析了金星上层霾在时间和空间上的变化,结果表明上层霾的消光系数存在短期(几个地球日)和长期(80个地球日)变化,并具有纬度依赖性。Luginin等8基于SPICAV IR数据,分析得出上层霾的气溶胶标高约为3.3 ± 0.7 km,单峰模式下霾粒子的有效半径reff = 0.54 ± 0.25 µm,并且极区的粒子半径比非极区的粒子半径小1.5~2倍。Takagi等4则利用SOIR太阳掩星观测数据研究了90 km以上金星霾层的消光特性,发现无论是低纬度还是高纬度地区,归一化消光系数都随着高度超过90 km而增加,这可能与SO和SO2气体的垂直廓线有关。

目前,国内对金星中上层大气的理论模拟较少,且缺乏完整的数据支持。本研究的重点是对2006年5月~2013年9月期间226次SPICAV SOIR太阳掩星数据进行分析,利用MODTRAN模型建立金星中上层大气透过率查找表,并对金星上层霾消光和微物理特性的时空分布进行全纬度长时间的定量分析,为后续金星中上层大气的辐射传输模拟提供数据支持。

1 SOIR仪器和数据

SOIR仪器是金星快车航天器上SPICAV/SOIR仪器套件的一个通道,其本质上是一款紧凑型星载高分辨率红外阶梯光栅光谱仪。它仅重6.5千克,体积为414 mm×254 mm×210 mm。SOIR旨在利用掩星方法探测金星大气层。Nevejans等

9、Bertaux等10和Mahieux等11详细描述了该仪器的功能和特性。这里仅进行简要描述。该仪器在红外波段(2.2~4.3 µm)下运行,对应波数范围4400 cm-1~2200 cm-1,其光谱分辨率在0.1~0.19 cm-1之间变化。AOTF的FWHM(半高全宽)约为24 cm-1。SOIR使用的衍射级从101到194不等。为了选择所需的阶梯光栅衍射级,声光可调谐滤波器AOTF位于仪器的入口狭缝之外。通过调谐AOTF并选择适当的射频激励以获得所需的衍射级波长,即可完成衍射级选择。

太阳掩星技术原理是将透过行星大气层观察太阳时测量的光谱除以行星大气层外参考光谱,以获得大气层的透过率。SOIR使用的太阳掩星技术可以推导出有关金星大气和云的垂直结构、成分等独特信息。SOIR在太阳掩星模式下运行,即仪器的视线(line of sight, LOS)始终指向太阳。当金星快车沿着其轨道移动时,LOS以连续的切线穿过金星大气层。由于航天器沿轨道运动,仪器对不同切点高度的行星大气层进行探测。可以观察到两种不同的情况:入口情况(切点高度随时间减小)和出口情况(切点高度随时间增加)。为简单起见,这里仅描述入口情况。出口情况可被视为过程相反的入口情况。图1显示了太阳掩星观测的几何结构。

图1  太阳掩星观测示意图

Fig 1  Schematic diagram of solar occultation observation

在一次掩星测量过程开始前,SOIR仪器开启,经过预冷却阶段以确保探测器温度达到约70 K的工作值,在视线穿过大气层之前开始记录太阳光谱。根据Hedin等

12的研究结果,金星大气层的高度可达 220 km。该程序开始得足够早,对应于远高于220 km的切点高度,以获得至少40个在大气层外拍摄的光谱。对低于220 km的切点高度,大气透过率的计算方法是将当前记录的光谱除以参考太阳光谱。

在一次掩星测量过程开始时,光路不会穿过金星大气层,切点高度Htg约为220 km。此时不存在大气对阳光的吸收,透过率Tobs等于1。随着金星快车仪器下行,光路越来越深入大气层,并发生两种吸收过程:金星上层霾引起的吸收和CO2等气体的吸收。随着仪器继续下行,阳光将完全被云层所遮挡,透过率约为0,此时一次掩星活动结束。

图2所示,太阳光谱已从观测光谱中去除。研究发现,无论高空观测如何,观测到的透过率的标准偏差都足够小,足以分离不同高度观测到的透过率。图3显示了本研究中使用的226次掩星测量随纬度和轨道号的分布情况。

图2  在一次掩星期间获得的光谱示例(2013年1月29日,所选衍射级order为149,轨道号为2475)。通过计算透过金星大气层看到的太阳光谱与在大气层上方测量的未衰减太阳光谱的比率来获得每个透过率。

Fig 2  Example of a spectrum obtained during an occultation (January 29, 2013, selected diffraction order order 149, orbit number 2475). Each transmittance is obtained by calculating the ratio of the solar spectrum seen through Venus' atmosphere to the unattenuated solar spectrum measured above the atmosphere.

图3  三角形(△)表示一次太阳掩星。观测到的掩星总数为226个。这些数据大部分是2006年至2013年在高纬度地区获得的。

Fig 3  The triangle (△) represents a solar occultation. The total number of observed occultations is 226. Most of these data were obtained from 2006 to 2013 at high latitudes.

2 分析方法

2.1 金星中上层大气气体透过率计算——MODTRAN建模

为了计算上层霾的消光廓线,需要从观测的透过率中去除金星大气分子(CO2、H2O、SO2等)造成的气体吸收效应。如式(1)所示,霾透过率Thaze为观测透过率Tobs除以气体透过率Tgas。其中气体透过率Tgas利用MODTRAN辐射传输模型进行计算。本研究中分子吸收计算未考虑散射效应。与可见波长相比,红外波长的光散射性较差。对于瑞利散射,在90 km高度处,CO2数密度约为1.0×1015 cm-3[

13,在3 µm处CO2的分子散射截面为5.0×10-30 cm2[14。已知金星半径为6051.8 km,上层霾分布在70~95 km范围8,上层霾以上的大气吸收微乎其微。金星大气层顶高度(TOA)可设置为100 km。由几何关系得,大气中的光路长度约为7.01×107 cm。因此,在90 km处的瑞利光学厚度估计为3.5×10-7,可忽略不计。

Thaze = TobsTgas (1)

MODTRAN

15是一种广泛使用的辐射传输(RT)代码,用于计算地球大气层的辐射传输。然而,MODTRAN中使用的辐射传输算法通常适用于任何层状大气,并且原则上可以应用于任何行星大气。Bernstein等16利用MODTRAN辐射传输模型建立了海王星上层大气模型,其中主要吸收气体是H2、CH4、C2H6和C2H2,以解释在海王星极冷的高层大气中形成的主成分为CH4和C2H6的云和霾。本研究通过修改MODTRAN模型参数,自定义大气廓线,建立适用于金星的中上层大气透过率查找表。其中金星大气的主要吸收气体是CO2,占金星大气总含量的96.5%,其吸收带位于SOIR的测量波段范围内,是建模首要考虑的吸收气体。除此之外,金星大气中含有许多微量气体,如H2O、CO、SO2、OCS、HCl等,也会对气体总透过率产生影响,建模考虑的吸收气体成分和气体混合比见表1,其中气体混合比数据来自Marcq等17的研究。

表1  金星大气吸收气体成分和混合比
Table 1  Venus atmosphere absorption gas composition and mixing ratio
金星大气层吸收气体成分Mixing Ratio
CO2 96.5%
H2O 31.18 ppmv
CO 21.53 ppmv
SO2 127.053 ppmv
HCl 0.4 ppmv
OCS 7.1 ppmv

金星的大气廓线参考VIRA

18(金星国际参考大气),从中获取金星中上层大气的温度和压强信息。建模从60 km开始,因为此时气体和气溶胶的对阳光的吸收过强,太阳光几乎不能透过金星大气。根据VIRA分层设置金星中上层大气的温度T和压强P,气体混合比在各层保持不变。采用相关k分布方法计算气体分子吸收截面,光谱分辨率0.1cm-1,小于SOIR仪器的光谱分辨率。气体参数取自MODTRAN内置的HITRAN 2008数据19。层间隔2 km,直至100 km,共21层。仪器狭缝函数设置为三角形函4,通过Voigt卷积计算得金星的中上层气体透过率。

观测几何选择如图4所示,选择“Slant Path to Space”模式,该模式适用于掩星观测,通过设置观测高度H1和切点高度Hmin来定义一条斜向光学路径。本研究中,设置H1 = 100 km,为大气建模的TOA高度。Hmin为SOIR仪器逐渐下行所测量的切点高度htg。完成这些设置后,即可建立金星的中上层大气透过率Tgas查找表。图5显示了90 km处MODTRAN计算的金星中上层气体透过率,波数范围对应于衍射级为149时SOIR仪器的波数范围。从图6可以看出,气体透过率折线应位于观测透过率折线上方,计算霾透过率时应尽量避开气体强吸收带,以减小误差。

图4  “Slant Path to Space”观测几何示意图

Fig 4  “Slant Path to Space” observation geometry diagram

图5  90 km处MODTRAN计算的金星中上层气体透过率(order 149)

Fig 5  Gas transmittance in the middle and upper layers of Venus calculated by MODTRAN at 90 km (order 149)

图6  气体透过率(Tgas,红色折线)、观测透过率(Tobs,蓝色折线)和霾透过率(Thaze,黄色折线)的示例,对应2013年1月17日获得的观测结果。(order 149,orbit 2463)

Fig 6  Examples of gas transmittance Tgas red broken line), observed transmittance (Tobs blue broken line) and haze transmittance (Thaze yellow broken line), corresponding to the observation results obtained on Jan 17, 2013. (order 149, orbit 2463)

2.2 金星上层霾消光特性的反演

每次太阳掩星测量对应特定的切点高度。如图1所示,在入口情况下,掩星活动期间切点高度htg随着时间的推移而减小。第j次测量对应切点高度htg,j。两次连续测量的切点高度htg,j-1htg之间的区域定义了大气层第j层(Lj)。τj对应于第j次观测路径上的水平光学厚度,光学路径包括层Lj和所有位于第j层上的大气层,定义为:

τj = -ln (Thaze,j) (2)

式(2)中,Thaze,j为(1)式中的Thaze,对应第j条观测路径。τj,in是第j层(Lj)的水平光学厚度,可以由剥洋葱法计算得出:

τj,in = τj-i = 1j-1 2dxi×ki (3)

式(3)中,ki是第i层的气溶胶消光系数,dxi是第i层的水平路径长度,Li层位于Lj层外部。由于切点两侧的对称性,第j层的局部消光系数kj可以通过τj,in除以2倍的第j层局部光学路径长度xj得到:

kj = τj,in2dxj (4)

重复上述过程,可以计算得出各层的气溶胶消光系数,即获取上层霾的消光廓线。图6展示了单轨消光廓线的计算结果。消光系数误差由式(5~7)确定:

Δτj = ΔThaze,jThaze,j (5)
Δτj,in = (Δτj)2+(i = 1j-1(2dxi×Δki))2 (6)
Δkj = Δτj,in2dxj (7)

式(5)中,ΔThaze,j为测量误差,该项可以从SOIR观测数据中获得,其范围在0.02%~0.12%,而由于误差通过式(6)进行逐层累积,Δkj的值会比ΔThaze,j大很多,约为35%。霾层的垂直光学厚度利用公式(8)计算:

τvert = kjdz (8)

图7  2011年6月29日(order 149 orbit 1895)剥洋葱法反演的消光廓线示例,消光系数平均误差为35.24%

Fig 7  Example of the extinction profile inverted by the onion peeling method on June 29, 2011 (order 149 orbit 1895), the average error of the extinction coefficient is 35.24%

2.3 金星上层霾微物理特性的反演

反映金星上层霾微物理特性的参数,如霾粒子的数密度等,可由2.2节获得的消光系数ki推导得出。在切点高度htg处,波长为λ的消光系数ki(λ,z)式(9)给出:

ki(λ,htg) = N(htg) σext(λ,r)n(r,htg)dr (9)

其中,σext(λ,r)是消光截面,N(htg)是粒子数密度,nr,htg是粒子谱分布。将霾粒子实际考虑为球形的气溶胶颗粒,其消光截面根据经典的Mie理论进行建模计算。本研究从 Hansen和Hovenier

3的研究结果中选择75% H2SO4的浓度值,对应的折射率的实部和虚部源值自Hummel20等人的研究结果,涵盖了SOIR仪器的光谱范围,如图8所示。

图8  75% H2SO4 的折射率在温度 T = 220 K 时的复折射率。SPICAV SOIR 光谱范围用阴影表示

Fig 8  Complex refractive index of the refractive index of 75% H2SO4 at temperature T = 220 K. SPICAV SOIR spectral range shaded

有几种分布可用于描述云和霾颗粒的尺寸分布。Hansen和Travis

21研究了使用不同尺寸分布的非偏振入射光单次散射对偏振的影响。他们发现测试的几种尺寸分布之间只有很小的差异。Hansen和Hovenier3利用伽马分布来描述金星云和霾颗粒,而Marcq等22则采用对数正态分布来描述。本研究中霾粒子谱分布采用对数正态单峰分布,其表达式如式(10)所示。根据Luginin等8的研究结果,上层霾中“模式2”粒子数密度很小,在90 km高度处仅0.1cm-3,可忽略不计,故假设上层霾中只存在“模式1”粒子,根据Kawabata等5的研究结果设置有效半径reff = 0.23 µm,有效方差veff = 0.18,有效半径reff与平均半径veff的转换关系及有效方差veff和标准差σg的转换关系如式(11)所示。

n(r) = constrexp(-(ln r-ln rg0)22ln2σg0) (10)
reff = rgexp(52ln2σg),νeff = exp(ln2σg)-1 (11)

在确定粒子谱分布之后,即可计算霾粒子的数密度,如式(12)所示:

N(htg) = ki(λ,htg)0 σext(λ,htg)n(r,htg)(dr) (12)

金星云顶高度定义为单位光学厚度处的大气高

2,即水平光学厚度τj = 1处的切点高度。霾层的标高可通过拟合指数函数来计算,如式(13)所示:

ki(λ,htg) = Ce-htg/H (13)

式(13)中,C为常数,htg为切点高度,H为霾层标高。

3 结果与讨论

八年的金星快车任务揭示了上层霾的短期和长期变化。SOIR仪器光谱范围为2.3~4.4 µm,对应衍射级101~190。Wilquet等

7发现消光系数几乎没有表现出光谱依赖性,即相邻衍射级之间消光系数非常相近。本文选取了2006~2013年8年间衍射级为149的SOIR数据,对应波长约为3 µm。在该衍射级附近观测数据较充分,利于数据分析处理。计算的消光系数结果如图9所示。为了总结规律,将研究区域划分为北极地区(60 °N~90 °N)、南极地区(60 °S~90 °S)、中纬度地区(30 °N~60 °N和30 °S~60 °S)、低纬度地区(30 °S~30 °N)。研究的金星大气层高度从67 km到92 km,涵盖了上层霾分布的范围。研究内容包括金星上层霾的消光系数、光学厚度、数密度、标高以及云顶高度随纬度和时间的变化。

图9  2006~2013年金星上层霾的消光廓线(蓝色折线代表北极地区,绿色折线代表南极地区,橙色折线表中纬度地区,红色折线代表低纬度地区,衍射级为149):图(a)~(h)分别表示2006~2013年金星上层霾的消光廓线

Fig 9  Extinction profile of the upper haze on Venus from 2006 to 2013 (the blue polyline represents the Arctic region, the green polyline represents the Antarctic region, the orange polyline represents the mid-latitude region, and the red polyline represents the low-latitude region, the diffraction order is 149):Figs (a) ~ (h) represent the extinction profiles of Venus’ upper haze from 2006 to 2013 respectively.

3.1 金星上层霾消光特性的时空分布

图9展示了8年间上层霾的消光廓线。从图9可看出,SOIR仪器在北极地区的观测比较密集,2006年由于金星快车刚进入轨道,故观测数较少。2007~2009年观测数较多,其余年份观测数大致相当。金星上层霾的消光系数总体趋势是随着高度增加,霾的消光系数减少。部分轨道消光系数并非随着高度增加而递减,消光廓线会在83 km等高度附近折叠。在北极地区高度70 km高度处,霾的消光系数约为10-2 km-1,而高度升至90 km时,消光系数约为10-4 km-1,相对降低了两个数量级。而在中低纬地区高度70 km处,消光系数约为10-2 km-1,当高度升高为90 km时,消光系数接近10-3 km-1,相对降低了一个数量级。不同区域间的消光系数变化也较大,中低纬地区消光系数总体大于极地地区消光系数,而低纬地区又略大于中纬地区,南极地区又略大于北极地区。在金星快车的前1000次绕行过程中(对应2006~2009年),低纬度地区霾的消光系数增加了约一个数量级。此外,还观察到高度90公里以上的地区也存在霾,其消光系数范围在10-4~10-5 km-1之间。

图10展示了8年间上层霾的消光系数随纬度、轨道号、高度的变化趋势。从图中可看出,8年间观测集中在北极地区,即60 °N~90 °N的范围内,主要包括前1000轨的测量。观察上层霾的消光系数随纬度的变化可发现,在相同高度下,中低纬地区的消光系数比北极地区大,上层霾的消光系数从赤道到北极地区总体呈现递减趋势,但在25 °N附近具有极大值。北极地区的消光系数从60 °N~90 °N也呈现递减趋势。南极地区观测数据较少,区域内不具有明显的趋势变化,但消光系数总体也小于中低纬地区。观察上层霾的消光系数随高度的变化可发现,消光系数随高度的增加而减少,相同高度下中低纬地区的消光系数大于极地地区,北极地区的消光系数在高度75 km处才到达10-3 km-1的量级,而中低纬地区到达该量级则在80 km以上。

图10  2006~2013年金星上层霾的消光系数随纬度、轨道号、高度的变化

Fig 10  Changes in the extinction coefficient of Venus’ upper haze with latitude, orbit number, and altitude from 2006 to 2013

为研究上层霾的消光系数在长时间尺度上的变化,且由于北极地区观测数据较充分,选择70 °N~80 °N和80 °N~90 °N的区域,分别在75 km和85 km高度处进行研究。消光系数的区域计算结果如图11所示。从图中可以看出,上层霾的消光系数随时间变化的规律不明显。在70 °N~80 °N区域内,霾的消光系数开始呈现上升趋势,在2007年中期达到最大值,之后消光系数开始减少,直至2008年底。2009~2010年消光系数呈现先上升后下降的趋势。在75 km处,平均消光系数为2.4×10-3 km-1,而在85 km处,平均消光系数为2.7×10-4 km-1,降低了约9倍。在80 °N~90 °N区域内,霾的消光系数开始也呈现上升趋势,但在2007年中有所下降,随后又有所增加。2008~2010年缓慢下降,此后霾的消光系数变化较平稳,在75 km处维持在10-3 km-1附近,在85 km处维持在10-4 km-1附近。此区域内,75 km处平均消光系数为1.5×10-3 km-1,而85 km处平均消光系数为2.0×10-4 km-1,降低了约8倍。

图11  2006~2013年金星北极地区上层霾的消光系数随时间的变化:(a) 研究区域在70 °N~80 °N,高度75 km处;(b) 研究区域在80 °N~90 °N,高度75 km处;(c) 研究区域在70 °N~80 °N,高度85 km处;(d) 研究区域在80 °N~90 °N,高度85 km处

Fig 11  Changes of the extinction coefficient with time in the upper haze in the Arctic region of Venus from 2006 to 2013:(a) The study area is between 70 °N~80 °N and an altitude of 75 km; (b) The study area is between 80 °N~90 °N and an altitude is 75 km; (c) The study area is between 70 °N~80° N, at an altitude of 85 km; (d) The study area is between 80 °N~90 °N, at an altitude of 85 km

这说明了随着高度增加,消光系数会随之减小,每上升10 km高度,平均消光系数降低约一个数量级。比较北极地区同一高度处数据可发现,纬度越高,平均消光系数越小。

此外,在金星快车任务刚开始时,低纬度的霾的数量呈现上升趋势。具体表现在前1000轨的观测中,上层霾的消光系数显而易见地增加。如图12所示,在低纬度区域(20 °S~20 °N)75 km高度处霾的消光系数从第141轨的9.8×10-4 km-1到第932轨的6.7×10-3 km-1,增大了约7倍。85 km高度处霾的消光系数从第141轨的1.5×10-4 km-1到第710轨的3.2×10-3 km-1,增大了约21倍。即2006年至2013年期间,除了低纬度霾在任务开始时大幅增加外,上层霾没有表现出明显的系统性或周期性变化。

图12  2006~2013年金星低纬度地区上层霾的消光系数随时间的变化:(a) 研究区域在20 °S~20 °N,高度75 km处;(b) 研究区域在20 °S~20 °N,高度85 km处

Fig 12  Changes in the extinction coefficient with time of the upper haze in the low latitudes of Venus from 2006 to 2013: (a) The research area is at 20 °S~20 °N, at an altitude of 75 km; (b) The research area is at 20 °S~20 °N, at an altitude of 85 km

为研究上层霾消光系数和霾层光学厚度的短期变化(如昼夜变化),将SOIR数据集划分为晨间掩星和夜间掩星两种类型,以金星本地时作为参考,观测时间在本地时6点附近的归为晨间掩星,在本地时18点附近的归为夜间掩星。图13展示了2006~2013年间226轨数据的纬度和昼夜分布。从图中可看出,晨间掩星和夜间掩星观测总数大致相当,分别为126轨和100轨,中低纬观测集中在6点和18点,即太阳升起和落下时,而高纬度观测时间则分布较广。

图13  2006~2013年间226轨数据的纬度和昼夜分布(红色圆代表晨间掩星,蓝色三角形代表夜间掩星)

Fig 13  Latitude and diurnal distribution of 226 orbit data from 2006 to 2013 (red circles represent morning occultations, blue triangles represent night occultations)

图14展示了80 km高度处金星上层霾的消光系数随昼夜和纬度的变化。可看出霾的消光系数在赤道附近较高,极区附近较低,从南极到北极消光系数总体呈现先增加后减小的趋势。晨间和夜间霾的平均消光系数分别为1.6×10-3 km-1、2.0×10-3 km-1,夜间霾的消光系数比晨间大25%。在低纬度地区,晨间和夜间霾的平均消光系数分别为3.95×10-3 km-1、3.93×10-3 km-1;在北极地区,晨间和夜间霾的平均消光系数分别为7.0×10-4 km-1、6.5×10-4 km-1;在30 °S~60 °S区域内,晨间和夜间霾的平均消光系数分别为2.1×10-3 km-1、3.8×10-3 km-1;在30 °N~60 °N区域内,晨间和夜间霾的平均消光系数分别为2.6×10-3 km-1、5.4×10-3 km-1;在南极地区,晨间和夜间霾的平均消光系数分别为1.3×10-3 km-1、2.0×10-3 km-1。即低纬度地区昼夜间霾的平均消光系数几乎没有变化,极区夜间霾的平均消光系数略大于晨间霾的平均消光系数,而中纬度夜间霾的平均消光系数约是晨间的2倍。

图14  2006~2013年间金星上层霾的消光系数随纬度和昼夜的变化(红色圆代表晨间掩星,蓝色三角形代表夜间掩星)

Fig 14  The extinction coefficient of Venus’ upper haze changes with latitude and day and night from 2006 to 2013 (red circles represent morning occultations, blue triangles represent night occultations)

在70~90 km区域内,霾层的垂直光学厚度随纬度和昼夜的变化如图15所示。可以发现霾层光学厚度的分布与消光系数的分布具有较高相似性,即霾层的垂直光学厚度在赤道附近较高,极区附近较低,从南极到北极光学厚度总体呈现先增加后减小的趋势。晨间和夜间霾层的平均垂直光学厚度分别为0.0279、0.0344,夜间霾层的平均垂直光学厚度比晨间大23%。

图15  2006~2013年间金星上层大气70~90 km处霾层垂直光学厚度随纬度和昼夜的变化(红色圆代表晨间掩星,蓝色三角形代表夜间掩星)

Fig 15  Changes in the vertical optical thickness of the haze layer at 70~90 km in the upper atmosphere of Venus with latitude and day and night from 2006 to 2013 (red circles represent morning occultations, blue triangles represent night occultations)

若以云顶高度(平均高度大于80 km)定义上霾层的起始高度,则计算的上霾层的垂直光学厚度如图16所示。可以看到,上霾层的垂直光学厚度随纬度和昼夜变化不大,数量级在10-2附近。其中晨间光学厚度为0.010 5,夜间光学厚度为0.010 2。

图16  2006~2013年间金星上霾层的垂直光学厚度随纬度和昼夜的变化(红色圆代表晨间掩星,蓝色三角形代表夜间掩星)

Fig 16  Vertical optical thickness of the haze layer on Venus changes with latitude and day and night from 2006 to 2013 (red circles represent morning occultations, blue triangles represent night occultations)

3.2 金星上层霾数密度的空间分布及敏感性分析

图17展示了2006~2013年间金星上层霾数密度的空间分布情况。可以看出,上层霾的数密度随高度增加而减小,在70 km处,上层霾的平均数密度约为330 cm-3;在80 km处,平均数密度减小至150 cm-3,在90 km处继续减小至40 cm-3。从南极到北极,上层霾的数密度呈现先增加后减小的趋势。在80 km高度处,南极地区上层霾的平均数密度约为150 cm-3,30 °S~60 °S地区平均数密度增大至160 cm-3,在低纬地区增大至180 cm-3;30 °N~60 °N地区平均数密度减小至160 cm-3,在北极地区则减小至120 cm-3

图17  2006~2013年间金星上层霾数密度的的空间分布(蓝色点代表北极地区,绿色点代表南极地区,橙色点表中纬度地区,红色点代表低纬度地区)

Fig 17  Spatial distribution of haze number density in the upper reaches of Venus from 2006 to 2013 (blue dots represent the Arctic region, green dots represent the Antarctic region, orange dots represent mid-latitude regions, and red dots represent low-latitude regions)

选取第141轨数据分析霾的数密度随有效半径reff和有效方差veff的变化趋势。从图18可以看出,霾的数密度随有效半径reff的增大而减小,随有效方差veff的增大而增大。原因是reff增加会导致粒子的消光截面增大,在相同消光系数下粒子数密度会减小;有效方差增大则会影响粒子谱分布,使得粒子谱中包含更多的小粒子,会导致总体数密度增加。霾粒子数密度对有效半径reff较敏感,即有效半径reff对霾的数密度影响较大,这一点可以从图19看出。当reff从0.15 µm增加至0.3 µm,霾的数密度从350cm-3降低至50cm-3以下。 reff = 0.15 µm时,霾粒子数密度在 veff = 0.1处为321 cm-3,在veff = 0.3处为379 cm-3。而当reff = 0.3 µm时,霾粒子数密度在 veff = 0.1处为12 cm-3,在veff = 0.3处为13 cm-3,几乎没有变化,可见霾粒子有效半径reff越大,上层霾的数密度越小,对有效方差veff的变化越不敏感。

图18  2006年9月9日(order 149 orbit 141)金星上层霾数密度随有效半径reff和有效方差veff的变化:(a) 数密度随有效半径reff的变化;(b) 数密度随有效方差veff的变化

Fig 18  Changes in haze number density in the upper haze of Venus with effective radius reff and effective variance veff on September 9, 2006 (order 149 orbit 141): (a) The number density changes with the effective radius reff; (b) The number density changes with the effective variance veff

图19  2006年9月9日(order 149 orbit 141)89.3 km处金星上层霾数密度的热图(横轴为有效半径reff,范围0.15~0.3 µm,间隔0.001 µm;纵轴为有效方差veff,范围0.1~0.3,间隔0.01)

Fig 19  Heat map of the upper haze number density of Venus at 89.3 km on September 9, 2006 (order 149 orbit 141) (the horizontal axis is the effective radius reff, the range is 0.15~0.3 µm, and the interval is 0.001 µm; the vertical axis is the effective variance veff, Range 0.1~0.3, interval 0.01)

3.3 金星云顶高度及霾层标高的时空分布

图20展示了2006~2013年间金星云顶高度随纬度的分布情况。云顶高度定义为单位光学厚度处的大气高度(τj = 1),因此与波长相关。本文研究结果为波长3 µm处的云顶高度。可以看出,低纬度地区的云顶高度较高,为82.7 ± 5.8 km;极区的云顶高度较低,北极地区云顶高度为73.3 ± 2.4 km,南极地区为79.5 ± 3.5 km。云顶高度在赤道附近最大,在30 °N附近开始减小,并在极点处达到最小值。这表明上层云在垂直方向上向极地收缩。云顶高度随年份变化特征不明显,2006~2007年在低纬地区测得的云顶高度较低,随后几年云顶高度有所增加,但在各自纬度区域内变化不大,推测与金星快车任务刚开始时低纬度霾的大量增加有关。

图20  2006~2013年间金星云顶高度随纬度的分布(年份用不同颜色表示,蓝色和红色分别表示低纬地区和极区的云顶高度,横线表示均值,竖线表示标准差)

Fig 20  Distribution of Venus cloud top height with latitude from 2006 to 2013 (years are represented by different colors, blue and red represent the cloud top height in low latitudes and polar regions respectively, the horizontal line represents the mean, and the vertical line represents the standard deviation)

由于SOIR仪器的垂直分辨率随着远离北极而变得越来越粗糙,反演仅针对北纬高纬度轨道,因为此时仪器的空间分辨率不会大于标

8。2006~2013年间金星北极地区的霾层标高随纬度的分布情况如图21所示。可以看出,北极地区的霾层标高分布在2~6 km范围内,区域内标高变化不大,根据消光廓线计算的北极地区霾层平均标高为4.0 ± 0.9 km。这表明在北极地区,上层霾分布均匀且稀薄,霾层可延伸至较高的高度。

图21  2006~2013年间北极地区金星霾层标高随纬度的分布(横线表示均值,竖线表示标准差)

Fig 21  Distribution of Venus haze layer elevation with latitude in the Arctic region from 2006 to 2013 (the horizontal line represents the mean, and the vertical line represents the standard deviation)

将各部分主要研究结果总结成表格,如表2所示。为了更好地进行对比,将相关研究成果及参考文献进行总结,如表3所示。

表2  本文主要研究结果总结
Table 2  Summary of the main research results of this article
研究内容观测数据集纬度范围研究结果
上层霾的消光系数的时空分布特征 2006~2013年SPICAV SOIR数据集 全纬度 每上升10 km高度,平均消光系数降低约一个数量级。
20 °S~20 °N 在金星快车的前1000次绕行过程中,低纬度地区霾的消光系数增加了约一个数量级。
70~90km上霾层平均垂直光学厚度 全纬度

0.0279(晨间)

0.0344(夜间)

上层霾数密度 全纬度

约340/cm3

(70km处)

约150/cm3

(80km处)

约40/cm3(90km处)

云顶高度

30 °S~30 °N

60 °N~90 °N

60 °S~90 °S

82.7 ± 5.8 km

73.3 ± 2.4 km

79.5 ± 3.5 km

霾层标高 60 °N~90 °N 4.0 ± 0.9 km
表3  相关研究成果及参考文献
Table 3  Summary of related research results and references
研究内容观测数据集纬度范围研究结果参考文献
上层霾的消光系数的时空分布特征 2007年8月19日SPICAV SOIR观测数据 70 °N 消光系数在 25 km高度范围内下降两个数量级以上。 Titov等2
2006~2010年SPICAV SOIR数据集 0°~40 °N 在金星快车的前1000次绕行过程中,低纬度地区霾的消光系数增加了一个数量级。 Wilquet7等人
70~90km上霾层平均垂直光学厚度 2006~2013年SPICAV SOIR数据集 全纬度

0.0284(晨间)

0.0351(夜间)

Takagi14
上层霾数密度 2006~2013年SPICAV IR数据集 全纬度

约500/cm3(70km处)

约50/cm3(90km处)

Luginin等8
云顶高度 2006~2010年SPICAV SOIR数据集

30 °S~30 °N

60 °N~90 °N

82 ± 1 km

73 ± 2 km

Wilquet7等人
霾层标高 2006~2013年SPICAV IR数据集 80 °N~90 °N 4.4± 1.0 km Luginin等8

4 结论

金星上层霾的时空分布对于大气中间层的化学和辐射平衡具有重要意义。选取2006~2013年间226轨SPICAV SOIR太阳掩星数据进行分析,首先利用MODTRAN进行金星中上层大气建模以去除气体吸收效应,后利用剥洋葱法计算了67 ~ 92 km 高度范围内金星上层霾的消光廓线。通过分析,得到如下结论:

(1)研究上层霾消光特性的时空分布发现,金星上层霾的消光系数总体趋势是随着高度增加,霾的消光系数减少。每上升10 km高度,平均消光系数降低约一个数量级。不同区域间的消光系数变化也较大,中低纬地区消光系数总体大于极地地区消光系数,而低纬地区又略大于中纬地区,南极地区又略大于北极地区。在金星快车的前1 000次绕行过程中(对应2006~2009年),低纬度地区霾的消光系数增加了约一个数量级。此外,还观察到高度90公里以上的地区也存在霾,其消光系数范围在10-4~10-5 km-1之间。研究上层霾的昼夜变化发现,低纬度地区昼夜间霾的平均消光系数几乎没有变化,极区夜间霾的平均消光系数略大于晨间霾的平均消光系数,而中纬度夜间霾的平均消光系数约是晨间的2倍。在70~90 km区域内,晨间和夜间霾层的平均垂直光学厚度分别为0.027 9、0.034 4。若以云顶高度定义上霾层的起始高度,上霾层的垂直光学厚度随纬度和昼夜变化不大,数量级在10-2附近。其中晨间霾层光学厚度为0.010 5,夜间霾层光学厚度为0.010 2。

(2)研究上层霾微物理特性的空间分布发现,数密度随高度增加而减小,在70 km处,上层霾的平均数密度约为340 cm-3;在80 km处,平均数密度减小至150 cm-3,在90 km处继续减小至40 cm-3。从南极到北极,上层霾的数密度呈现先增加后减小的趋势。在80 km高度处,南极地区上层霾的平均数密度为150 cm-3,30 °S~60 °S地区平均数密度增大至160 cm-3,在低纬地区增大至180 cm-3;30 °N~60 °N地区平均数密度减小至160 cm-3,在北极地区则减小至120 cm-3。霾粒子有效半径reff越大,上层霾的数密度越小,对有效方差veff的变化越不敏感。

(3)此外,还利用消光廓线反演了金星云顶高度和霾层标高。研究发现低纬度地区的云顶高度较高,为82.7 ± 5.8 km;极区的云顶高度较低,北极地区云顶高度为73.3 ± 2.4 km,南极地区为79.5 ± 3.5 km。云顶高度在赤道附近最大,在30 °N附近开始减小,并在极点处达到最小值,云顶高度随年份变化特征不明显。选取北极地区的观测数据进行分析,得出北极地区平均霾层标高为4.0 ± 0.9 km。

References

1

Esposito L WKnollenberg R GMarov M Y Aet al. 16. the clouds and hazes of venus[J]. Venus1983484. [百度学术] 

2

Titov D VIgnatiev N IMcGouldrick Ket al. Clouds and hazes of Venus[J]. Space Science Reviews20182141-61. [百度学术] 

3

Hansen J EHovenier J W. Interpretation of the polarization of Venus[J]. Journal of Atmospheric Sciences1974314): 1137-1160. [百度学术] 

4

Takagi SMahieux AWilquet Vet al. An uppermost haze layer above 100 km found over Venus by the SOIR instrument onboard Venus Express[J]. Earth, Planets and Space2019711): 1-12. [百度学术] 

5

Kawabata KCoffeen D LHansen J Eet al. Cloud and haze properties from Pioneer Venus polarimetry[J]. Journal of Geophysical Research: Space Physics198085A13): 8129-8140. [百度学术] 

6

Sato MTravis L DKawabata K. Photopolarimetry analysis of the Venus atmosphere in polar regions[J]. Icarus19961242): 569-585. [百度学术] 

7

Wilquet VDrummond RMahieux Aet al. Optical extinction due to aerosols in the upper haze of Venus: Four years of SOIR/VEX observations from 2006 to 2010[J]. Icarus20122172): 875-881. [百度学术] 

8

Luginin MFedorova ABelyaev Det al. Aerosol properties in the upper haze of Venus from SPICAV IR data[J]. Icarus2016277154-170. [百度学术] 

9

Nevejans DNeefs EVan Ransbeeck Eet al. Compact high-resolution spaceborne echelle grating spectrometer with acousto-optical tunable filter based order sorting for the infrared domain from 2.2 to 4.3 μm[J]. Applied Optics20064521): 5191-5206. [百度学术] 

10

Bertaux J LVandaele A CKorablev Oet al. A warm layer in Venus' cryosphere and high-altitude measurements of HF, HCl, H2O and HDO[J]. Nature20074507170): 646-649. [百度学术] 

11

Mahieux ABerkenbosch SClairquin Ret al.In-flight performance and calibration of SPICAV SOIR onboard Venus express[J].Applied Optics20084713):2252-65.DOI:10.1364/AO.47.002252. [百度学术] 

12

Hedin A ENiemann H BKasprzak W Tet al. Global empirical model of the Venus thermosphere[J]. Journal of Geophysical Research: Space Physics198388A1): 73-83. [百度学术] 

13

Mahieux ARobert SPiccialli Aet al. The SOIR/Venus Express species concentration and temperature database: CO2, CO, H2O, HDO, H35Cl, H37Cl, HF individual and mean profiles[J]. Icarus2023405115713. [百度学术] 

14

Seiko Takagi. Study of the Venus’ upper haze[D]. The University of Tokyo2013. [百度学术] 

15

Berk ACooley T WAnderson G Pet al. MODTRAN5: A reformulated atmospheric band model with auxiliary species and practical multiple scattering options[C]//Remote Sensing of Clouds and the Atmosphere IX. SPIE2004557178-85. [百度学术] 

16

Bernstein L SBerk ASundberg R L. Application of MODTRAN to extra-terrestrial planetary atmospheres[C]//2007 AMOS Technology ConferenceMaui, HI. 2007. [百度学术] 

17

Marcq EMills F PParkinson C Det al. Composition and chemistry of the neutral atmosphere of Venus[J]. Space Science Reviews20182141-55. [百度学术] 

18

Seiff ASchofield J TKliore A Jet al. Models of the structure of the atmosphere of Venus from the surface to 100 kilometers altitude[J]. Advances in Space Research1985511): 3-58. [百度学术] 

19

Rothman L SGordon I EBarbe Aet al. M. imekov, MAH Smith, K. Sung, SA Tashkun, J. Tennyson, RA Toth, AC Vandaele, JV Auwera[J]. J. Quant. Spectrosc. Radiat. Transf2009110533-572. [百度学术] 

20

Hummel J RShettle E PLongtin D R. AD-A210 110[J]. 1988. [百度学术] 

21

Hansen J ETravis L D. Light scattering in planetary atmospheres[J]. Space science reviews1974164): 527-610. [百度学术] 

22

Marcq EJessup K LBaggio Let al. Climatology of SO2 and UV absorber at Venus' cloud top from SPICAV-UV nadir dataset[J]. Icarus2020335113368. [百度学术]