星载毫米波遥感作为重要的对地观测手段, 其在自然灾害监测预警、 全球数值天气预报等领域发挥着重要作用^[1]。 相比可见光、 红外遥感, 毫米波遥感具有全天候全天时工作等优势, 其在云雨天气条件下, 可穿透云雨层, 探测大气内部结构变化。 但其同时也有空间分辨率较差、 光谱通道较少、 光谱分辨率偏低等缺点。 目前毫米波遥感器只部署在极地轨道卫星上, 静止轨道卫星尚未有毫米波遥感器搭载^[2]。 相比静止轨道卫星, 极地轨道卫星的重返周期要长很多。 极轨卫星一天两次的过境观测频率不能满足对于需要高重复周期监测的目标及自然灾害等观测需求。

目前, 国内外极地轨道卫星上已大量搭载毫米波遥感器^{[2, 3]}。 我国自主发展的长观测序列毫米波遥感器, 主要从风云三号极轨气象卫星开始。 在风云三号A星(FY3A, FengYun3A)和FY3B上搭载有我国第一代业务化的毫米波辐射探测器MWHS, 其工作频率设置在183.31 GHz的水汽吸收通道, 用于探测全球大气水汽的分布信息。 同时我国也已规划在后续风云静止轨道卫星上部署静止轨道毫米波遥感器^[4]。 由于其轨道高度较高, 因此, 在传统极地轨道卫星所选毫米波探测通道基础上, 静止轨道毫米波探测拟加入更高频段的探测通道, 通过缩短遥感器工作波长, 改善遥感器空间分辨率。 目前拟选择425 GHz的氧气吸收通道, 主要用于大气中温度信息的探测。 相关研制单位已开展地面样机试制工作, 试制样机同时包含183和424 GHz系列毫米波通道^[5]。

毫米波辐射传输建模是对毫米波遥感观测资料开展相关科学研究的基础。 目前, 毫米波辐射传输模型主要有多层毫米波辐射传输模型和快速毫米波辐射传输模型两大类。 前者用于分析计算, 后者通常用于数值天气预报模式等应用。 两类毫米波辐射传输模型中, 均使用毫米波大气吸收模型来计算毫米波辐射传输模型中晴空大气对毫米波的吸收及晴空大气自身毫米波发射^[6]。 因此, 毫米波大气吸收模型的精确性对毫米波辐射传输模型的精度具有决定性作用。 毫米波大气吸收模型传统上以逐线吸收模型为主, 常见的主要有MPM模型^[7]、 ROS模型^[8]、 HITRAN模型^[9]等。 逐线吸收模型均基于吸收线库进行模拟计算。

目前由于静止轨道毫米波遥感仍在前期探索中, 对静止轨道所可能使用的毫米波频段探测通道性能的具体分析较少, 对几种常见逐线吸收模型的模拟性能分析比较也较少。 本文利用MPM, ROS和HITRAN三种逐线大气吸收模型, 分析不同吸收模型在拟新增的静止轨道424 GHz工作频段处, 大气吸收模拟特性的差异。 同时, 还利用三种模型对海陆交界地区的183.31 GHz毫米波辐射模拟性能进行分析。

根据静止轨道毫米波地面样机试制所选用的183.31和424 GHz通道, 本研究拟分别对183.31和424 GHz通道毫米波辐射特征进行分析。

针对424 GHz通道, 由于目前没有相近通道的在轨实际观测进行数据验证, 因而本研究对424 GHz通道在不同温度、 气压、 湿度条件下的毫米波大气吸收特征进行分析比较。

针对183.31 GHz通道, 由于FY3B/MWHS上搭载有类似观测通道, 且其在轨运行时间较长, 积累较多观测数据, 因而选择其183.31 GHz所在三通道(通道参数如表1所示)数据, 验证逐线毫米波大气吸收模型对毫米波辐射模拟的影响。 (卫星数据来源为中国气象局国家卫星气象中心数据下载网站)183.31 GHz验证工作需要FY3B卫星观测同时, 地面也有探空仪进行大气参数观测, 利用此大气参数输入进毫米波辐射传输模型, 模拟得到毫米波辐射值, 该值与卫星观测差异可反映模拟模型的性能。

