摘要
金星上层霾消光和微物理特性的变化会影响金星大气的化学和辐射平衡,为研究其时空分布,选取2006~2013年间金星快车SPICAV SOIR仪器的太阳掩星数据进行分析。首先利用MODTRAN建模去除金星中上层大气的吸收效应,后利用剥洋葱法计算67~92 km高度处金星上层霾的消光廓线。研究表明:1)金星上层霾的消光系数总体趋势是随着高度增加而减小。不同区域间的消光系数变化较大,低纬度霾在任务开始时大幅增加,昼夜间霾的平均消光系数略有变化。霾层的垂直光学厚度在1
金星有着太阳系四颗类地行星中最浓厚的大气层,其中超过96%都是CO2。其表面被厚云所覆盖。金星主云层分布在距地表47~70 km高度处,在其上方和下方有较薄的霾。早期的地面旋光观测表明,云顶的颗粒呈球形,半径约为1 µm。Venera 9号和Venera 10号着陆器以及四个先锋金星(Pioneer Venus)轨道器上的仪器可同时对金星云的垂直结构和霾粒的尺寸分布进行原位测量。Esposito等
了解上层霾的变化对于金星中间层的化学和辐射平衡具有重要意
2005年欧空局发射的金星快车(VEX)航天器上的几台仪器已经在广泛的波长范围内积累了金星大气层的新数据,其上的金星大气特征调查光谱仪(SPICAV/SOIR)能够研究主云层上方的霾。它具有三个工作模式:紫外(UV)、红外(IR)、太阳掩星红外(SOIR),可提供三个光谱范围(200~300 nm,0.65~1.7 µm和2.2~4.3 µm)的掩星测量。Wilque
目前,国内对金星中上层大气的理论模拟较少,且缺乏完整的数据支持。本研究的重点是对2006年5月~2013年9月期间226次SPICAV SOIR太阳掩星数据进行分析,利用MODTRAN模型建立金星中上层大气透过率查找表,并对金星上层霾消光和微物理特性的时空分布进行全纬度长时间的定量分析,为后续金星中上层大气的辐射传输模拟提供数据支持。
SOIR仪器是金星快车航天器上SPICAV/SOIR仪器套件的一个通道,其本质上是一款紧凑型星载高分辨率红外阶梯光栅光谱仪。它仅重6.5千克,体积为414 mm×254 mm×210 mm。SOIR旨在利用掩星方法探测金星大气层。Nevejans等
太阳掩星技术原理是将透过行星大气层观察太阳时测量的光谱除以行星大气层外参考光谱,以获得大气层的透过率。SOIR使用的太阳掩星技术可以推导出有关金星大气和云的垂直结构、成分等独特信息。SOIR在太阳掩星模式下运行,即仪器的视线(line of sight, LOS)始终指向太阳。当金星快车沿着其轨道移动时,LOS以连续的切线穿过金星大气层。由于航天器沿轨道运动,仪器对不同切点高度的行星大气层进行探测。可以观察到两种不同的情况:入口情况(切点高度随时间减小)和出口情况(切点高度随时间增加)。为简单起见,这里仅描述入口情况。出口情况可被视为过程相反的入口情况。

图1 太阳掩星观测示意图
Fig 1 Schematic diagram of solar occultation observation
在一次掩星测量过程开始前,SOIR仪器开启,经过预冷却阶段以确保探测器温度达到约70 K的工作值,在视线穿过大气层之前开始记录太阳光谱。根据Hedin等
在一次掩星测量过程开始时,光路不会穿过金星大气层,切点高度约为220 km。此时不存在大气对阳光的吸收,透过率等于1。随着金星快车仪器下行,光路越来越深入大气层,并发生两种吸收过程:金星上层霾引起的吸收和CO2等气体的吸收。随着仪器继续下行,阳光将完全被云层所遮挡,透过率约为0,此时一次掩星活动结束。
如

图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.
为了计算上层霾的消光廓线,需要从观测的透过率中去除金星大气分子(CO2、H2O、SO2等)造成的气体吸收效应。如
, | (1) |
MODTRA
金星大气层吸收气体成分 | 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 |
金星的大气廓线参考VIR
观测几何选择如

图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 气体透过率(,红色折线)、观测透过率(,蓝色折线)和霾透过率(,黄色折线)的示例,对应2013年1月17日获得的观测结果。(order 149,orbit 2463)
Fig 6 Examples of gas transmittance red broken line), observed transmittance ( blue broken line) and haze transmittance ( yellow broken line), corresponding to the observation results obtained on Jan 17, 2013. (order 149, orbit 2463)
每次太阳掩星测量对应特定的切点高度。如
, | (2) |
, | (3) |
, | (4) |
重复上述过程,可以计算得出各层的气溶胶消光系数,即获取上层霾的消光廓线。
, | (5) |
, | (6) |
, | (7) |
, | (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.2节获得的消光系数推导得出。在切点高度处,波长为的消光系数由
, | (9) |
其中,是消光截面,是粒子数密度,是粒子谱分布。将霾粒子实际考虑为球形的气溶胶颗粒,其消光截面根据经典的Mie理论进行建模计算。本研究从 Hansen和Hovenie

图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和Travi
, | (10) |
, | (11) |
在确定粒子谱分布之后,即可计算霾粒子的数密度,如
, | (12) |
金星云顶高度定义为单位光学厚度处的大气高
, | (13) |
八年的金星快车任务揭示了上层霾的短期和长期变化。SOIR仪器光谱范围为2.3~4.4 µm,对应衍射级101~190。Wilquet等

图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.

图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 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 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轨数据的纬度和昼夜分布(红色圆代表晨间掩星,蓝色三角形代表夜间掩星)
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 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 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 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)

