留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

FY-3C无线电掩星折射率廓线的反演及验证

徐晓华 张纪满 罗佳 高攀

徐晓华, 张纪满, 罗佳, 高攀. FY-3C无线电掩星折射率廓线的反演及验证[J]. 武汉大学学报 ● 信息科学版, 2022, 47(1): 36-44. doi: 10.13203/j.whugis20190361
引用本文: 徐晓华, 张纪满, 罗佳, 高攀. FY-3C无线电掩星折射率廓线的反演及验证[J]. 武汉大学学报 ● 信息科学版, 2022, 47(1): 36-44. doi: 10.13203/j.whugis20190361
XU Xiaohua, ZHANG Jiman, LUO Jia, GAO Pan. Retrieval and Validation of FY-3C Radio Occultation Refractivities[J]. Geomatics and Information Science of Wuhan University, 2022, 47(1): 36-44. doi: 10.13203/j.whugis20190361
Citation: XU Xiaohua, ZHANG Jiman, LUO Jia, GAO Pan. Retrieval and Validation of FY-3C Radio Occultation Refractivities[J]. Geomatics and Information Science of Wuhan University, 2022, 47(1): 36-44. doi: 10.13203/j.whugis20190361

FY-3C无线电掩星折射率廓线的反演及验证

doi: 10.13203/j.whugis20190361
基金项目: 

国家自然科学基金 41774033

国家自然科学基金 41774032

国家重点研发计划 2018YFC1503502

武汉大学地球空间环境与大地测量教育部重点实验室开放基金 18-02-10

重庆市基础研究与前沿探索专项 cstc2018jcyjAX0015

详细信息
    作者简介:

    徐晓华,博士,教授,主要从事GNSS大气探测及应用的研究。xhxu@sgg.whu.edu.cn

    通讯作者: 罗佳,博士,教授。jluo@sgg.whu.edu.cn
  • 中图分类号: P228

Retrieval and Validation of FY-3C Radio Occultation Refractivities

Funds: 

The National Natural Science Foundation of China 41774033

The National Natural Science Foundation of China 41774032

the National Key Research and Development Program of China 2018YFC1503502

the Open Project Fund of the Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, Wuhan University 18-02-10

the Chongqing Research Program of Basic Research and Frontier Technology cstc2018jcyjAX0015