表1

Table 1

表1(Table 1)

表1 FY3B MWHS通道参数 Table 1 Channel parameters of FY3B MWHS 通道序号 中心频率/GHz 带宽/MHz 极化方式 1 150 1 000 V 2 150 1 000 H 3 183.31± 7.0 2 000 H 4 183.31± 3.0 1 000 H 5 183.31± 1.0 500 H

表1 FY3B MWHS通道参数 Table 1 Channel parameters of FY3B MWHS

海陆交界地区同时兼具陆地和海洋气候特征, 其大气波导特征参数, 具有一定的代表性。 上海探空站具有下午14点加密观测能力, FY3B卫星的过境时间通常为当地时间13:30— 14:30, 因而选择上海探空站观测可有较多的加密观测资料与FY3B卫星观测形成时空匹配。 (探空观测数据来源为中国气象局大气探测中心数据下载网站)

毫米波逐线大气吸收模型用于模拟大气中各气态分子对毫米波辐射的吸收作用。 通过对大气中各主要气态成分分子对毫米波的共振与非共振吸收计算, 综合得到整个大气对毫米波辐射的吸收作用。

对于一束频率为f平面波在大气中传播了距离z后, 其电场相位和幅度的综合响应可由式(1)表示

$ E(z)=\exp \left[i k z\left(1+N \times 10^{-6}\right)\right\rfloor E(0)$(1)

其中, E(0)为初始电场, k=2π f/c为自由空间中的电磁波传播波数, c为电磁波在真空中的传播速度。 大气媒质的光谱特性由复折射率式(2)表示

N=N0+N'+iN″(2)

式(2)中, N₀是一个与工作频率无关的分量, N'与N″均与电磁波的工作频率有关的分量, 前者表示电磁波的传播延时性质, 后者为N的虚部, 表示了其对于电磁波能量的吸收特性。 大气吸收系数由式(3)可计算

α=0.1820fN″(3)

大气中的干空气折射率由式(4)计算得到,

ND=Nd+∑kSkFk+Nn(4)

其中, k为氧气吸收谱线数量, N_d为非色散分量, N_n为非共振分量, 式(4)右端中间累加项为依据微波光谱吸收线库所计算得到的干空气共振折射分量。 其中S_k为氧气谱线强度, F_k为氧气谱线线型函数。 大气中水汽引起的复折射率可由式(5)计算,

NV=Nv+∑tStFt+Nc(5)

其中, t为水汽吸收谱线数量, N_v为非色散分量, N_c为水汽贡献非共振分量, 式(5)右端中间累加项为依据微波光谱吸收线库所计算得到的水汽共振折射分量。 其中, S_t为水汽谱线强度, F_t为水汽谱线线型函数。

对于不同逐线吸收模型, 在选择吸收线型库, 计算共振折射分量时不同, 最终计算出的大气波导中毫米波的吸收效应也不同。

MPM模型选用44线的氧气吸收线库, 和35线的水汽吸收线库, ROS模型选用45线的氧气吸收线库和15线的水汽吸收线库, HITRAN模型则分别选择氧气和水汽吸收线库中吸收线强度前50位的吸收线组成吸收线库, 进行大气吸收系数计算。 (具体计算过程见参考文献[7])

影响毫米波大气辐射特征的主要因素有大气温度、 大气压强、 大气湿度, 通过固定其中两个参数, 分析毫米波大气吸收系数随另一参数的变化而变化。

2.1 毫米波大气吸收的温度敏感性分析

三种逐线吸收模型对大气温度的敏感性通过固定其他大气参数, 模拟大气吸收系数随大气温度变化而得到。 设置模拟工作频率为424 GHz, 大气压强设定为1 013.25 hPa, 大气湿度设定为0 g· m^-3, 假设大气物理温度由200 K逐渐增大至310 K。 针对毫米波探测通道, 三种逐线吸收模型模拟出的大气吸收系数随大气温度变化如图1所示。 对于424 GHz通道, MPM与ROS模型较为接近, HITRAN模型在温度较低(即海拔高度较高)时, 明显偏低, 随着大气温度增大, HITRAN模型模拟值接近MPM模型, 而ROS模型值略高于其他两种模型。

图1

Fig.1

	Figure Option ViewDownloadNew Window
	图1 大气温度对424 GHz通道大气吸收系数模拟影响Fig.1 Impacts of atmospheric temperature on atmospheric absorption coefficient simulation of Channel 424 GHz

2.2 毫米波大气吸收的气压敏感性分析

三种逐线吸收模型对大气压强的敏感性通过固定其他大气参数, 模拟大气吸收系数随大气压强变化而得到。 假设模拟工作频率为424 GHz, 大气物理温度设定为273.15 K, 大气湿度设定为0 g· m^-3, 假设大气压强由1 hPa逐渐增大至1 050 hPa。 针对毫米波探测通道, 三种逐线吸收模型计算出的大气吸收系数随气压变化如图2所示。 对于424 GHz通道, MPM与ROS模型在气压较低时模拟结果较为接近, 随着气压增大二者差异增大, MPM模型逐渐小于ROS模型, HITRAN模型在整个模拟的气压范围内均明显低于另两种模型, 随着气压增大, 其模拟结果接近MPM模型。

图2

Fig.2

	Figure Option ViewDownloadNew Window
	图2 大气压强对424 GHz通道大气吸收系数模拟影响Fig.2 Impacts of atmospheric pressure on atmospheric absorption coefficient simulation of Channel 424 GHz

2.3 毫米波大气吸收的湿度敏感性分析

三种逐线吸收模型对大气湿度的敏感性通过固定其他大气参数, 模拟大气吸收系数随大气湿度变化而得到。 假设模拟工作频率为424 GHz, 大气物理温度设定为273.15 K, 大气压强设定为1 013.15 hPa, 假设大气湿度由0 g· m^-3逐渐增大至1 g· m^-3。 针对毫米波探测通道, 三种逐线吸收模型模拟出的大气吸收系数随大气湿度变换如图3所示。 对于424 GHz通道, 三种模型在湿度较低时, 差异明显, ROS模型模拟值最大, MPM模型次之, HITRAN模型最小, 随着大气湿度增大, MPM模型与ROS模型间的差异减小, 两者模拟结果趋于一致, 而HITRAN模型随着湿度增大, 与另两种模型间的差异逐渐增大, 模拟值明显低于MPM模型与ROS模型。

图3

Fig.3

	Figure Option ViewDownloadNew Window
	图3 大气湿度对424 GHz通道大气吸收系数模拟影响Fig.3 Impacts of atmospheric humidity on atmospheric absorption coefficient simulation of Channel 424 GHz

综上三种情况分析, 在424 GHz频段, MPM模型与ROS模型计算的毫米波大气吸收系数随大气温度、 大气压强和大气湿度的变化特征更相似, 而HITRAN模型的数值则普遍小于MPM及ROS模型, 但三个模型随大气温度、 大气压强和大气湿度变化的趋势则相同。

静止轨道毫米波观测是未来对地观测发展的重要方向。 笔者引入三种逐线吸收模型, 评估未来拟在静止轨道应用的424 GHz毫米波探测通道对大气温度、 湿度、 气压的敏感性, 并分析不同逐线吸收模型的变化特征。 同时, 引入多层毫米波辐射传输模型, 结合FY3B/MWHS真实观测, 评估在多层辐射传输模型中, 三种逐线吸收模型模拟毫米波辐射的性能。 MPM模型的模拟效果最优, 对FY3B/MWHS的三个毫米波通道均为如此, HITRAN模型的模拟误差最大, 显示MPM模型采用的光谱吸收线库及计算精度最高。 需要指出的是, 183.31 GHz通道为水汽吸收通道, 前述模拟分析结果更多显示的是水汽吸收的毫米波通道特征。 同时, 上海地区水汽含量中等, 对于高湿区或干旱区的情况有待进一步讨论, 且探空观测的大气垂直信息本身及进行的插值运算也带入误差。 另外, 对于毫米波氧气吸收通道的特征分析还有待于后续进一步开展。