图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轨数据分析霾的数密度随有效半径和有效方差的变化趋势。从

图18 2006年9月9日(order 149 orbit 141)金星上层霾数密度随有效半径和有效方差的变化:(a) 数密度随有效半径的变化;(b) 数密度随有效方差的变化
Fig 18 Changes in haze number density in the upper haze of Venus with effective radius and effective variance on September 9, 2006 (order 149 orbit 141): (a) The number density changes with the effective radius ; (b) The number density changes with the effective variance

图19 2006年9月9日(order 149 orbit 141)89.3 km处金星上层霾数密度的热图(横轴为有效半径,范围0.15~0.3 µm,间隔0.001 µm;纵轴为有效方差,范围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 , the range is 0.15~0.3 µm, and the interval is 0.001 µm; the vertical axis is the effective variance , Range 0.1~0.3, interval 0.01)

图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仪器的垂直分辨率随着远离北极而变得越来越粗糙,反演仅针对北纬高纬度轨道,因为此时仪器的空间分辨率不会大于标

图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)
将各部分主要研究结果总结成表格,如
研究内容 | 观测数据集 | 纬度范围 | 研究结果 |
---|---|---|---|
上层霾的消光系数的时空分布特征 | 2006~2013年SPICAV SOIR数据集 | 全纬度 | 每上升10 km高度,平均消光系数降低约一个数量级。 |
20 °S~20 °N | 在金星快车的前1000次绕行过程中,低纬度地区霾的消光系数增加了约一个数量级。 | ||
70~90km上霾层平均垂直光学厚度 | 全纬度 |
0.0279(晨间) 0.0344(夜间) | |
上层霾数密度 | 全纬度 |
约340/c (70km处) 约150/c (80km处) 约40/c | |
云顶高度 |
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 |
研究内容 | 观测数据集 | 纬度范围 | 研究结果 | 参考文献 |
---|---|---|---|---|
上层霾的消光系数的时空分布特征 | 2007年8月19日SPICAV SOIR观测数据 | 70 °N | 消光系数在 25 km高度范围内下降两个数量级以上。 |
Titov等 |
2006~2010年SPICAV SOIR数据集 | 0°~40 °N | 在金星快车的前1000次绕行过程中,低纬度地区霾的消光系数增加了一个数量级。 |
Wilque | |
70~90km上霾层平均垂直光学厚度 | 2006~2013年SPICAV SOIR数据集 | 全纬度 |
0.0284(晨间) 0.0351(夜间) |
Takag |
上层霾数密度 | 2006~2013年SPICAV IR数据集 | 全纬度 |
约500/c 约50/c |
Luginin等 |
云顶高度 | 2006~2010年SPICAV SOIR数据集 |
30 °S~30 °N 60 °N~90 °N |
82 ± 1 km 73 ± 2 km |
Wilque |
霾层标高 | 2006~2013年SPICAV IR数据集 | 80 °N~90 °N | 4.4± 1.0 km |
Luginin等 |
金星上层霾的时空分布对于大气中间层的化学和辐射平衡具有重要意义。选取2006~2013年间226轨SPICAV SOIR太阳掩星数据进行分析,首先利用MODTRAN进行金星中上层大气建模以去除气体吸收效应,后利用剥洋葱法计算了67 ~ 92 km 高度范围内金星上层霾的消光廓线。通过分析,得到如下结论:
(1)研究上层霾消光特性的时空分布发现,金星上层霾的消光系数总体趋势是随着高度增加,霾的消光系数减少。每上升10 km高度,平均消光系数降低约一个数量级。不同区域间的消光系数变化也较大,中低纬地区消光系数总体大于极地地区消光系数,而低纬地区又略大于中纬地区,南极地区又略大于北极地区。在金星快车的前1 000次绕行过程中(对应2006~2009年),低纬度地区霾的消光系数增加了约一个数量级。此外,还观察到高度90公里以上的地区也存在霾,其消光系数范围在1
(2)研究上层霾微物理特性的空间分布发现,数密度随高度增加而减小,在70 km处,上层霾的平均数密度约为340 c
(3)此外,还利用消光廓线反演了金星云顶高度和霾层标高。研究发现低纬度地区的云顶高度较高,为82.7 ± 5.8 km;极区的云顶高度较低,北极地区云顶高度为73.3 ± 2.4 km,南极地区为79.5 ± 3.5 km。云顶高度在赤道附近最大,在30 °N附近开始减小,并在极点处达到最小值,云顶高度随年份变化特征不明显。选取北极地区的观测数据进行分析,得出北极地区平均霾层标高为4.0 ± 0.9 km。
References
Esposito L W, Knollenberg R G, Marov M Y A, et al. 16. the clouds and hazes of venus[J]. Venus, 1983, 484. [百度学术]
Titov D V, Ignatiev N I, McGouldrick K, et al. Clouds and hazes of Venus[J]. Space Science Reviews, 2018, 214: 1-61. [百度学术]
Hansen J E, Hovenier J W. Interpretation of the polarization of Venus[J]. Journal of Atmospheric Sciences, 1974, 31(4): 1137-1160. [百度学术]
Takagi S, Mahieux A, Wilquet V, et al. An uppermost haze layer above 100 km found over Venus by the SOIR instrument onboard Venus Express[J]. Earth, Planets and Space, 2019, 71(1): 1-12. [百度学术]
Kawabata K, Coffeen D L, Hansen J E, et al. Cloud and haze properties from Pioneer Venus polarimetry[J]. Journal of Geophysical Research: Space Physics, 1980, 85(A13): 8129-8140. [百度学术]
Sato M, Travis L D, Kawabata K. Photopolarimetry analysis of the Venus atmosphere in polar regions[J]. Icarus, 1996, 124(2): 569-585. [百度学术]
Wilquet V, Drummond R, Mahieux A, et al. Optical extinction due to aerosols in the upper haze of Venus: Four years of SOIR/VEX observations from 2006 to 2010[J]. Icarus, 2012, 217(2): 875-881. [百度学术]
Luginin M, Fedorova A, Belyaev D, et al. Aerosol properties in the upper haze of Venus from SPICAV IR data[J]. Icarus, 2016, 277: 154-170. [百度学术]
Nevejans D, Neefs E, Van Ransbeeck E, et 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 Optics, 2006, 45(21): 5191-5206. [百度学术]
Bertaux J L, Vandaele A C, Korablev O, et al. A warm layer in Venus' cryosphere and high-altitude measurements of HF, HCl, H2O and HDO[J]. Nature, 2007, 450(7170): 646-649. [百度学术]
Mahieux A , Berkenbosch S , Clairquin R ,et al.In-flight performance and calibration of SPICAV SOIR onboard Venus express[J].Applied Optics, 2008, 47(13):2252-65.DOI:10.1364/AO.47.002252. [百度学术]
Hedin A E, Niemann H B, Kasprzak W T, et al. Global empirical model of the Venus thermosphere[J]. Journal of Geophysical Research: Space Physics, 1983, 88(A1): 73-83. [百度学术]
Mahieux A, Robert S, Piccialli A, et al. The SOIR/Venus Express species concentration and temperature database: CO2, CO, H2O, HDO, H35Cl, H37Cl, HF individual and mean profiles[J]. Icarus, 2023, 405: 115713. [百度学术]
Seiko Takagi. Study of the Venus’ upper haze[D]. The University of Tokyo, 2013. [百度学术]
Berk A, Cooley T W, Anderson G P, et al. MODTRAN5: A reformulated atmospheric band model with auxiliary species and practical multiple scattering options[C]//Remote Sensing of Clouds and the Atmosphere IX. SPIE, 2004, 5571: 78-85. [百度学术]
Bernstein L S, Berk A, Sundberg R L. Application of MODTRAN to extra-terrestrial planetary atmospheres[C]//2007 AMOS Technology Conference, Maui, HI. 2007. [百度学术]
Marcq E, Mills F P, Parkinson C D, et al. Composition and chemistry of the neutral atmosphere of Venus[J]. Space Science Reviews, 2018, 214: 1-55. [百度学术]
Seiff A, Schofield J T, Kliore A J, et al. Models of the structure of the atmosphere of Venus from the surface to 100 kilometers altitude[J]. Advances in Space Research, 1985, 5(11): 3-58. [百度学术]
Rothman L S, Gordon I E, Barbe A, et al. M. imekov, MAH Smith, K. Sung, SA Tashkun, J. Tennyson, RA Toth, AC Vandaele, JV Auwera[J]. J. Quant. Spectrosc. Radiat. Transf, 2009, 110: 533-572. [百度学术]
Hummel J R, Shettle E P, Longtin D R. AD-A210 110[J]. 1988. [百度学术]
Hansen J E, Travis L D. Light scattering in planetary atmospheres[J]. Space science reviews, 1974, 16(4): 527-610. [百度学术]
Marcq E, Jessup K L, Baggio L, et al. Climatology of SO2 and UV absorber at Venus' cloud top from SPICAV-UV nadir dataset[J]. Icarus, 2020, 335: 113368. [百度学术]