Page 250 - 《高原气象》2025年第6期
P. 250
高 原 气 象 44 卷
1658
线性插值方法插值至掩星观测的中心经纬度位置, θ 的关系有 Snell 定律: r·n·sinθ = A (常数), 通常用
然后将掩星资料插值到 ERA5再分析资料的垂直气 Abel 积分逆变换把弯曲角序列反演得到大气折射
压高度上, 以保持水平和垂直方向的一致性(孟恬 指数剖面:
等, 2021; 殷延安和杨胜朋, 2023)。这种严格的时 é ê ê 1 ∞ α (a ) ù ú ú ú
ê
ê
空匹配方法确保了所有对比分析基于相同的时空 n ( x ) = exp ê ê ∫ x a - x 2 da ú ú ú û (6)
π
ë
2
条件, 从而提高了评估结果的准确性和可信度。 a
6
2. 2 观测算子 大气折射率为N = (n - 1) × 10 ; 海拔h = n (a ) -
2. 2. 1 折射率局地观测算子 R, 其中 R 是切点处的曲率半径, a 是前面定义的碰
中性大气中, 大气折射率是水汽、 气压和温度 撞参数。
的函数, 模式背景场往往不包含云中的液态水项和 文中所有图均是基于国家测绘地理信息局标
冰水项, 基于大气压强计算折射率N如下: 准地图服务网站下载的审图号为 GS(2016)1611 号
P
5
N = 77.6 ⋅ + 3.73 × 10 ⋅ P w (1) 的世界地图制作, 底图无修改。
T T 2
式中: P 为大气压强(单位: hPa); P 为水汽压(单 3 结果分析
W
位: hPa); T为温度(单位: K)。
3. 1 最低下探高度
ERA5 不直接提供水汽压及折射率信息, 需要
本研究利用 2024年 5月 6 -9日期间, 全球范围
基于相对湿度及温度格网资料逐步求解, 采用
内的 FY-3E/GNOS 掩星与 COSMIC-2 掩星资料进行
Goff-Gratch 方程计算水汽压。气象格网资料的高
了对比分析, 并将两者的折射率和弯角廓线与对应
度和平面上分别采用三次样条和双线性插值计算
的 ERA5的再分析数据进行对比。根据 2023年 6月
掩星廓线对应的大气廓线资料。
12 日 FY-3E 卫星掩星附加相位统计沿纬度分布的
在对比分析中, 折射率和弯角采用相对偏差
最低下探深度如图1所示。
RE(Relative Error)和相对偏差标准差 STD(Stan‐
图 1为 2023年 6月 12日 FY-3E卫星的最低下探
dard Deviation of Relative Error)来验证精度, 分别
高度分布。结果表明, FY-3E 卫星在全球范围内具
由公式(2)、 公式(4)得到。以折射率为例:
N - N 有较强的低层大气探测能力, 大部分廓线最低的探
RE i = o,i b,i × 100% (2) 测深度都在 0. 5 km 以下。赤道和热带地区的下探
N b,i
1 N N o,i - N b,i 深度相较于高纬地区明显偏高, 这与该区域内丰富
MRE = ∑ N b,i × 100% (3) 的水汽含量以及复杂的气象条件密切相关。水汽
N
i = 1
1 i = 1( N - N ) 2 在掩星信号传播过程中会导致无线电信号的散射
N
STD = ∑ o,i b,i × 100% - MRE (4) 和吸收, 从而产生模糊散焦效应。此外, 超折射现
N - 1 N
b,i
象使得大气层中温度和湿度的强梯度会导致无线
2. 2. 2 弯角观测算子 电波的折射角度异常增大, 进而影响信号的传输路
当 GNSS 信号穿过大气层时, 由于介质密度不 径, 使得实际下探深度受到限制。
均匀, 信号传播路径会发生一定的偏折而产生弯曲 3. 2 FY-3E与COSMIC-2掩星资料廓线对比分析
角 α, 根据通过几何光学原理从相位延迟序列中获 通过时空匹配, 2024 年 5 月 6 -9 日期间共获得
得这种弯曲角, 可基于 Snell 定律和 Abel 积分公式 474 条有效廓线, 为清晰展示数据特征, 特选取具
反演大气折射率剖面(Zou and Zeng, 2006)。在局 有代表性的单根廓线进行可视化。从图 2中可以看
部球对称的假设下, 该弯曲角可与碰撞参数 a(定义 出, FY-3E 和 COSMIC-2 的对比揭示了二者在 10~
为折射中心点到信号射线的垂直距离)建立起唯一 40 km 高度范围内的整体一致性, 体现了 FY-3E/
的对应关系。在去除电离层效应并统计优化后, 使 GNOS 掩星数据和 COSMIC-2 的一致性。值得注意
用Abel积分变化将弯曲角转换为大气折射率。 的是, FY-3E/GNOS 掩星数据资料的弯角和折射率
∞ 1 dn adx
α(a) = -2a ∫ n dx ⋅ (5) 廓线呈现明显的平滑特征, 这表明 FY-3E 掩星数据
2
a x - a 2 的处理过程中, 可能使用了较粗的采样频率或者较
式中: n为折射率; a为信号路径近地点离地心的距 强的滤波算法以减少观测数据中的噪声, 提升信号
离, 也称碰撞参数。在地球大气介质局部球对称的 清晰度。然而, 这一滤波过程可能导致了过度平
假设下, 球面大气分层的半径 r与折射率 n、 折射角 滑, 进而模糊或丢失了某些关键的大气层结构细节。