More Information
    Author Bio:

    XU Xiaohua, PhD, professor, specializes in GNSS atmospheric and application. E-mail: xhxu@sgg.whu.edu.cn

    Corresponding author: LUO Jia, PhD, professor. E-mail: jluo@sgg.whu.edu.cn
  • 摘要: 2013年, 中国发射了首颗进行全球导航卫星系统(global navigation satellite system,GNSS)无线电掩星观测的气象卫星FY-3C,其掩星数据产品已由国家卫星气象中心(national satellite meteorological center,NSMC)发布。基于FY-3C附加相位数据,通过几何光学反演和全谱反演相结合的方法自主反演获取大气折射率廓线,并利用NSMC提供的折射率产品和欧洲中期天气预报中心(European Centre for Medium-Range Weather Forecasts, ECMWF)的气候再分析数据集(ECMWF reanalysis-interim, ERA-Interim)对反演结果进行验证。2017年全年FY-3C掩星折射率自主反演的成功率为82%,与NSMC 83.5%的反演成功率基本一致。以ERA-Interim再分析资料为参考,自主反演折射率相对偏差绝对值在25 km以下, 小于1.5%,略优于NSMC产品,但偏差标准差比NSMC产品略大,达到了3.6%。自主反演和NSMC的折射率廓线在高纬地区的相对偏差均低于中低纬地区,且自主反演折射率的相对偏差总体上在夏季高于冬季。
  • 图  1  FY-3C掩星折射率反演流程图

    Figure  1.  Flowchart for Retrieving the FY-3C Radio Occultation Refractivity Profiles

    图  2  FY-3C折射率廓线反演及验证个例

    Figure  2.  Example for the Retrieval and Validation of the Refractivity Profile from FY-3C Radio Occultation Event

    图  3  2017年折射率廓线数对比

    Figure  3.  Comparison of the Number of the Refractivity Profiles in 2017

    图  4  NSMC和自主反演折射率廓线统计验证

    Figure  4.  Statistical Validations of the Refractivity Profiles Provided by NSMC and Those Retrieved Independently

    图  5  FY-3C掩星NSMC折射率产品不同纬度带的统计比较

    Figure  5.  Statistical Comparison of NSMC FY-3C Refractivity Products over Different Latitudes

    图  6  FY-3C掩星自主反演折射率廓线不同纬度带的统计比较

    Figure  6.  Statistical Comparison of FY-3C Refractivity Profiles Retrieved Independently over Different Latitudes

    图  7  NSMC的FY-3C掩星折射率产品相对偏差不同季节的纬度-高度分布

    Figure  7.  Latitude-Height Distribution of the Relative Biases of FY-3C Radio Occultation Refractivity Products Provided by NSMC in Different Seasons

    图  8  自主反演的FY-3C掩星折射率廓线相对偏差不同季节的纬度-高度分布

    Figure  8.  Latitude-Height Distribution of the Relative Biases of FY-3C Radio Occultation Refractivity Products Retrieved Independently in Different Seasons

  • [1] Fischbach F F. A Satellite Method for Pressure and Temperature Below 24 km[J]. Bulletin of the American Meteorological Society, 1965, 46(9): 528-532 doi:  10.1175/1520-0477-46.9.528
    [2] Hardy K R, Hajj G A, Kursinski E R. Accuracies of Atmospheric Profiles Obtained from GPS Occultations[J]. International Journal of Satellite Communications, 1994, 12(5): 463-473 doi:  10.1002/sat.4600120508
    [3] Syndergaard S. Modeling the Impact of the Earth's Oblateness on the Retrieval of Temperature and Pressure Profiles from Limb Sounding[J]. Journal of Atmospheric and Solar-Terrestrial Physics, 1998, 60(2): 171-180 doi:  10.1016/S1364-6826(97)00056-4
    [4] 廖蜜, 张鹏, 毕研盟, 等. 风云三号气象卫星掩星大气产品精度的初步检验[J]. 气象学报, 2015, 73(6): 1131-1140 https://www.cnki.com.cn/Article/CJFDTOTAL-QXXB201506011.htm

    Liao Mi, Zhang Peng, Bi Yanmeng, et al. A Preliminary Estimation of the Radio Occultation Products Accuracy from the Fengyun-3C Meteorologica Satellite[J]. Acta Meteorologica Sinica, 2015, 73(6): 1131-1140 https://www.cnki.com.cn/Article/CJFDTOTAL-QXXB201506011.htm
    [5] Liao M, Zhang P, Yang G L, et al. Preliminary Validation of the Refractivity from the New Radio Occultation Sounder GNOS/FY-3C[J]. Atmospheric Measurement Techniques, 2016, 9(2): 781-792 doi:  10.5194/amt-9-781-2016
    [6] 徐晓华, 朱洲宗, 罗佳. 基于IGRA2探空数据和COSMIC掩星资料对FY-3C掩星中性大气产品的质量分析[J]. 武汉大学学报·信息科学版, 2020, 45(3): 384-393 https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH202003009.htm

    Xu Xiaohua, Zhu Zhouzong, Luo Jia. Quality Analysis of the Neutral Atmospheric Products from FY-3C Radio Occultation Based on IGRA2 Radiosonde Data and COSMIC RO Products[J]. Geomatics and Information Science of Wuhan University, 2020, 45 (3): 384-393 https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH202003009.htm
    [7] Dee D P, Uppala S M, Simmons A J, et al. The ERA-Interim Reanalysis: Configuration and Performance of the Data Assimilation System[J]. Quarterly Journal of the Royal Meteorological Society, 2011, 137(656): 553-597
    [8] Lohmann M S. Analysis of Global Positioning System(GPS)Radio Occultation Measurement Errors Based on Satellite de Aplicaciones Cientificas-C (SAC-C)GPS Radio Occultation Data Recorded in Open-Loop and Phase-Locked-Loop Mode[J]. Journal of Geophysical Research: Atmospheres, 2007, 112(D9), DOI:  10.1029/2006JD007764
    [9] Jensen A S, Lohmann M S, Benzon H H, et al. Full Spectrum Inversion of Radio Occultation Signals [J]. Radio Science, 2003, 38(3): 1040-1055
    [10] 徐晓华, Zhang Kefei, Fu Erjiang, 等. 澳大利亚上空COSMIC掩星廓线的反演[J]. 武汉大学学报·信息科学版, 2008, 33(8): 800-804 http://ch.whu.edu.cn/article/id/1685

    Xu Xiaohua, Zhang Kefei, Fu Erjiang, et al. Inversion of Atmospheric Profiles with COSMIC Radio Occultation Data over Australia[J]. Geomatics and Information Science of Wuhan University, 2008, 33 (8): 800-804 http://ch.whu.edu.cn/article/id/1685
  • [1] 王耀兴, 张秋昭, 沈震.  COSMIC掩星反演湿温廓线质量的时空分析 . 武汉大学学报 ● 信息科学版, 2021, 46(6): 887-894. doi: 10.13203/j.whugis20190051
    [2] 朱洲宗, 徐晓华, 罗佳.  利用FY-3C折射率对大气边界层高度的反演与分析 . 武汉大学学报 ● 信息科学版, 2021, 46(3): 395-401. doi: 10.13203/j.whugis20190271
    [3] 张绍成, 黄龙强, 殷飞, 侯威震.  星基反演TEC与地基GNSS TEC数据融合可行性分析 . 武汉大学学报 ● 信息科学版, 2020, 45(4): 557-564. doi: 10.13203/j.whugis20180293
    [4] 徐晓华, 朱洲宗, 罗佳.  利用IGRA2探空数据和COSMIC掩星资料对FY-3C掩星中性大气产品进行质量分析 . 武汉大学学报 ● 信息科学版, 2020, 45(3): 384-393. doi: 10.13203/j.whugis20180490
    [5] 徐晓华, 刘树纶, 罗佳.  利用COSMIC掩星折射指数分析全球大气边界层顶结构变化 . 武汉大学学报 ● 信息科学版, 2018, 43(1): 94-100. doi: 10.13203/j.whugis20160183
    [6] 安家春, 杜玉军, 屈小川, 杨剑.  掩星大气反演中的电离层二阶项效应的分析 . 武汉大学学报 ● 信息科学版, 2015, 40(11): 1440-1445. doi: 10.13203/j.whugis20150228
    [7] 江 鹏, 叶世榕, 何书镜, 刘炎炎.  自适应Kalman滤波用于GPS层析大气湿折射率 . 武汉大学学报 ● 信息科学版, 2013, 38(3): 299-302.
    [8] 自适应Kalman滤波用于GPS层析大气湿折射率 . 武汉大学学报 ● 信息科学版, 2013, 38(3): 299-.
    [9] 徐晓华, 汪海洪.  不同季节GPS掩星廓线精度的比较研究 . 武汉大学学报 ● 信息科学版, 2010, 35(6): 639-643.
    [10] 刘经南, 赵莹, 张小红.  GNSS无线电掩星电离层反演技术现状与展望 . 武汉大学学报 ● 信息科学版, 2010, 35(6): 631-635.
    [11] 赵莹, 张小红.  COSMIC掩星观测数据反演电离层电子密度廓线 . 武汉大学学报 ● 信息科学版, 2010, 35(6): 644-648.
    [12] 徐晓华, 罗佳.  COSMIC掩星折射指数廓线的统计验证 . 武汉大学学报 ● 信息科学版, 2009, 34(2): 214-217.
    [13] 张双成, 叶世榕, 万蓉, 陈波.  基于Kalman滤波的断层扫描初步层析水汽湿折射率分布 . 武汉大学学报 ● 信息科学版, 2008, 33(8): 796-799.
    [14] 范磊, 符养, 杜晓勇, 李黄.  雾灵山山基掩星观测反演误差分析 . 武汉大学学报 ● 信息科学版, 2008, 33(1): 89-92.
    [15] 徐晓华, ZhangKefei, FuErjiang, 罗佳.  澳大利亚上空COSMIC掩星廓线的反演 . 武汉大学学报 ● 信息科学版, 2008, 33(8): 800-804.
    [16] 徐晓华, 李征航, 罗佳.  GPS掩星数据反演中的地球扁率影响改正 . 武汉大学学报 ● 信息科学版, 2005, 30(6): 502-505.
    [17] 徐晓华, 李征航, 罗佳.  利用GPS/MET与CHAMP掩星资料反演地球中性大气参数的比较 . 武汉大学学报 ● 信息科学版, 2004, 29(3): 214-217.
    [18] 徐晓华, 李征航, 罗佳.  利用GPS掩星资料反演地球中性大气参数折射角方法研究 . 武汉大学学报 ● 信息科学版, 2003, 28(5): 589-592.
    [19] 余模智, 林银森, 仲思东.  透镜折射率的高精度非接触测量方法 . 武汉大学学报 ● 信息科学版, 1994, 19(2): 173-177.
    [20] 余模智.  玻璃微珠折射率测定 . 武汉大学学报 ● 信息科学版, 1994, 19(4): 371-373.
  • 加载中
