LEO-LEO红外激光掩星CO2浓度测量技术研究 下载: 1024次
1 引言
自工业革命以来,人类活动加剧,随之而来的酸雨、温室效应等环境问题逐渐显露并日趋严峻。研究表明,以CO2为主的温室气体是温室效应的主要推力,大气中CO2浓度的变化对气候的影响备受气候、环境等自然学科的重视。而高时空分辨率的CO2浓度数据是气候科学研究的重要需求,也是研究大气全球变化、人类活动对大气环境影响的重要参数,对于提高数值天气预报的精度以及研究日地能量传输过程具有重要意义。
星载积分路径差分吸收(IPDA)激光雷达具有探测精度高、测量区域广且不受气溶胶干扰等优点[1],刘继桥课题组[2-3]针对CO2 IPDA激光雷达开展了广泛研究,目前已研制出原理样机,其对CO2体积混合比柱浓度的理论测量误差可达到10-6。但是IPDA激光雷达获取的是整层大气CO2柱浓度的平均值,无法提供大气CO2的垂直变化数据。为了兼顾全球CO2浓度的测量精度和垂直分辨能力,可以考虑将IPDA激光雷达技术应用到新的大气探测方式中。通过卫星组网构成掩星星座,对地球大气进行临边探测,是获得大气参量廓线的有效手段。目前利用无线电信号进行探测的全球导航卫星系统无线电掩星技术GNSS RO(Global Navigation Satellite System Radio Occultation,GRO),已实现对大气热力学参量廓线的测量,其数据产品已被广泛应用于数值天气预报和全球气候监测[4]。鉴于GRO技术在大气热力学参量廓线测量中的成功运用,欧洲的ACCURATE(Atmospheric Climate and Chemistry in the UTLS Region And Climate Trends Explorer)[5]提出了在近地轨道(LEO)卫星构成的掩星星座上搭载红外激光发射机和接收机的设想,采用主动临边探测的方式进行大气成分和风速的测量,即激光差分吸收与掩星探测相结合的红外激光掩星(LIO)技术。Clifford等[6]也提出了类似的方案,即将搭载有激光发射机的卫星与多个空间反射镜组网进行临边探测,并针对对流层顶到平流层底高度(UTLS, 5~35 km)区域水汽的探测进行仿真研究。
LIO技术结合了激光差分吸收与掩星探测的优点,主动式激光差分吸收光谱技术可达到较高的探测精度,同时LEO-LEO掩星探测方式可以实现对全球高覆盖率的测量,获得垂直分辨能力(约1 km)的探测数据。本文主要研究了利用LIO技术测量全球CO2浓度廓线的方法。建立针对UTLS区域CO2浓度廓线的探测仿真模型,通过计算分析优化选择激光的工作波数,利用选择的工作波数对LIO技术的探测性能进行仿真分析。
2 基本原理
2.1 激光路径的确定
如
式中:
在发生掩星事件时,两颗LEO卫星的位置坐标
结合(1)~(3)式可实时对
2.2 浓度廓线的反演
on-line和off-line激光在穿过大气层时被CO2分子吸收,传输路径上的差分吸收光学厚度(DAOD)表示为
式中:
式中:Δ
对(5)式应用Abel积分变化可以得到
根据Beer-Lambert定律,
式中:
式中:
2.3 链路衰减模型
激光信号从发射机到接收机整个传输过程中的损耗主要由三部分组成:自由空间传播损耗、大气损耗和接收机光学系统损耗。由激光雷达方程可知,接收机接收到的激光信号功率为
式中:
大气衰减主要包含CO2吸收、气溶胶吸收、其他无关气体(H2O、CO2、O3、N2O、CO、CH4、O2等)吸收、大气分子瑞利散射。瑞利散射产生的等效吸收截面为
式中:
计算过程中的大气温度、湿度、压强采用美国1976标准大气模型[8]获得,红外波段的折射率采用Edlén方程[9]获得。利用HITRAN 2016数据库拟合分子的吸收截面,结合(5)式对光学厚度进行计算。
3 激光波数的优化选择
3.1 on-line波数选择
目前,学术领域针对CO2浓度的主动探测主要集中在1.5 μm和2 μm的红外波段。在临边探测方式下,考虑到激光信号穿过大气层的路径较长,在1.5 μm波段处,空气分子的瑞利散射和气溶胶米散射对信号的影响很大,而2 μm波段处受到的太阳辐射和大气热辐射的影响较小,因此拟将2 μm作为CO2浓度廓线探测的优选波段[5]。
off-line处的光学厚度很小,DAOD主要由on-line光学厚度决定。如果所选择的on-line处的CO2吸收过强,则当掩星切点高度较低时,on-line信号就会因穿过大气的路程较长而具有较大的光学厚度,从而被强烈吸收,难以测量;相反,如果on-line处的吸收过弱,则当掩星切点高度较大时,on-line信号穿过大气的路程较小,并且由于高空中的CO2含量减少,导致光学厚度太小,到达接收机的on-line和off-line信号的差别就比较微弱,不利于减小测量误差。所以,应合理选择on-line波数,使UTLS高度范围(5~35 km)具有较高的探测精度。在切点高度
式中:Δ
利用2.3节中所述链路衰减模型,分别对不同on-line波数下系统的相对随机误差进行仿真计算,得到了5 km和35 km高度处相对随机误差与波数的关系曲线,如
图 2. 5 km和35 km高度处相对随机误差随on-line波数的变化
Fig. 2. RRE as a function of on-line wavenumber at altitudes of 5 km and 35 km
由
3.2 off-line波数选择
为了达到良好的差分效果,off-line波数应选在CO2吸收尽可能弱的位置。为了消除大气宽带效应的影响,同时降低光源研制的技术难度,off-line与on-line之间的波数间隔不能过大。除此之外,还应考虑水汽对on-line和off-line的影响,在选择off-line波数时,应尽可能使on-line和off-line处水汽的光学厚度相近。由于大气中的水汽含量随海拔高度的升高而迅速减小,对CO2测量的干扰主要在低空,仿真计算切点高度为5 km时,激光信号传输路径上CO2和H2O在4771 cm-1附近的吸收光谱如
由仿真结果可知,在选定的on-line波数处,H2O的光学厚度为0.013,将产生约1.5%的浓度测量误差,影响较大。可以将off-line选在与on-line H2O光学厚度接近的4772.0240 cm-1波数处,以减小水汽的干扰。考虑到在不同高度处温度、压强的差异,水汽在不同波数处的吸收谱线将有不同程度的展宽,故5~35 km各高度上水汽对CO2差分光学厚度的干扰有所差异。水汽对CO2差分光学厚度测量结果造成的误差可表示
对所选定的on-line(4771.6215 cm-1)和off-line(4772.0240 cm-1),仿真计算5~35 km不同掩星切点高度处水汽干扰产生的相对系统误差(RSE)如
4 探测性能分析
仿真中所采用的激光掩星系统参数和LEO-LEO轨道参数设置分别如
表 1. 系统仿真参数
Table 1. System simulation parameters
|
表 2. LEO-LEO轨道参数
Table 2. LEO-LEO orbital parameters
|
利用卫星工具包STK生成卫星位置坐标,结合(1)~(3)式仿真计算了一次掩星事件中,激光信号扫描高度与时间的关系以及不同高度下的扫描速度,结果如
图 5. 掩星事件中激光信号的扫描特点。(a)切点高度随时间的变化; (b)不同高度下的垂直扫描速度
Fig. 5. Scanning characteristics of laser signals in occultation events. (a) Tangent point altitude as a function of time; (b) vertical scanning speed at different altitudes
由
图 6. 5~35 km不同高度处测量得到的CO2浓度的随机误差
Fig. 6. RRE of CO2 concentration profile at different altitudes from 5 km to 35 km
图 7. 掩星事件全球分布图。(a)一个月; (b)一季度; (c)一年
Fig. 7. Global distributions of occultation events. (a) One month; (b) one quarter; (c) one year
5 结论
红外激光掩星技术是一种重要的全球大气CO2浓度测量手段,具有高浓度测量精度和高全球覆盖率的优点,同时具备垂直分辨能力,可以成为IPDA大气柱浓度测量方法的重要补充。通过对CO2浓度测量随机误差进行理论分析和仿真计算,得到了UTLS区域边界CO2廓线测量随机误差随工作波数的变化,基于UTLS区域整体测量随机误差水平最小和水汽干扰误差最小的原则,选定on-line和off-line波数分别为4771.6215 cm-1和4772.0240 cm-1。探测性能仿真结果显示,所选工作波数对UTLS区域CO2浓度测量的垂直分辨率为0.6~1.4 km,相对测量误差小于0.8%,最高探测精度在10 km高度附近可达到0.229%(对应的体积分数误差为9.16×10-7)。研究结果表明,LIO技术在CO2浓度廓线探测方面具有广阔应用前景,对开展相关设计具有重要的参考价值。
[1] 胡顺星, 陈亚峰, 刘秋武, 等. 差分吸收激光雷达系统探测背景大气SO2和NO2[J]. 中国激光, 2018, 45(9): 0911009.
[2] 谢杨易, 刘继桥, 姜佳欣, 等. 使CO2浓度测量误差减小的星载激光雷达波长优化[J]. 红外与激光工程, 2014, 43(1): 88-93.
Xie Y Y, Liu J Q, Jiang J X, et al. Wavelengths optimization to decrease error for a space-borne lidar measuring CO2 concentration[J]. Infrared and Laser Engineering, 2014, 43(1): 88-93.
[3] 穆永吉, 李蕊, 万渊, 等. 星载激光雷达系统杂散光分析与抑制[J]. 中国激光, 2018, 45(5): 0510005.
[4] 王树志, 朱光武, 白伟华, 等. 风云三号C星全球导航卫星掩星探测仪首次实现北斗掩星探测[J]. 物理学报, 2015, 64(8): 089301.
Wang S Z, Zhu G W, Bai W H, et al. For the first time fengyun3 C satellite-global navigation satellite system occultation sounder achieved spaceborne Bei Dou system radio occultation[J]. Acta Physica Sinica, 2015, 64(8): 089301.
[5] ProschekV, KirchengastG, SchweitzerS, et al. ACCURATE: greenhouse gas profiles retrieval from combined IR-laser and microwave occultation measurements[C]∥EGU General Assembly Conference 2010, May 2-7, 2010, Vienna, Austria. [S.l.: s.n.], 2010: 14352
[6] CliffordD, HoffmannA, WeitnauerC, et al. A novel mission concept for upper air water vapour observations: active limb sounding with a constellation of retroreflectors[C]∥American Geophysical Union, Fall Meeting 2011. [S.l.: s.n.], 2011: A21D- 0115.
[7] Fjeldbo G, Kliore A J, Eshleman V R. The neutral atmosphere of Venus as studied with the mariner V radio occultation experiments[J]. The Astronomical Journal, 1971, 76: 123-140.
[8] Anderson GP, Clough SA, Kneizys FX,et al. AFGL atmospheric constituent profiles (0-120 km)[R]. [S.l.: s.n.], 1986.
[9] Edlén B. The refractive index of air[J]. Metrologia, 1966, 2(2): 71-80.
[10] Sofieva VF, KyröläE. Abel integral inversion in occultation measurements[M] ∥Kirchengast G, Foelsche U, Steiner A K. Occultations for probing atmosphere and climate. Berlin, Heidelberg: Springer, 2004: 77- 85.
Article Outline
李文冬, 刘继桥, 朱亚丹, 陈卫标. LEO-LEO红外激光掩星CO2浓度测量技术研究[J]. 中国激光, 2019, 46(8): 0810001. Wendong Li, Jiqiao Liu, Yadan Zhu, Weibiao Chen. LEO-LEO Infrared Laser Occultation Technique to Measure Atmospheric Carbon Dioxide Concentration[J]. Chinese Journal of Lasers, 2019, 46(8): 0810001.