Loading web-font TeX/Math/Italic
 注册 登录 English Version
国产遥感卫星历史数据再定标共性技术 | 浏览量 : 791 下载量: 506 CSCD: 0 更多指标
  • 导出

  • 分享

  • 收藏

  • 专辑

    • FY-3C微波湿度计-Ⅱ数据的再定标模型

    • Recalibration Model of MWHTS’ historical data onboard FY-3C Satellite

    • 王振占

      1

      肖雨伟

      12

      王可昕

      12

      张升伟

      1
    • 2023年27卷第10期 页码:2283-2294   

      纸质出版日期: 2023-10-07

    • DOI: 10.11834/jrs.20221446     

    扫 描 看 全 文

  • 引用

    阅读全文PDF

  • 王振占,肖雨伟,王可昕,张升伟.2023.FY-3C微波湿度计-Ⅱ数据的再定标模型.遥感学报,27(10): 2283-2294 DOI: 10.11834/jrs.20221446.
    Wang Z Z,Xiao Y W,Wang K X and Zhang S W. 2023. Recalibration Model of MWHTS’ historical data onboard FY-3C Satellite. National Remote Sensing Bulletin, 27(10):2283-2294 DOI: 10.11834/jrs.20221446.
  •  
  •  
    论文导航

    摘要

    星载微波湿度计观测的地球表面和大气的亮温是数值天气预报和气候研究的重要基础数据源。由于微波湿度计系统响应特性差异、在轨工作条件的变化,以及星上定标基准的差异,不同卫星微波湿度计、同一微波湿度计数据在长时间范围内都存在一定的偏差。如果这些偏差不进行修正,则会对长时间序列应用(如气候研究等)具有较大的影响,甚至可能导致错误的结果。由于微波辐射目前还不存在绝对基准,无法通过地面绝对基准的传递进行星上数据的绝对定标,因此需要利用其他相对定标方法对湿度计的稳定性及其变化特征进行评价。本文基于OMB(模拟和观测数差)定标的思想,结合微波湿度计的系统响应特性和星上定标基本原理,建立了微波湿度计定标偏差的校正模型—再定标模型,从而实现历史数据的一致性定标。初步研究结果表明:使用再定标模型后,FY-3C微波湿度计-II定标结果的时间序列稳定,原始OMB偏差中的各种异常波动现象得到良好修正。由于该模型是基于微波湿度计的系统响应特性建立的,因此这一方法同样适用于其他FY-3系列卫星微波湿度计数据的定标偏差校正,实现FY-3系列卫星微波湿度计数据的一致性定标。

    Abstract

    The microwave radiation data of the earth’s surface and atmosphere observed by the spaceborne microwave radiometer are a basic source for numerical weather prediction and climate research. Certain bias exists between the data of different or the same radiometers over a long period of time due to the difference in the response characteristics and calibration absolute reference between each sounder and variation of the on-orbit working conditions. If these biases are not corrected, then they will have significant implications for long time series applications (such as climate studies) and may even lead to erroneous results. In addition, absolute calibration of onboard data cannot be performed because the absolute reference for calibration is unobtainable at present. Therefore, other relative calibration methods should be used to evaluate the stability and variation characteristics of the radiometer. In this paper, a correction model for the calibration biases of the microwave humidity sounder is established, considering the principle of system responses and on-orbit calibration, which is also called re-calibration methods, based on OMB technique (observed brightness temperature (O) minus simulated brightness temperature (B)). Using the re-calibration methods, L1 level observation data of the FY-3C microwave humidity sounder were analyzed and the re-calibration coefficients were calculated. Then, the re-calibration coefficients were verified using two data sets. The residual and validation results show that the time series of the recalibration findings is stable, and some abnormal fluctuations in the original OMB biases have been well corrected after re-calibration, indicating that the re-calibration methods are valid. The method is based on the system responses of microwave humidity sounder. Therefore, it is also applicable to other FY-3 series satellite microwave humidity sounders for calibration bias correction. When using the same Radiative Transfer (RT) model, this re-calibration methods can achieve consistent calibration of FY-3 series satellite microwave humidity sounder data.

    关键词

    FY-3C微波湿度计; 定标偏差; OMB定标; 系统响应特性; 星上定标; 再定标模型; 卫星数据的一致性

    Keywords

    FY-3C microwave humidity sounders; re-calibration methods; system responses; consistency of satellite data

    1 引言

    长期的、一致的微波探测仪观测数据集对全球气候监测、数值天气预报等具有重大的影响。自2008年起,中国已经发射了风云三号(FY-3)A、B、C和D 4颗气象卫星,FY-3A和FY-3B上搭载的为微波湿度计MWHS(Microwave Humidity Sounder),FY-3C和FY-3D上搭载的为微波湿度计-II MWHTS(Microwave Humidity and Temperature Sounder),4颗卫星的湿度计已积累了超过10年的有效历史观测数据,利用这些数据,可以建立中国长期一致和稳定的大气湿度观测数据集(

    李娇阳 等,2019)。但是由于4颗卫星上的4台微波温度探测仪的系统响应特性和定标方法不尽相同,需要对这些有效历史观测数据进行一致性定标。

    交叉定标法(常用星下点同步交叉定标法SNOs(Simultaneous Nadir Overpasses)、微波探测仪观测亮温减去基于辐射传输模型计算的模拟亮温OMB(Observation Minus Background brightness temperature statistics)、双差分法DD(Double Difference technique))等是建立长期一致观测数据集常用的方法。

    目前,国内外已进行了一定的计算与研究。利用SNOs对NOAA(National Oceanic Atmospheric Adminstration)系列卫星的MHS(Microwave Humidity Sounder)和AMSU-B(the Advanced Microwave Sounding Unit-B)进行一致性定标(

    John等,2012)。利用Megha-Tropiques(Indo-French joint satellite)卫星上的Saphir(Sounder for Probing Vertical Profiles of Humidity)探测仪对NOAA-18、19以及Metop-A(Meteorological Operational-A)上的湿度计进行交叉定标,并使用DD方法对几台探测仪进行交叉定标(Jones等,2013)。利用SEVIRI(Spinning Enhanced Visible and Infrared Imager)探测仪作为参考探测仪,对NOAA-18和Metop-A的MHS进行长期一致的交叉定标(Buehler等,2020)。利用OMB、DD法对Metop-A、NOAA-19上的MHS进行交叉定标,使用NWP(Numerical Weather Prediction)模型进行亮温模拟(Saunders等,2013)。利用自然目标(南极和热带海洋)、区域平均(60°S到60°N)、SNOs 3种交叉定标方法对NOAA-15、16的MHS观测数据进行一致性定标(John等,2013)。

    SNOs是直接建立定标辐射计观测亮温与参考辐射计亮温之间的函数关系,方法较为简单,但具有纬度局限性,而且本身容易受到参考辐射计亮温稳定性的影响,定标结果不容易区分两类辐射计数据各自的定标精度及其稳定性。OMB和DD交叉定标法则是采用长时间序列数据的平均或者直方图进行分析,这种方法在一定程度上忽略了对辐射计系统的真实响应特性的分析,也没有考虑辐射传输模拟的正确性对于定标结果的影响。OMB的方法是直接把背景场参数带入微波辐射传输模型,获得卫星同步观测的亮温,进而与卫星数据进行比较,从而校正定标偏差的过程,这种方法虽然受到背景场参数、辐射传输模型精度等影响,但它可以在全球尺度内进行大量匹配数据的比对,具有很好的统计稳定性,因此被广泛用于在轨数据的定标。

    本文基于OMB方法的累积,通过分析不同卫星的系统响应特性,结合辐射传输方程的计算,确定测量电压与亮温之间的系统响应特性,进一步分析系统响应的时空变化特征,从而利用模拟亮温随测量电压之间的函数关系建立湿度计定标模型。

    2 研究方法或原理

    2.1 微波辐射计特性基本原理

    FY-3C发射于2013年9月,运行轨道的高度为836 km(Gou等,2015)。其上搭载的微波湿度计为第二代微波湿度计(MWHTS),包括15个探测通道,均使用双边带接收机。

    表1给出了其15个通道上下两个接收机的中心频点以及对应的极化方式,该极化方式随着扫描角θ的变化而变化。在θ为0°是极化是单一极化,其他角度都是混合极化。上下两个边带接收机接收辐射亮温的原理可表示为

    TB, p=vd2vd1fdown(v)TB, p(v)dv+vu2vu1fup(v)TB, p(v)dv
    (1)

    式中,fdown(v)fup(v)为下边带和上边带的通道响应函数,vd1vd2为下边带的频率边界,vu1vu2为上边带的频率边界。P表示极化方式。

    表1  MWHTS通道特点
    Table 1  Channel Characteristics of MWHTS
    中心频点/GHz通道极化方式备注
    87.75 1 θ 上边带
    90.25 1 θ 下边带
    113.75 9 90-θ 上边带
    115.75 8 90-θ 上边带
    116.25 7 90-θ 上边带
    117.65 6 90-θ 上边带
    117.95 5 90-θ 上边带
    118.45 4 90-θ 上边带
    118.55 3 90-θ 上边带
    118.67 2 90-θ 上边带
    118.83 2 90-θ 下边带
    118.95 3 90-θ 下边带
    119.05 4 90-θ 下边带
    119.55 5 90-θ 下边带
    119.85 6 90-θ 下边带
    121.25 7 90-θ 下边带
    121.75 8 90-θ 下边带
    123.75 9 90-θ 下边带
    148.76 10 θ 上边带
    151.24 10 θ 下边带
    176.31 15 90-θ 上边带
    178.81 14 90-θ 上边带
    180.31 13 90-θ 上边带
    181.51 12 90-θ 上边带
    182.31 11 90-θ 上边带
    184.31 11 90-θ 下边带
    185.11 12 90-θ 下边带
    186.31 13 90-θ 下边带
    187.81 14 90-θ 下边带
    190.31 15 90-θ 下边带
    icon 下载:  导出CSV icon 下载:  下载表格图片

    FY-3C微波湿度计-II采用两个反射面天线,在与卫星前进方向垂直的竖直面内进行交轨扫描,实现±53.35°扫描角的观测范围,其中,中心频率为89 GHz的通道1和中心频率为118.75 GHz的通道2—9共用扫描天线,天线波束宽度约为2°;中心频率为150 GHz的通道10和中心频率为183.31 GHz的通道11—15共用扫描天线,天线波束宽度约为1.1°(

    Zhang等,2014)。由于天线反射面在旋转,而馈源固定不动,因此微波湿度计接收的亮温是两个极化亮温的组合(Weng等,2003),即:

    TB=cos2θTBp+sin2θTBq
    (2)

    式中,θ是扫描角,当观测天底点时θ=0°TBpTBq表示两个正交极化的亮温。每个通道的亮温是这两个正交极化亮温的组合。因此,通道1和10的亮温可表示为

    TB, other=cos2θTBh+sin2θTBv
    (3)

    其他通道的亮温为

    TB, 1, 10=cos2θTBv+sin2θTBh
    (4)

    2.2 微波湿度计的天线温度和亮温

    微波湿度计计测量的天线温度TA可以表示为3种形式:

    TA=TH-TCVH-VCVA+VHTC-VCTHVH-VC+μρACρAH(TH-TC)2
    (5)
    TA=ρAC(TH-TC)+TC+μρACρAH(TH-TC)2
    (6)
    TA=ρACTH-ρAHTC+μρACρAH(TH-TC)2
    (7)

    式中,

    ρAC=VA-VCVH-VC
    (8)
    ρAH=VA-VHVH-VC
    (9)

    式中,VAVHVC分别为观测地球场景电压、热源电压和冷源电压,THTC分别是热源和冷源的温度,单位为K。u为非线性系数,单位为m2Srcm-1mW-1,通常μ是接收机温度的函数,在发射前通过热真空实验,测量不同仪器温度下的系统响应获得。在轨定标进行非线性修正时通过查表来获得实际工作温度下的μ值。当考虑天线方向图和天线辐射效率对于天线温度贡献的情况下,卫星接收的亮温TB, p可以简化为

    TB, p=aTA, p+bTA, q+cTant+d
    (10)

    式中,TA, pTA, q表示表示同极化和交叉极化的天线温度,Tant表示天线的物理温度,abcd为系数。对于大气探测通道,通常只有一个极化方式,且这两种极化的天线温度非常接近,因此从实际数据处理的角度出发,上式可以进一步简化为

    TB=aTA+bTant+c
    (11)

    式(11)忽略了交叉极化贡献。把式(5)—(7)代入式(11)即可得到测量亮温的3种表达形式(为了简化篇幅这里从略),作为在轨定标偏差校正的基础。下文将以分析在轨数据特征及其与辐射传输计算偏差为基础,构建定标偏差校正模型—再定标模型。

    2.3 FY-3C微波湿度计-Ⅱ在轨数据及其OMB结果分析

    使用OMB方法进行数据质量评价和定标误差分析主要是基于数据应用的考虑:一方面,由于微波辐射计的观测亮温至今没有绝对的参考基准,通过星上两点定标的结果可能存在一定的偏差。通过辐射传输模型和背景场计算的亮温具有统计稳定性,为卫星观测数据提供了可进行比较的参考基准,并且微波辐射传输模型计算的亮温作为参考基准可以实现不同载荷之间参考基准的一致性传递;另一方面,微波辐射计数据的应用是基于微波辐射传输模型的,利用OMB方法得到的再定标数据和所使用的微波辐射传输模型之间是无偏的,可以直接用于基于同一模型的后续反演和同化应用,不需要再进行偏差修正。如果不使用当前微波辐射传输模型进行反演和同化,还是需要进行偏差校正的,这是因为不同微波辐射传输模型之间仍然存在偏差。总之,利用OMB方法获得的再定标数据是稳定的,相对当前微波辐射传输模型来说是无偏的,相对其他微波辐射传输模型可能还存在稳定的偏差。

    为了计算模拟背景亮温,这里使用MPM93(

    Liebe,1989Liebe等,1993)逐线计算大气吸收系数模型和FASTEM5(English和Hewison,1998Liu等,2011Bormann等,2012)的海面发射率模型。计算的输入背写数据来自ERA5的每小时的再分析数据集(Dee等,2011)。在模拟两个正交极化不同扫描角亮温以后,根据式(2)—(4)计算对应通道的亮温。在积分过程中,通道响应函数设定为简单的矩形,上下边带的权重设定为0.5(王振占,2023)。后续计算将带入真实权重函数进行仿真。此处选取2013年11月到2016年12月卫星湿度计的L1级数据产品作为分析数据。

    根据再定标前OMB值随时间(扫描线)的变化结果,因为有些通道的变化具有共性,图1只展示了通道7、14的结果作为代表进行说明再定标前OMB值随时间(扫描线)变化特征。从图1可以看出,以通道14为代表的通道4、6、11、13、14,其再定标前OMB值随时间变化趋势存在明显的波动,以通道7为代表的其余通道总体OMB值随时间呈现稳定的变化趋势。对于通道4、6、11、13、14的OMB值随时间变化波动的原因需要识别,并建立校正算法。

    fig
    icon 下载:  | 高精图 | 低精图

    图1  通道7、14再定标前观测亮温与模拟亮温的差值

    Fig. 1  Biases of observed and simulated brightness temperature before re-calibration in 7 and 14 channels

    对这些数据进行进一步分析,发现引起偏差波动的主要原因是仪器温度的变化和接收机自动增益控制电压(AGC)的变化。图2展示了FY-3C微波湿度计在2013年—2017年期间接收机射频和中频温度的变化。可见在这个时间段,仪器温度有几个较为明显的跳变,为卫星系统操控的结果。部分通道(如13、14通道)的偏差波动与仪器温度变化有一定的相关性。

    fig

    图2  2013年—2017年期间FY-3C微波湿度计仪器射频和中频温度的变化

    Fig. 2  Time varying instrument RF and IF of FY-3C MWHTS during 2013—2017

    icon 下载:  原图 | 高精图 | 低精图

    15个通道AGC随着时间的变化如图3所示。其中,通道6的AGC在2014年—2015年期间有一个大的跳变;对于其它通道,在发射初期(2013年9月)和2015年初期也都有一个较为明显的AGC变化;AGC的变化与通道4和5偏差的跳变具有相关性。AGC变化导致接收机非线性系数发生变化,而接收机的非线性系数又是仪器工作温度的函数,因此可以从仪器温度和定标源温度的变化分析定标偏差的变化原因,从而建立再定标算法。

    fig

    图3  2013年—2016年期间FY-3C微波湿度计在15个通道在轨AGC的变化

    Fig. 3  Time varying AGC of MWHTS 15 channels during 2013-2017

    icon 下载:  原图 | 高精图 | 低精图

    3 再定标模型的建立及结果分析

    3.1 再定标模型

    结合式(10)和(11),定标偏差dTB

    dTB=TBsim-TBobs
    (12)

    式中,TBsim表示模拟亮温,TBobs表示观测亮温。当综合考虑冷热定标源误差、非线性误差以及天线温度的贡献时,可得到如下的表达形式:

    dTB=aρAC+bρAH+ρACρAH(c(TH-TC)TIF+d(TH-TC)+eTIF+f)+gTP+h
    (13)

    式中,a,b,c,d,e,f,g,h为系数。如果忽略冷空亮温的变化,同时由于微波湿度计的天线和定标源处于相同的位置,其物理温度不是独立的,因此与热源温度合并考虑为一个影响因素,则式(13)可以简化为

    dTB=TBsim-TBobs=aVA-VCVH-VC+bTIF+c
    (14)

    式中,a,b,c为再定标方程的系数,用来解释由于热源温度、接收机温度变化产生的定标误差。

    为了提高算法的稳定性和代表性,使用82天的L1数据与ECMWF的ERA5每小时的再分析数据集进行匹配。匹配标准:使用南北纬60°的晴空海面亮温,根据湿度计L1产品中的地表类型为海洋进行筛选,晴空判断标准为积分液态水含量小于0.1,海面亮温大于275 K以去除海冰的影响;时间上取观测轨道首尾时间上最为接近的EC数据,时间匹配标准为小于±30 min,图4为FY-3C的L1级观测数据与ERA5每小时在分析数据集经过时间匹配后两者的时间差,由图4可以看出,经数据匹配后得到的数据级在时间差范围(±0.5小时内)呈整体均匀分布,保证了卫星数据观测时间与ERA5时间的无偏性;在全球经纬度的范围内,ERA5数据是分辨率为0.25°×0.25°的网格形式,将L1级数据的经纬度对应到此网格中,经纬度的匹配标准小于±0.125°。

    fig

    图4  FY-3C L1级观测数据与ERA5每小时再分析数据集进行时间匹配后两者的时间差

    Fig. 4  The time difference between L1 level observed brightness temperature and ERA5 hourly reanalysis data set after time matching

    icon 下载:  原图 | 高精图 | 低精图

    3.2 再定标结果分析

    对匹配好的数据进行分析,发现落在一个ERA网格通常有1—4次观测点,为计算数据的均匀性和有效性,选择了在一个网格内有3—4次观测点作为基准观测数据,并要求这些观测数据之间的标准差小于0.3 K。利用这些匹配的网格点的ERAS数据进行模拟亮温计算产生最终的再定标算法系数。经第2节分析,AGC的变化可能对再定标结果造成影响,同时也计算了输入数据在不同AGC下的定标系数。

    图5为再定标前后部分通道(通道7和14)观测亮温和模拟亮温之间的定标偏差和标准差随热源温度和时间的变化示意图,其中图5(a)(b)中“产品”图标代表着再定标前OMB值及其标准差,“区分”图标代表着在不同AGC下计算的多组再定标系数得到的再定标偏差和标准差,“不区分”图标代表着所有数据计算的一组再定标系数得到的再定标偏差和标准差;图5(a)(b)分别为通道7和通道14再定标前后OMB值随热源温度的变化、再定标前后OMB值随时间(天数)的变化及再定标前后OMB值的标准差随时间(天数)的变化。

    fig
    icon 下载:  | 高精图 | 低精图

    图5  再定标前后通道7、14观测亮温和模拟亮温之间的定标偏差、标准差随热源温度和时间的变化

    Fig. 5  The variation of calibration deviation and standard deviation of observed and simulated brightness temperature before and after re-calibration in ch7 and ch14

    图5(a)中间图可以看出,对于通道7,区分和不区分AGC下再定标结果均较好,OMB值的大小和波动都得到校正,但是区分AGC的再定标结果优于不区分AGC通道6、11、12结果与通道7的结果相似,在区分AGC的条件下,再定标后的偏差小于不区分AGC的定标偏差,并且偏差的波动得到更好的改善,限于篇幅原因,结果图不再列出。因此,对于通道6、7、11、12,这些通道AGC的变化对定标结果有一定的影响,在定标系数中需要考虑AGC的影响。

    图5(b)中图可以看出,对于通道14,区分和不区分AGC下计算的定标结果均较好,再定标后的偏差大小和波动都得到改善,两种情况下计算的结果基本没有差别。通道1、13、14、15再定标结果与通道14结果相似;通道2、3、8、9、10在不区分AGC下的再定标结果略好,限于篇幅原因,这些通道的结果图不再列出。因此,对于受AGC的影响不明显的通道,选择不区分AGC下的定标系数进行再定标。

    对于通道4、5,再定标后的OMB值大小得到校正,但是OMB值仍然存在波动,图6为再定标前后通道4观测亮温和模拟亮温之间的定标偏差和标准差随热源温度和时间的变化示意图。从图6中的中间时间变化图可以看出,无论区分和不区分AGC的定标结果,在30—50对应的时间范围,OMB值仍然存在相对较大的波动。通过分析原因发现,这可能是在星上定标时使用的非线性系数不合适。L1级观测亮温是L0级原始观测数据经过星上定标得到的,在星上定标过程中需要进行非线性修正,而在波动处该通道的中频温度发生变化,非线性系数u对中频温度的变化不敏感,导致在该处进行星上定标时非线性修正量不准确,因此导致L1级观测数据存在人为误差,进而造成的此处波动,后续会对非线性系数进行进一步讨论。通道5的在第27—30时间段OMB同样值存在相对较大的波动,波动原因与通道4一致,此处不进行详细分析。

    fig

    图6  再定标后通道4观测亮温和模拟亮温之间的定标偏差、标准差随热源温度和时间的变化

    Fig. 6  The variation of calibration deviation and standard deviation of observed and simulated brightness temperature before and after re-calibration in ch4

    icon 下载:  原图 | 高精图 | 低精图

    虽然通道4、5再定标结果不是很理想,但其余通道通过考虑定标源温度、仪器温度和AGC的变化,可以实现OMB偏差的大幅改善,实现了数据的稳定定标。为分析该定标模型的有效性,4、6、7、11、12通道采用L1数据区分AGC的再定标系数;其他通道采用L1数据不区分AGC的系数。再定标系数见表2

    表2  FY-3C MWHTS各通道再定标模型系数
    Table 2  The parameters of re-calibration processing of FY-3C MWHTS
    通道AGC回归系数a回归系数b回归系数c
    1 不区分 11.798 -0.0065994 -5.1504
    2 不区分 -60.853 -0.17126 96.373
    3 不区分 -7.8505 -0.049257 19.537
    4 5.0012 -11.231 -0.060992 27.164
    4 5.3114 1.4821 0.038085 -10.365
    5 不区分 -19.59 -0.03603 25.909
    6 3.0769 -0.42783 0.021256 -4.9395
    6 43.087 -0.78161 0.18674 -52.052
    6 3.2234 0.71813 0.062367 -17.871
    7 4.6886 -0.095851 0.049648 -13.227
    7 5.0012 -0.030544 0.031727 -8.7581
    8 不区分 1.4908 0.075928 -21.798
    9 不区分 3.2397 0.098696 -30.393
    10 不区分 -18.183 -0.034933 30.65
    11 3.5946 0.6611 0.12259 34.314
    11 3.5165 0.11654 0.15852 -44.488
    11 3.7167 1.3483 0.17011 -48.846
    12 4.0635 -1.2457 -0.053651 16.714
    12 4.105 0.40627 0.015025 -4.3331
    13 不区分 -8.0574 0.29274 -78.266
    14 不区分 -16.285 0.46035 -119.83
    15 不区分 -11.968 0.005388 11.352
    icon 下载:  导出CSV icon 下载:  下载表格图片

    图7为通道7、14再定标后OMB偏差随时间(扫描线)的变化。由图7(b)可以看出,通道14偏差波动得到改善。

    fig
    icon 下载:  | 高精图 | 低精图

    图7  通道7、14再定标后观测亮温与模拟亮温的结果比较

    Fig. 7  Comparison of observed and simulated brightness temperature after recalibration in ch7 and ch14

    为了定标结果的稳定性,我们下面主要分析算法的残差特性。图8为通道7、14再定标前后的OMB偏差直方图,从图8可以看出,定标前后的偏差直方图存在明显差异:如通道7虽然呈现一定的正态分布,但是它们明显存在两个分布特征,直方图的分布至少是两个状态的叠加,这对应于这些通道AGC的变化;通道14定标之前明显不是正态分布,存在两个峰值。而经过再定标以后,通道7、14的OMB值呈现理想的正态分布,证明了该方法的有效性。对于其他通道,OMB值也都呈现理想的正态分布,这里不在一一列出。

    fig
    icon 下载:  | 高精图 | 低精图

    图8  匹配观测点为3—4次时,通道7、14再定标前后的OMB偏差直方图比较

    Fig. 8  Comparison of the OMB bias histogram before and after re-calibration processing in ch7 and ch14

    4 再定标结果验证及分析

    4.1 利用同时间样本中的其他数据进行验证

    下面利用同一网格内匹配点只有1—2次的数据进行验证,以判断算法系数是否可以用于一般的匹配情况。一般情况中0.25°网格内的不同匹配数据可能包含位置差异、大气条件差异和计算误差等,而且数据之间的离散度越大,说明这些差异越大。因此,可以利用这些数据验证再定标算法系数是不是产生稳定的效果。由于匹配点1—2个的情况下无法判读像元内数据的稳定性,不能有效判断降水等在像元内是否存在,所以这些数据的离散程度略微增大。图9为当匹配观测点为1—2次时通道7、14再定标前后OMB偏差直方图。由图9可以看出虽然这两个通道的离散度略微增大,但是明显看出定标系数对于这些通道是完全适用的,很好地实现了定标偏差和异常数据的校正,对于其他通道特征与通道7、14相似。

    fig
    icon 下载:  | 高精图 | 低精图

    图9  匹配观测点为1—2次时通道7、14再定标前后OMB偏差直方图比较

    Fig. 9  Comparison of the OMB bias histogram before and after re-calibration processing in ch7 and ch14

    4.2 利用未覆盖时间段的数据进行验证

    为了检验定标结果,进一步选择其他独立的数据进行验证。在此,选择2017年—2019年期间60天的数据进行算法鲁棒性的验证。图10(a)为通道7再定标前后观测亮温与模拟亮温之间的偏差和标准差随天数的变化,通道7再定标前OMB值本身很小,再定标完之后OMB值相对来说变化不大。图10(b)为通道14再定标前后观测亮温与模拟亮温之间的偏差和标准差随天数的变化,可以看出利用再定标系数,OMB值变小,波动也得到改善。对于其他通道的再定标模型也都很好地实现了定标偏差的校正,使得定标偏差稳定,在所分析的跨度3年的60天范围内,定标偏差的波动较小。图11部分通道经过再定标之后的直方图,这两个通道都呈现理想的正态分布,其他通道也是如此,再次证明了该方法的有效性。

    fig
    icon 下载:  | 高精图 | 低精图

    图10  利用2017年—2019年期间60天数据进行再定标的验证结果

    Fig. 10  Validated results of the re-calibration processing using sixty days’ data during 2017—2019

    fig
    icon 下载:  | 高精图 | 低精图

    图11  利用2017年—2019年期间60天数据进行再定标前后OMB偏差直方图比较

    Fig. 11  Comparison of the OMB bias histogram before and after recalibration processing using sixty days’ data during 2017—2019

    4.3 定标残差的分布特性统计

    由于微波湿度计采用交轨扫描的形式进行观测,因此在98个观测位置上的偏差是否存在位置相关性一致都是大家关心的问题(

    Goldberg等,2001Tian等,2018)。这里对于再定标前后观测亮温与模拟亮温之间的偏差和标准差的分布特性进行统计。图12为FY-3C的MWHTS的所有通道再定标后观测亮温和模拟亮温之间的偏差随扫描点的变化趋势。在选取数据时,选取了3—4次观测点进行了模型的建立,位于边缘扫描点的数据个数不满足条件,造成了边缘位置的匹配的统计数据缺失。由图12可以看出,对于C星微波湿度计的所有通道,第16—88个扫描点的偏差得到了较好的修正,但再定标前后的偏差都随着扫描点的变化而变化,具体原因会在之后的工作中进行分析。图13为15个通道再定标前后观测亮温和模拟亮温之间的标准差随扫描点的变化趋势,再定标前后OMB的标准差变化很小,趋势也不变。

    fig

    图12  匹配观测点为3—4次时,所有通道再定标前后OMB偏差随扫描点的变化

    Fig. 12  Varying of OMB bias with scanning point before and after recalibration of all channels using 3—4 observed points

    icon 下载:  原图 | 高精图 | 低精图
    fig

    图13  匹配观测点为3—4次时,所有通道再定标前后OMB标准差随扫描点的变化

    Fig. 13  Varying of OMB std with scanning point before and after recalibration of all channels using 3—4 observed points

    icon 下载:  原图 | 高精图 | 低精图

    5 结论

    本文利用微波辐射计通道特性,对FY-3C湿度计历史数据进行再定标分析。基于OMB方法,直接利用辐射计系统响应特性,结合辐射传输方程来进行定标偏差分析,在考虑模拟和观测数据时空匹配的有效性时,不但依赖其统计特性的稳定性,而且在2.667 s定标周期内的需要通过湿度计实时观测的电压数据来评判时空匹配的有效性。所以本文基于对不同卫星的非线性参数、热源温度和电压以及匹配点的数据电压特性的综合考虑,建立新型的历史数据再定标算法。再定标后的OMB值明显减小,在统计特征上,OMB的值基本符合正态分布特性。从2013年—2017年4年期间选出代表性的82天数据偏差变化来看,通道4和13的最大偏差不超过1 K;其他通道再定标后的观测亮温和模拟亮温的最大偏差不超过0.5 K。除通道4、5这两个通道的波动本身就很小外,其余通道偏差的波动得到改善。从2017年—2019年3年期间抽取的代表性60天验证数据来看,通道3、5、6、7、10、11、12、15再定标后的观测亮温亮温和模拟亮温最大偏差不超过0.5 K,其余通道最大误差不超过1 K。对于再定标结果仍存在轻微波动的通道,后续会进一步分析非线性系数对再定标模型的影响。

    致谢

    致谢:感谢国家卫星气象中心提供的风云三号湿度计L1和L0数据,感谢ECMWF提供的ERA5每小时再分析数据。

    参考文献(References)

    Bormann N, Geer A and English S J. 2012. Evaluation of the microwave ocean surface emissivity model FASTEM-5 in the IFS. Bracknell: European Centre for Medium Range Weather Forecasts [DOI: 10.21957/l936s9df] [百度学术] 

    Buehler S A, Prange M, Mrziglod J, John V O, Burgdorf M and Lemke O. 2020. Opportunistic constant target matching—a new method for satellite intercalibration. Earth and Space Science, 7(5): e2019EA000856 [DOI: 10.1029/2019EA000856] [百度学术] 

    Dee D P, Uppala S M, Simmons A J, Berrisford P, Poli P, Kobayashi S, Andrae U, Balmaseda M A, Balsamo G, Bauer P, Bechtold P, Beljaars A C M, Van De Berg L, Bidlot J, Bormann N, Delsol C, Dragani R, Fuentes M, Geer A J, Haimberger L, Healy S B, Hersbach H, Hólm E V, Isaksen L, Kållberg P, Köhler M, Matricardi M, Mcnally A P, Monge-Sanz B M, Morcrette J J, Park B K, Peubey C, De Rosnay P, Tavolato C, Thépaut J N and Vitart F. 2011. The ERA-interim reanalysis: configuration and performance of the data assimilation system. Quarterly Journal of the Royal Meteorological Society, 137(656): 553-597 [DOI: 10.1002/qj.828] [百度学术] 

    English S J and Hewison T J. 1998. Fast generic millimeter-wave emissivity model//Proceedings of SPIE 3503, Microwave Remote Sensing of the Atmosphere and Environment. Beijing: SPIE [DOI: 10.1117/12.319490] [百度学术] 

    Goldberg M D, Crosby D S and Zhou L H. 2001. The limb adjustment of AMSU-A observations: methodology and validation. Journal of Applied Meteorology, 40(1): 70-83 [DOI: 10.1175/1520-0450(2001)040<0070:TLAOAA>2.0.CO;2] [百度学术] 

    Guo Y, Lu N M, Qi C L, Gu S Y and Xu J M. 2015. Calibration and validation of microwave humidity and temperature sounder onboard FY-3C satellite. Chinese Journal of Geophysics, 58(1): 20-31 [百度学术] 

    郭杨, 卢乃锰, 漆成莉, 谷松岩, 许健民. 2015. 风云三号C星微波湿温探测仪的定标和验证. 地球物理学报, 58(1): 20-31 [DOI: 10.6038/cjg20150103] [百度学术] 

    John V O, Allan R P, Bell W, Buehler S A and Kottayil A. 2013. Assessment of intercalibration methods for satellite microwave humidity sounders. Journal of Geophysical Research: Atmospheres, 118(10): 4906-4918 [DOI: 10.1002/jgrd.50358] [百度学术] 

    John V O, Holl G, Buehler S A, Candy B, Saunders R W and Parker D E. 2012. Understanding intersatellite biases of microwave humidity sounders using global simultaneous nadir overpasses. Journal of Geophysical Research: Atmospheres, 117(D2): D02305. [DOI: 10.1029/2011JD016349] [百度学术] 

    Jones L, Datta S, Santos-Garcia A, Wang J R, Payne V, Viltard N and Wilheit T. 2013. Radiometric intercalibration of the microwave humidity sounder on NOAA-18, MetOp-A, and NOAA-19 using SAPHIR on megha-tropiques//IEEE International Geoscience and Remote Sensing Symposium. Melbourne: IEEE: 1151-1154 [DOI: 10.1109/IGARSS.2013.6721369] [百度学术] 

    Zhang S W, He J Y, Wang Z Z and Wang X. 2014. In-orbit performance of microwave humidity and temperature sounder (MWHTS) of the Chinese FY-3C meteorological satellite. General Assembly and Scientific Symposium (URSI GASS), Beijing, China. 1-4[doi: 10.1109/URSIGASS.2014.6929598] [百度学术] 

    Li J Y, Wang Z Z, Gu S Y and Zhang S W. 2019. Common re-calibration technology for spaceborne microwave atmospheric humidity sounder. Remote Sensing Technology and Application, 34(6): 1212-1220 [百度学术] 

    李娇阳, 王振占, 谷松岩, 张升伟. 2019. 星载微波大气湿度探测仪再定标共性技术分析. 遥感技术与应用, 34(6): 1212-1220 [DOI: 10.11873/j.issn.1004-0323.2019.6.1212] [百度学术] 

    Liebe H J. 1989. MPM—An atmospheric millimeter-wave propagation model. International Journal of Infrared and Millimeter Waves, 10(6): 631-650 [DOI: 10.1007/BF01009565] [百度学术] 

    Liebe H J, Hufford G A and Cotton M G. 1993. Propagation modeling of moist air and suspended water/ice particles at frequencies below 1000 GHz//AGARD Conference Proceedings. Boulder: National Telecommunications and Information Administration: 542 [百度学术] 

    Liu Q H, Weng F Z and English S J. 2011. An improved fast microwave water emissivity model. IEEE Transactions on Geoscience and Remote Sensing, 49(4): 1238-1250 [DOI: 10.1109/TGRS.2010.2064779] [百度学术] 

    Saunders R W, Blackmore T A, Candy B, Francis P N and Hewison T J. 2013. Monitoring satellite radiance biases using NWP models. IEEE Transactions on Geoscience and Remote Sensing, 51(3): 1124-1138 [DOI: 10.1109/TGRS.2012.2229283] [百度学术] 

    Tian X X, Zou X L and Yang S P. 2018. A limb correction method for the microwave temperature sounder 2 and its applications. Advances in Atmospheric Sciences, 35(12): 1547-1552 [DOI: 10.1007/s00376-018-8092-8] [百度学术] 

    Wang Z ZXu H WDuan Y QWang W YDing JHe W M and Zhang S W. 2023. Effect analysis of spectral response function of microwave humidity and temperature sounder onboard the FY-3D satellite. National Remote Sensing Bulletin272):394-405 [百度学术] 

    王振占许皓文段永强王文煜丁甲何文明张升伟.2023.FY-3D星微波湿温探测仪通道响应函数的影响分析.遥感学报272): 394-405 [DOI: 10.11834/jrs.20220450] [百度学术] 

    Weng F Z, Zhao L M, Ferraro R R, Poe G, Li X F and Grody N C. 2003. Advanced microwave sounding unit cloud and precipitation algorithms. Radio Science, 38(4): 8068 [DOI: 10.1029/2002RS002679] [百度学术] 

    文章被引用时,请邮件提醒。
    提交

    相关文章

    暂无数据

    相关作者

    暂无数据

    相关机构

    暂无数据
    0