图(8)
计量
  • 文章访问数:  248
  • HTML全文浏览量:  82
  • PDF下载量:  35
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-12-10
  • 刊出日期:  2022-01-05

FY-3C无线电掩星折射率廓线的反演及验证

doi: 10.13203/j.whugis20190361
    基金项目:

    国家自然科学基金 41774033

    国家自然科学基金 41774032

    国家重点研发计划 2018YFC1503502

    武汉大学地球空间环境与大地测量教育部重点实验室开放基金 18-02-10

    重庆市基础研究与前沿探索专项 cstc2018jcyjAX0015

    作者简介:

    徐晓华,博士,教授,主要从事GNSS大气探测及应用的研究。xhxu@sgg.whu.edu.cn

    通讯作者: 罗佳,博士,教授。jluo@sgg.whu.edu.cn
  • 中图分类号: P228

摘要: 2013年, 中国发射了首颗进行全球导航卫星系统(global navigation satellite system,GNSS)无线电掩星观测的气象卫星FY-3C,其掩星数据产品已由国家卫星气象中心(national satellite meteorological center,NSMC)发布。基于FY-3C附加相位数据,通过几何光学反演和全谱反演相结合的方法自主反演获取大气折射率廓线,并利用NSMC提供的折射率产品和欧洲中期天气预报中心(European Centre for Medium-Range Weather Forecasts, ECMWF)的气候再分析数据集(ECMWF reanalysis-interim, ERA-Interim)对反演结果进行验证。2017年全年FY-3C掩星折射率自主反演的成功率为82%,与NSMC 83.5%的反演成功率基本一致。以ERA-Interim再分析资料为参考,自主反演折射率相对偏差绝对值在25 km以下, 小于1.5%,略优于NSMC产品,但偏差标准差比NSMC产品略大,达到了3.6%。自主反演和NSMC的折射率廓线在高纬地区的相对偏差均低于中低纬地区,且自主反演折射率的相对偏差总体上在夏季高于冬季。

