GNSS掩星反演：从弯曲角到温度廓线的数学之旅
深入推导Abel变换在掩星反演中的应用，探讨球对称假设的局限性、多路径效应的物理机制，以及全谱反演（FSI）和无线电全息方法在低切线高度处的精度改进。
利用GPS、北斗等导航卫星信号穿过地球大气层时的折射效应，反演从地面至60公里高空的温度、湿度、气压与电离层电子密度垂直廓线。全球每日产生超过10,000条高精度大气廓线，被誉为"GPS气象学"的革命性突破。
GNSS无线电掩星（GNSS Radio Occultation, GNSS-RO）是一种利用全球导航卫星系统（GPS、GLONASS、北斗、Galileo）信号穿过地球大气层时发生的折射现象，反演大气折射率、温度、湿度及电离层电子密度垂直廓线的被动遥感技术。它被誉为20世纪末大气探测领域最具革命性的突破之一。
掩星探测的基本物理原理源于大气折射：当导航卫星信号以切线高度穿越地球大气层时，大气密度梯度使信号路径发生弯曲（折射角通常为0.1°–2°）。低轨卫星（LEO）上的高精度GNSS接收机持续跟踪导航卫星的载波相位变化，通过测量信号的多普勒频移偏移量，结合几何光学或波动光学方法，可反演出从地面至约60公里高度的大气折射率廓线。在假设大气处于流体静力平衡状态的条件下，折射率可进一步分解为干空气分量与湿空气分量，从而导出温度与湿度廓线。
与传统探空气球相比，GNSS掩星具有全天候、全球覆盖、高垂直分辨率、无需定标、长期稳定性好等独特优势。截至2025年，全球在轨运行的掩星探测卫星超过20颗，包括美国的COSMIC-2（FORMOSAT-7）6星星座、欧洲MetOp-SG系列、中国风云三号E/F/G星搭载的GNOS-II掩星探测仪，以及商业公司Spire的立方星星座等。这些系统每日联合提供超过10,000条折射率廓线，已成为世界气象组织（WMO）全球观测系统（GOS）和数值天气预报资料同化的核心数据源之一。
在数值天气预报领域，欧洲中期天气预报中心（ECMWF）和美国国家环境预报中心（NCEP）的研究表明，GNSS掩星资料对全球中期预报的相对贡献率约为5–8%，在数据稀疏的南半球海洋和极地地区贡献尤为突出。在气候监测领域，由于掩星探测不依赖辐射定标且不受仪器漂移影响，其长期序列被IPCC第六次评估报告（AR6）列为监测全球平流层温度趋势的关键数据集之一。
从几何光学到全谱反演，解析GNSS信号穿越大气时的物理过程与数据处理方法
当导航卫星（高度约20,200 km）相对于低轨卫星（LEO，高度约500–800 km）逐渐沉入地球大气层后方时，信号切线高度从约100 km持续下降至地面附近。大气折射率n随高度递减导致信号路径弯曲，弯曲角α与大气折射率梯度通过Abel积分变换相关联。切线高度处的折射率n与弯曲角α的关系由Abel变换给出，这是掩星反演的数学基础。
在球对称大气假设下，弯曲角α(a)与折射率n(r)满足Abel积分对：n(r) = exp[(1/π)∫(a到∞) α(a')/√(a'²-a²) da']。实际反演中，首先通过多普勒频移测量计算弯曲角，再进行Abel积分得到折射率廓线。对于低切线高度（<10 km），由于大气水平非均匀性和多路径效应，需引入全谱反演（FSI）或无线电全息方法以提高精度。
大气折射率N = (n-1)×10⁶ = 77.6(P/T) + 3.73×10⁵(e/T²)，其中P为气压（hPa），T为温度（K），e为水汽压（hPa）。在10 km以上干区，可忽略水汽项，直接由折射率计算温度。10 km以下湿区需结合数值预报模式背景场或联合探空资料进行水汽分离。电离层探测则利用双频（L1/L2或B1/B2）信号消除电离层延迟，反演电子密度Ne。
掩星接收机需具备高灵敏度（跟踪信号低至–160 dBm）、高采样率（50–100 Hz）和多系统多频点能力。天线采用全向相位中心稳定设计，以接收从地球边缘升落过程中的导航信号。中国风云三号GNOS-II接收机可同时跟踪GPS L1/L2、北斗B1/B2和Galileo E1/E5a信号，单星每日可产生约500条中性大气廓线和1,000条电离层廓线。
单颗LEO卫星每日仅能产生约100–200次掩星事件，且地理分布受轨道倾角限制。通过多星组网（如COSMIC-2的6颗卫星分布在6个轨道面），可将日产量提升至5,000条以上，并改善时空覆盖均匀性。低倾角轨道（24°）卫星在热带地区掩星密度高，太阳同步轨道卫星则在极区覆盖更优，两者互补形成全球观测网络。
掩星折射率廓线通过变分同化系统（ECMWF 4D-Var、NCEP GSI）直接同化，或转化为弯曲角廓线输入模式。由于掩星资料具有无偏、高垂直分辨率的特点，在模式中的权重系数较高。研究表明，同化掩星资料可使南半球500 hPa高度场预报技巧提升6–10小时，对热带气旋路径预报也有显著正贡献。
基于真实任务数据的可视化分析，从掩星几何到全球观测网络
示意图：GPS/北斗导航卫星（~20,200 km）信号以切线高度穿越大气层，被低轨卫星（LEO，~700 km）接收。随着相对几何关系变化，切线高度从~100 km扫描至地面，信号弯曲角α与大气折射率梯度通过Abel变换关联。蓝色弧线示意信号在大气中的实际弯曲路径，红色虚线为无大气时的直线路径。
数据来源：UCAR COSMIC项目办公室（2024年统计）与中国气象局卫星气象中心。全球主要掩星探测系统每日联合产生约10,000–12,000次中性大气掩星事件。COSMIC-2（6星，24°倾角）在热带地区密度最高；FY-3E/F/G（太阳同步轨道）在高纬度覆盖更优；MetOp-SG提供极轨补充。南半球海洋上空的掩星密度显著高于传统探空站网。
数据来源：基于COSMIC-2实测资料（2024年平均）与美国标准大气1976（US Standard Atmosphere 1976）对比绘制。折射率N（左轴，单位N-unit）在对流层底约300–400 N，随高度指数递减至平流层顶的约1 N。温度廓线（右轴）在对流层顶（约10–12 km）出现极小值后，在平流层下部因臭氧吸收形成显著逆温。掩星探测可分辨100 m尺度的温度垂直结构。
任务参数来源：UCAR、CNSA、EUMETSAT、NOAA公开资料（截至2025年）。图中展示COSMIC-2（6星，720 km，24°倾角）、风云三号E/F/G星（836 km，98.75°倾角，太阳同步）、MetOp-B/C（817 km，98.7°倾角）及Spire商业星座（~500 km，多种倾角）的轨道配置与掩星日产量对比。轨道高度与倾角的组合决定了掩星事件的地理分布特征。
从理论预言到星座组网，三十年间GNSS掩星技术从科学验证走向业务化全球观测
美国科学家Robert A. Phinney和Donald L. Anderson在研究行星无线电掩星观测时，首次系统阐述了利用Abel变换从信号弯曲角反演大气折射率廓线的数学方法。该理论为后来的GNSS掩星探测奠定了坚实的数学物理基础，但当时GPS系统尚未建成，技术条件不成熟。
理论奠基美国加州大学大气遥感合作研究所（UCAR）联合喷气推进实验室（JPL）启动GPS/MET计划，旨在验证利用GPS信号进行大气掩星探测的可行性。MicroLab-1卫星（质量68 kg）开始研制，搭载BlackJack双频GPS接收机，轨道高度735 km，倾角70°。
计划启动MicroLab-1由Pegasus火箭发射入轨。1995年10月，该卫星首次成功反演出大气温度和水汽廓线，验证了GPS掩星探测技术的可行性。GPS/MET在1995–1997年间共获取约40,000次掩星事件，温度反演精度在平流层达到1–2 K，开启了"GPS气象学"（GPS Meteorology）新纪元。
里程碑德国CHAMP卫星（2000年7月发射）搭载BlackJack接收机，首次实现了掩星探测的持续业务化运行，在10年寿命内获取超过600,000条廓线。阿根廷/美国合作的SAC-C卫星（2000年11月发射）同时搭载掩星接收机，标志着掩星探测从科学试验向长期气候监测转型。
业务化萌芽由台湾"国家太空中心"（NSPO）与美国UCAR合作发射的COSMIC（Constellation Observing System for Meteorology, Ionosphere and Climate）6星星座成功入轨。6颗卫星分布在6个轨道面（倾角72°，高度800 km），每日提供约2,000–3,000条中性大气廓线，首次实现了掩星探测的"星座化"，显著改善了全球覆盖的时空分辨率。
星座组网欧洲气象卫星应用组织（EUMETSAT）的MetOp-A极轨气象卫星搭载GRAS（GNSS Receiver for Atmospheric Sounding）接收机投入业务运行。MetOp系列采用太阳同步轨道（高度817 km，倾角98.7°），在极区掩星密度高，与COSMIC形成轨道互补。GRAS采用双天线设计（前向/后向），提高了掩星捕获率。
欧洲业务化中国风云三号C星（FY-3C）成功发射，搭载首台国产GNOS（GNSS Occultation Sounder）掩星探测仪，可同时接收GPS和北斗信号。这是中国首次在气象卫星上开展GNSS掩星探测业务试验，标志着中国进入掩星探测"俱乐部"。FY-3C在轨期间验证了国产掩星载荷的技术可行性。
中国起步COSMIC-2由SpaceX Falcon Heavy火箭一次性发射6颗卫星，轨道高度720 km，倾角24°（低倾角设计），专门优化热带地区掩星密度。搭载TriG接收机（跟踪GPS、GLONASS、Galileo和北斗），每日产生约5,000条中性大气廓线和10,000条电离层廓线，将全球掩星日产量提升至新量级。
精度跃升世界首颗民用晨昏轨道气象卫星FY-3E发射，搭载升级型GNOS-II掩星探测仪。该仪器在GNOS基础上增加了北斗三号全球信号（B1C/B2a）和Galileo E5a跟踪能力，单星日产量提升至约500条中性大气廓线。FY-3E的晨昏轨道（05:30/17:30升交点地方时）填补了全球观测时间窗的空白。
中国升级美国Spire Global公司运营的立方星掩星星座规模扩大至超过150颗卫星，成为目前全球最大的商业掩星数据提供商。Spire卫星采用低轨（~500 km）和多种倾角组合，每日产生超过5,000条廓线，并通过商业合作模式向ECMWF、NOAA等机构提供数据。商业航天的介入大幅降低了掩星探测成本，推动了数据民主化。
商业航天风云三号F星（上午轨道，2023年8月发射）和G星（降水测量，2024年发射）均搭载GNOS-II接收机，与E星形成三星组网。加上COSMIC-2、MetOp-B/C和Spire星座，全球每日掩星廓线产量首次突破10,000条，接近WMO全球探空站网（约1,300站每日2–4次施放）的日观测总量，标志着掩星探测从"补充观测"迈向"核心观测"时代。
万级时代从科学验证卫星到业务化星座，各任务的技术参数与科学产出一览
| 任务名称 | 所属机构 | 发射时间 | 卫星数 | 轨道高度 | 倾角 | 跟踪系统 | 日产量 | 状态 |
|---|---|---|---|---|---|---|---|---|
| GPS/MET | UCAR / NASA | 1995 | 1 | 735 km | 70° | GPS L1/L2 | ~100 | 已结束 |
| CHAMP | GFZ / DLR | 2000 | 1 | 454 km | 87° | GPS L1/L2 | ~200 | 已结束 |
| SAC-C | CONAE / NASA | 2000 | 1 | 705 km | 98° | GPS L1/L2 | ~150 | 已结束 |
| COSMIC (FM-3) | NSPO / UCAR | 2006 | 6 | 800 km | 72° | GPS L1/L2 | ~2,500 | 已结束 |
| MetOp-A/B/C GRAS | EUMETSAT | 2006/2012/2018 | 3 | 817 km | 98.7° | GPS L1/L2 | ~600 | 业务运行 |
| FY-3C GNOS | CMA / NSMC | 2013 | 1 | 836 km | 98.8° | GPS + 北斗 | ~300 | 业务运行 |
| COSMIC-2 (FM-7) | NSPO / NOAA / UCAR | 2019 | 6 | 720 km | 24° | GPS/GLO/GAL/北斗 | ~5,000 | 业务运行 |
| FY-3E/F/G GNOS-II | CMA / NSMC | 2021/2023/2024 | 3 | 836 km | 98.75° | GPS/GAL/北斗三 | ~1,500 | 业务运行 |
| Spire LEMUR | Spire Global | 持续发射 | 150+ | ~500 km | 多种 | GPS/GLO/北斗 | ~5,000+ | 业务运行 |
| MetOp-SG | EUMETSAT | 2025+ | planned 6 | 817 km | 98.7° | GPS/GAL/北斗 | ~3,000 | 计划 |
GNSS掩星探测栏目深度长文，从物理原理到工程实践
深入推导Abel变换在掩星反演中的应用，探讨球对称假设的局限性、多路径效应的物理机制，以及全谱反演（FSI）和无线电全息方法在低切线高度处的精度改进。
对比分析COSMIC-2低倾角星座与FY-3太阳同步轨道的掩星时空采样特征，探讨多星座联合观测对数值预报同化的增益，以及轨道高度对电离层探测精度的影响。
系统评估COSMIC-2和FY-3E掩星廓线在全球不同气候区的精度表现，重点分析低对流层偏差来源、平流层长期稳定性，以及电离层残差对中性大气反演的影响。