English Abstract

徐晓华, 张纪满, 罗佳, 高攀. FY-3C无线电掩星折射率廓线的反演及验证[J]. 武汉大学学报 ● 信息科学版, 2022, 47(1): 36-44. doi: 10.13203/j.whugis20190361
引用本文: 徐晓华, 张纪满, 罗佳, 高攀. FY-3C无线电掩星折射率廓线的反演及验证[J]. 武汉大学学报 ● 信息科学版, 2022, 47(1): 36-44. doi: 10.13203/j.whugis20190361
XU Xiaohua, ZHANG Jiman, LUO Jia, GAO Pan. Retrieval and Validation of FY-3C Radio Occultation Refractivities[J]. Geomatics and Information Science of Wuhan University, 2022, 47(1): 36-44. doi: 10.13203/j.whugis20190361
Citation: XU Xiaohua, ZHANG Jiman, LUO Jia, GAO Pan. Retrieval and Validation of FY-3C Radio Occultation Refractivities[J]. Geomatics and Information Science of Wuhan University, 2022, 47(1): 36-44. doi: 10.13203/j.whugis20190361
  • 无线电掩星技术最早被应用于火星和太阳系其他行星的大气探测。20世纪60年代,Fischbach[1]提出了利用掩星技术进行地球大气探测的思想;Hardy等[2]进行了理论研究及模拟计算,但由于地球周围没有符合条件的观测星体,该计划一直搁浅。20世纪80年代,随着全球定位系统(global positioning system,GPS)的建立及发展,美国大学大气研究联盟主持实施了GPS气象(GPS meteorology,GPS/MET)掩星探测实验任务,此后,各国和地区相继实施了以气象、电离层与气候星座观测系统(constellation observing system for meteorology ionosphere and climate,COSMIC)为代表的多项掩星任务。掩星提供的大气参数廓线具有高精度、高垂直分辨率、高稳定性、全天候、全球分布等优点,目前,已被各大数值天气预报中心同化到数值天气预报模式中,对天气预报和气候研究产生了积极的影响[3]。相比之下,中国无线电掩星探测计划起步较晚,2013年9月,国家气象局主持发射了FY-3C极轨气象卫星,该低轨卫星(low earth orbit,LEO)搭载的全球导航卫星掩星探测仪(global navigation satellite system occultation sounder,GNOS)除了对GPS卫星进行掩星观测,也兼容北斗导航卫星系统(BeiDou navigation satellite system,BDS)。该掩星探测仪每天可观测的GPS掩星事件为500次左右,GPS与北斗掩星事件共计1 000次左右[4]。作为中国首颗进行全球导航卫星系统(global navigation satellite system,GNSS)掩星观测的气象卫星,FY-3C卫星的掩星数据产品对于中国后续掩星计划实施具有重要的参考价值。

    虽然国家卫星气象中心(national satellite meteorological center,NSMC)发布了自2014年以来的FY-3C掩星数据产品,但国内外关于FY-3C掩星中性大气资料处理及验证的文章相对较少。廖蜜等[4]对NSMC提供的2013-10—2013-11的产品的质量检验结果表明,FY-3C掩星中性大气产品的精度基本达到预定目标,但存在一定改进空间;廖蜜等[5]进一步介绍了NSMC进行FY-3C掩星中性大气反演的主要流程,其中,低对流层区域的反演采用了第2类正则变换方法,并以欧洲中期天气预报中心(European Centre for Medium-Range Weather Forecasts,ECMWF)的气候再分析数据集(ECMWF reanalysis-interim,ERA-Interim)为参考,对NSMC 2013-10—2013-11的FY-3C GPS掩星及BDS掩星测试数据进行了检核。以上工作中选用时间段较短,且该时间段数据未正式公布。徐晓华等[6]以IGRA2探空资料和COSMIC掩星资料为参照,对2015—2017年的FY-3C掩星折射率、温度及比湿等中性大气产品的精度进行分析,并进一步对精度随季节、纬度带、昼夜和海陆分布的变化进行了讨论。但该工作中仅对NSMC公布的产品进行了核验,并未涉及基于附加相位数据的自主反演。

    本文基于NSMC提供的FY-3C掩星附加相位数据,实现了大气折射率廓线的自主反演。与NSMC所采用的反演方法不同,本文的自主反演过程中,通过结合全谱反演(full spectrum inversion,FSI)和几何光学反演方法获取弯曲角廓线,以ERA-Interim再分析资料为参考,分别对自主反演和NSMC公布的折射率产品进行了统计验证,并对折射率廓线在不同纬度带的误差特性及不同季节的相对偏差特征进行了分析。

    • NSMC公布了FY-3C的各级掩星数据产品。其中,1级产品为精密轨道数据,2级产品为中性大气及电离层的附加相位数据,3级产品包含大气折射率、温度、密度、湿度等大气参数廓线和电子密度廓线。本文基于附加相位数据进行了大气弯曲角和折射率廓线的反演。需要说明的是,目前NSMC仅公布了与GPS掩星相关的数据产品,因此,本文的反演及验证均是针对GPS掩星探测数据进行。由于大气密度、干气压及干温廓线均由折射率廓线计算,故以折射率廓线作为对比验证对象。本文将2017年自主反演的和NSMC公布的FY-3C大气折射率产品的质量进行了比较,NSMC产品的探测高度区间为地表至40~60 km,垂直分辨率为150~300 m。

      本文的质量评估以ECMWF提供的ERA-Interim再分析数据为参考,所使用的再分析产品的水平分辨率为0.75°×0.75°,时间分辨率为6 h,垂直方向为37个气压层,对应高度区间在0~50 km左右[7]。再分析数据并没有直接提供折射率参数,首先,利用再分析资料中的比湿、气压计算获取水汽压;然后,由水汽压、气压、温度进一步得到折射率。水汽压计算式为:

      e=Pqε+1-εq ]]>

      式中,P为大气压强,单位为hPaq为比湿,单位为g/kgε=Rd/RvRdRv分别表示干、湿空气气体常数,ε0.622。得到折射率为:

      N=77.6PT+3.73×105eT2 ]]>

      式中,T为大气的绝对温度,单位为K

    • FY-3C卫星上搭载的GNOS接收机提供GNSS卫星的双频相位及振幅数据。在无线电波穿过电离层及对流层时,由于其介质垂直密度的变化导致电波信号发生不同程度的弯曲,载波相位发生延迟。在低对流层,由于水汽含量丰富、多路径效应明显,掩星信号易发生失锁或跟踪错误,GNOS接收机采用开环跟踪模式解决该问题。由FY-3C附加相位数据文件出发,结合CIRA86+Q大气模型,进行FY-3C折射率自主反演的具体流程如图 1所示。

      图  1  FY-3C掩星折射率反演流程图

      Figure 1.  Flowchart for Retrieving the FY-3C Radio Occultation Refractivity Profiles

      图 1中,CIRA86+Q模型是国际空间研究委员会利用大量的地基、卫星等手段测量的风场和温度,给出的纬度范围为80°S~80°N、高度范围为0~120 km的风场和温度模型。反演关键步骤为:(1)在利用开环相位模型确定开环跟踪转换时间的基础上,结合L1模型附加相位、L1振幅数据,通过降频转换、移除调制导航资料及连接相邻相位以解决L1附加相位序列的半周跳、周跳等问题[8],实现对L1附加相位进行重构,并进而得到L1多普勒频移序列;(2)L2附加相位序列的质量一般比L1差,结合基于L2模型附加相位得到的L2模型多普勒频移,实现对L2多普勒频移的重构;(3)在上对流层,由重构后的双频多普勒频移出发,结合卫星位置信息,采用经典几何光学反演法计算大气弯曲角廓线,并通过双频组合进行弯曲角电离层影响修正;(4)对于低对流层,不同于NSMC所采用的第‍2类正则变换法,本文在自主反演中使用了FSI方法[9]。该方法由重构后的L1多普勒频移出发,直接从完整的复信号进行单一的Fourier变换得到瞬时频率;(5)两种方法反演的弯曲角廓线的拼接高度根据L1多普勒频移与模型多普勒频移的差异来确定,当采样点多普勒频移小于6 Hz,且与该点的多普勒频移模型值之差小于1 Hz时,则确定该点为弯曲角的拼接高度;(6)基于CIRA86+Q模型进行弯曲角上边界改正,利用Abel积分变换获取折射率廓线[10]

    • 在将自主反演的和NSMC提供的折射率产品分别与ERA-Interim再分析资料进行对比分析之前,需将掩星折射率廓线与ERA-Interim再分析资料进行时空匹配。本文采用的时间匹配为±3 h,空间匹配将ERA-Interim再分析资料采用三次样条插值后通过同一高度的4点均值方法匹配至相应掩星数据点。

      在对比分析中,折射率相对偏差δN、折射率平均相对偏差δN¯、折射率相对偏差标准差σN分别为:

      δNROi=NRO-NECMWFNECMWFδN¯ROi=1ni=1nδNROiσNi=1ni=1n(δNROi-δN¯ROi)2 ]]>

      式中,i为折射率廓线事件编号;n为折射率廓线事件总数量;RO=W,该处表示了两类廓线类型,故用RO(radio occultation)代替WN表示自主反演的或NSMC提供的折射率廓线事件。

    • 一次FY-3C掩星事件折射率廓线反演的个例如图 2所示。这是一次下降掩星事件,发生的时间为2017-01-30T23:56:00,位置为159.95°E,48.83°N发生掩星的GPS卫星为G16。该掩星事件的附加相位数据文件为“FY3C_GNOSX_GBAL_L1_20170130_2356_AEG16_MS.NC”。

      图  2  FY-3C折射率廓线反演及验证个例

      Figure 2.  Example for the Retrieval and Validation of the Refractivity Profile from FY-3C Radio Occultation Event

      图 2(a)给出了该事件的双频初始附加相位、L1模型附加相位和利用开环数据重构的L1附加相位随时间变化的序列,以及L1初始附加相位与重构附加相位的差值序列;图 2(b)给出了L1的振幅数据及开环时间;图‍2(c)利用重构后的双频附加相位得到的L1、L2附加多普勒频移序列、L1模型多普勒频移序列和L1附加多普勒频移与模型多普勒频移的差值序列;图‍2(d)以lg形式给出了L1、L2弯曲角廓线及经过电离层改正、优化后的弯曲角廓线和经过全谱反演算法反演得到的弯曲角廓线;图‍2‍(e)给出了自主反演的折射率廓线(W)、该掩星事件对应的NSMC折射率产品(N)和时空匹配的ERA-Interim折射率廓线(ERA)对比;图 2‍(f)则以ERA-Interim折射率廓线为参考,自主反演结果和NSMC的折射率产品各自的相对偏差廓线。

      图 2(a)中L1的重构附加相位序列是以L1原始附加相位序列为基础,结合L1模型附加相位序列及图 2(b)给出的L1振幅序列和开环时间得到的。由图 2(a)中的L1重构附加相位与原始附加相位的差值序列可见,在40 s之后,二者差异逐渐增大。由图 2(c)可见,基于重构L1附加相位序列得到的L1多普勒频移与L1模型多普勒频移的差值在-8~2.5 Hz之间变化。在反演中,当该差值绝对值超过5 Hz时,该点多普勒被认为是异常值,需进行剔除后内插。由图 2(a)可知,初始L2附加相位序列止于31 s,故图 2(c)中的L2多普勒频移序列在31 s之后是基于模型重构得到的。由图 2(d)可见,几何光学反演方法和FSI反演方法在高度20 km左右进行弯曲角拼接,获取的弯曲角廓线没有出现多值问题,解决了低对流层多路径效应的影响。由图 2(e)可知,自主反演廓线和NSMC产品都与ERA-Interim再分析折射率廓线有很好的一致性。图 2(f)表明,以ERA-Interim再分析数据为参考,自主反演的折射率廓线与NSMC结果的相对偏差基本一致。在高度40 km以下,自主反演的折射率相对偏差在±2%以内,在8~12 km及30~40 km的自主反演结果相对偏差小于NSMC结果。

    • 2017年参与反演的FY-3C的掩星事件(F)、NSMC公布的折射率产品廓线‍(N)及本文自主反演的折射率廓线(W)的数量柱状图如图 3所示。

      图  3  2017年折射率廓线数对比

      Figure 3.  Comparison of the Number of the Refractivity Profiles in 2017

      图 3可知,2017年NSMC公布的反演成功率为83.5%,自主反演的成功率为82%,总体反演成功率基本一致。图 4为以ERA-Interim再分析为参考,2017年NSMC提供的和自主反演的FY-3C折射率廓线的统计验证结果。其中,图 4(a)4(b)分别为NSMC产品和自主反演产品的统计结果。

      图  4  NSMC和自主反演折射率廓线统计验证

      Figure 4.  Statistical Validations of the Refractivity Profiles Provided by NSMC and Those Retrieved Independently

      图 4(a)可见,NSMC折射率产品在25 km以下折射率相对偏差绝对值小于2%,相对偏差标准差绝对值小于3%,25 km以上相对偏差及标准差逐渐增大,在35 km处相对偏差达到3%。由图 4(b)可见,自主反演的折射率廓线在10 km以下存在明显负偏差,折射率相对偏差绝对值在25 ‍km以下小于1.5%,略优于NSMC的相对偏差统计结果,但相对偏差标准差比NSMC产品大,其绝对值达到了3.6%。自主反演和NSMC反演过程中在低对流层采用不同的反演方法,以及质量控制标准的不同,都可能带来统计验证结果的差异。图 4(a)图 4(b)中的绿色点线表明参与统计比较的廓线数量随着高度的增加逐渐增多,在15 km以上趋于稳定,表明自主反演和NSMC的折射率产品的探测深度在低对流层基本一致,但自主反演成功的廓线数略低于NSMC产品数量。

      图 5图 6分别给出了以ERA-Interim折射率为参考,FY-3C掩星的NSMC折射率产品和本文的自主反演折射率廓线在不同纬度带的统计比较情况。图 5图 6中5个子图分别对应5‍个纬度带的统计比较,即低纬度带(30°N~30°‍S)、北半球中纬度带(30°N~60°N)、南半球中纬度带(30°‍S~60°S)、北半球高纬度带(60°N~90°‍N)和南半球高纬度带(60°S~90°S)。

      图  5  FY-3C掩星NSMC折射率产品不同纬度带的统计比较

      Figure 5.  Statistical Comparison of NSMC FY-3C Refractivity Products over Different Latitudes

      图  6  FY-3C掩星自主反演折射率廓线不同纬度带的统计比较

      Figure 6.  Statistical Comparison of FY-3C Refractivity Profiles Retrieved Independently over Different Latitudes

      图‍‍5‍(a)可见,在低纬度地区,NSMC折射率产品在低对流层的相对偏差达到2.5%,相对偏差标准差最大值为5%;在10~20 km高度区间,相对偏差及相对偏差标准差均小于2%;而在20 km以上,相对偏差及相对偏差标准差逐渐增大。由图 5(b)图 5(c)可见,在两个半球的中纬度地区,NSMC折射率产品在低对流层的相对偏差及相对偏差标准差均明显低于低纬度地区,其相对偏差标准差绝对值在2.2%以内;在20 km以上,相对偏差及相对偏差标准差逐渐增大,其中,南半球中纬度地区相对偏差及相对偏差标准差均略大于北半球中纬度地区。由图 5(d)图 5(e)可见,NSMC折射率产品在两个半球高纬度地区的相对偏差及标准差均明显低于低纬度地区及高纬度地区。

      图 6(a)可见,在低纬度地区,自主反演折射率廓线在2 km以上的低对流层高度区间,相对偏差绝对值小于1%,相对偏差标准差绝对值小于4%,优于NSMC公布的该区域折射率产品结果。由图 6(b)图 6‍(c)可见,在两个半球的中纬度地区,自主反演折射率廓线在低对流层高度区间存在较明显的负偏差,相对偏差标准差绝对值小于3%;在20 km以上区域,相对偏差及标准差逐渐增大,且南半球中纬度地区偏差明显大于北半球中纬度地区。由图 6(d)图 6(e)可见,在两个半球的高纬度地区,高度在10 km以下,自主反演的折射率廓线同样存在较明显的负偏差;高度在20 km以下,相对偏差标准差绝对值小于2%,略大于NSMC折射率产品的相应统计结果。

      对比图 5图 6可见,在低纬度地区,自主反演折射率廓线的质量总体上略优于NSMC产品;在中纬度地区,自主反演折射率廓线和NSMC产品的质量基本一致;在高纬度地区,NSMC产品质量略优于自主反演产品。

    • 对FY-3C掩星的NSMC折射率产品和自主反演的折射率廓线的相对偏差特性进行不同季节的对比分析。按照南北半球的不同季节划分,北半球3月—5月、6月—8月、9月—11月、12‍月至次年2月,南半球9月—11月、12‍月至次年2‍月、3月—5月、6月—8月分别定义为春、夏、秋、冬4‍个不同季节。

      NSMC折射率产品相对于ERA-Interim再分析资料的相对偏差在2017年的春、夏、秋、冬不同季节随纬度和高度变化的情况如图7(a)~7(d)所示。

      图  7  NSMC的FY-3C掩星折射率产品相对偏差不同季节的纬度-高度分布

      Figure 7.  Latitude-Height Distribution of the Relative Biases of FY-3C Radio Occultation Refractivity Products Provided by NSMC in Different Seasons

      图 7(a)和‍7‍(c)可见,在春、秋两季,NSMC折射率产品在20 km高度以下的相对偏差在-‍1.6%~0.8%,其中,南北纬60°以上相对偏差低于其他纬度值;高度在20 km以上,相对偏差逐渐增大,在中低纬度地区30 km高度左右,相对偏差达到4%。由图 7(b)7(d)可见,NSMC折射率产品在夏季2 ‍km高度以下为-‍1.6%~-0.4%,在高度2~15 km有-‍0.4%~0.4%的相对偏差,在冬季有-0.8%~0.4%的相对偏差,夏季的折射率负偏差比冬季更显著。

      自主反演的FY-3C掩星折射率廓线相对偏差不同季节的纬度-高度分布如图 8所示。自主反演的FY-3C掩星折射率廓线相对于ERA-Interim再分析资料的相对偏差在2017年的春、夏、秋、冬4个不同季节随纬度和高度变化的情况如图8(a)~8(d)所示。

      图  8  自主反演的FY-3C掩星折射率廓线相对偏差不同季节的纬度-高度分布

      Figure 8.  Latitude-Height Distribution of the Relative Biases of FY-3C Radio Occultation Refractivity Products Retrieved Independently in Different Seasons

      图 8(a)8(c)可见,在春秋两季,自主反演折射率相对偏差在2 km高度以下为-‍2.4%~-0.8%,2~15 km高度区间为-0.8%~0%,15 km高度以上为0.4%~4.0%。由图 8(b)8‍(d)可见,自主反演折射率相对偏差总体上在夏季明显高于冬季,低对流层高度区间夏季的相对偏差为-2.4%~-‍0.8%,而2~10 km高度区间冬季的相对偏差接近0%。

      对比图 7图 8可见:(1)在4个季节中,自主反演结果和NSMC结果,中低纬地区都比高纬地区偏差更大;在25 km以上,中低纬地区正偏差可能达到4%。(2)自主反演折射率廓线的偏差在2~15 km高度区间略优于NSMC的折射率产品。

    • 本文从NSMC提供的FY-3C掩星附加相位数据出发,采用几何光学反演和全谱反演算法相结合的方法,反演得到2017年全年的掩星大气折射率廓线。以ERA-Interim再分析数据为参考,对NSMC公布的和本文自主反演的FY-3C掩星折射率产品分别进行了统计验证,对两类折射率产品的相对偏差及其标准差进行了对比分析,并对不同纬度带和高度区间的误差特性及相对偏差的纬度-高度分布的季节变化进行了讨论,得到如下结论:

      1)对2017年全年的FY-3C掩星数据,自主反演折射率廓线的成功率为82%,与NSMC产品的反演成功率83.5%基本一致。

      2)虽然自主反演和NSMC分别采用了全谱反演和第2类正则变换方法获取低对流层的弯曲角廓线,但最终得到的折射率廓线的精度基本一致。

      3)以ERA-Interim再分析为参考,在25 km高度以下,NSMC的FY-3C掩星折射率产品相对偏差和其标准差分别小于2%和3%;自主反演的FY-3C折射率廓线的相对偏差和其标准差略低于NSMC产品。

      4)以ERA-Interim再分析为参考,NSMC提供的和自主反演的FY-3C掩星折射率廓线在高纬地区的偏差均低于中低纬地区的偏差,低纬地区偏差最大。FY-3C掩星折射率廓线偏差纬度-高度分布的季节变化明显,而自主反演折射率的相对偏差总体上在夏季高于冬季。

参考文献 (10)

目录

    /

    返回文章
    返回