基于FY-4A卫星探测区域模式背景误差和观测误差估计的长波红外通道选择研究

中国新一代静止气象卫星风云四号A星（FY-4A）上搭载的干涉式大气垂直探测仪GIIRS（Geostationary Interferometric Infrared Sounder）是世界上首台在静止轨道运行的高光谱红外探测仪器,拥有1650个通道,其中长波通道689个,中波通道961个。由于计算机存储能力、数据传输和变分同化时效性等限制,目前很难在业务环境中同化全部通道,并且多通道卫星信息往往存在一定的空间相关和光谱相似,故在实际应用中必须降低高光谱资料的通道维数,去除数据冗余和观测相关性,挑选出对特定目标起主要作用的通道子集,进而利用有限的通道来最大限度地提供观测信息。考虑GRAPES（Global/Regional Assimilation and Prediction Enhanced System）数值预报系统估计得到的FY-4A观测范围内的背景误差协方差矩阵,并且在GIIRS观测误差后验重估计的基础上,结合通道权重函数峰值和信息熵分步迭代法对其长波通道进行了通道选择研究和试验。结果表明,在给定选择通道个数时,权重函数法和信息熵方法得到的通道组合进行温度和湿度误差分析时误差整体比单一方法小,在给定通道个数（40个）少于单一方案（50个）时,误差分析结果相当。GIIRS最优通道选择的研究为该资料的同化应用建立了基础。      

GNSS/PWV与风云四号A星GIIRS水汽廓线融合应用研究

中国新一代地球静止气象卫星风云四号A星（FY-4A）搭载的干涉式大气垂直探测仪（Geostationary Interferometric Infrared Sounder, GIIRS）以红外高光谱干涉分光方式探测三维大气温湿结构,取得了在静止轨道上探测大气的突破性进展。地基全球导航卫星系统（Global Navigation Satellite System,GNSS）是一种连续监测大气可降水量（Precipitable Water Vapor,PWV）的有效手段,基于2018年6—8月中国地基GNSS站监测的PWV和FY-4A/GIIRS水汽廓线的业务产品以及常规无线电探空资料,开展GNSS/PWV与FY-4A/GIIRS水汽廓线快速融合应用,以提高卫星资料反演大气水汽廓线的精度。结果表明:与常规无线电探空相比,FY-4A/GIIRS水汽廓线产品在大气低层均方根误差（Root Mean Square Error,RMSE）为4.5 g/kg,700 hPa为2.4 g/kg,500 hPa以上因水汽含量较低RSME小于1.5 g/kg。GNSS/PWV与FY-4A/GIIRS水汽廓线融合后,FY-4A/GIIRS水汽廓线误差整层RMSE减小20%,从近地层到600 hPa RMSE平均减小20%—25%,尤其是850—700 hPa改善最明显,极大改善了卫星水汽反演资料的可用性。对一次多系统影响的暴雨天气过程应用分析表明,GNSS/PWV和FY-4A/GIIRS融合产品可获得高时、空密度的大气水汽廓线,对强降水的临近预报有非常重要的支撑作用。      

FY-4A GIIRS高时间分辨率温湿度廓线反演及其在台风中的应用

针对2020年第9号台风"美莎克"期间FY-4A 高光谱红外干涉式大气垂直探测仪GIIRS每15 min一次的目标区跟踪加密观测资料,用三维卷积神经网络算法反演的全天空大气温度、湿度廓线分析了台风处于生命史不同阶段时暖心结构和湿度场结构的演变特征。结果表明:卷积神经网络的深度机器学习算法可以用来反演全天空的三维大气温度和湿度垂直廓线,不光适用范围广(晴空和有云视场)、反演精度高,而且反演速度快。利用静止卫星平台高时间分辨率的特性,反演得到的温度、湿度廓线可以细致追踪台风处于发展、成熟和登陆等阶段时暖心结构和湿度场的时空演变特征。台风从发展阶段(热带风暴和强热带风暴)到成熟阶段至登陆消亡时,暖心首先出现在对流层中高层较薄的区域,随着台风强度的加强,深厚的暖心结构明显、强度增加,水平面积增大且垂直往下延伸。由于对流云中强上升气流的输送水汽正距平区逐渐上传至300 hPa,台风最强时密闭云区与四周下沉气流区比湿差高到8 K·kg^-1。暖心结构和高湿度中心随着台风登陆而逐渐消失。     

冬奥赛场卫星和地基微波辐射计气温廓线观测对比和验证

利用探空站数据对北京和张家口冬奥赛场周边地区的地基微波辐射计和FY 4A大气垂直探测仪资料进行验证分析。选取2020年全年FY 4A大气垂直探测仪资料以及2020年12月至2021年3月期间对冬奥赛事有重要影响的寒潮前后、雾霾和沙尘暴这3种不同天气现象下探测得到的气温廓线数据进行个例分析。结果表明：在晴空条件下,地基微波辐射计和FY 4A大气垂直探测仪探测大气垂直气温的精度较高,平均相关系数达到0.97,低层大气（500 hPa以下）较高层大气探测结果的一致性较好;大气污染对探测精度产生一定影响,其中PM2.5产生的影响较小,PM10的提高对FY 4A大气垂直探测仪的探测产生较大影响,尤其是发生沙尘暴时,星载探测仪无法对低层大气进行探测。经过对比和验证,卫星探测作为补充探测手段,可以与地基微波辐射计互相补充,尤其是在空间覆盖和时间分辨率上具备一定优势,但在有云和沙尘暴的天气条件下无法对低层大气开展探测;地基微波辐射计可以对一个地点的大气垂直参数开展不间断的探测,与探空站数据的一致性较高。卫星探测和地基微波辐射计均可以为冬奥赛事提供高时间分辨率的探测数据,为数值天气预报提供有力的数据支撑。     

“珍珠”台风卫星红外通道亮温的数值模拟

热带地区云的覆盖率高, 卫星辐射资料特别是红外通道的数值模拟容易受到云的影响, 提高台风等热带天气系统卫星亮温资料的模拟能力, 有助于热带地区数值模式的检验以及提高卫星资料直接同化的水平。首先, 利用天气研究和预报模式(Weather Research and Forecasting Model, 简称WRF) 对 “珍珠” 台风的三维气象场, 包括云结构特征, 进行了模拟, 并诊断分析了总云量、云顶气压等物理量; 然后, 结合快速辐射传输模式 (fast radiative transfer model for TOVS, 简称RTTOV) 8.7版对 “珍珠” 台风的高分辨率红外辐射探测仪3型 (High Resolution Infrared Radiation Sounder, 简称HIRS/3) 红外辐射亮温进行了模拟, 并且比较了有云和无云条件下的模拟结果。结果表明, 如果不考虑云的影响, RTTOV模式无法模拟出台风的结构特征; 结合三维云特征的模拟量, RTTOV模式可以较好地模拟出HIRS/3第4～10、13～19通道的辐射亮温。加入云结构特征对峰值能量高度在云层以下的红外通道的模拟有正面影响, 该方案既可以作为检验模式模拟的一种手段, 同时也表明红外资料在热带天气系统的资料同化领域具有潜在的应用能力。     

面向资料同化的FY-4A 卫星GIIRS探测仪偏差特征分析和偏差订正

干涉式大气垂直探测仪（Geostationary Interferometric Infrared Sounder,简称GIIRS）是国际上第一部对地静止卫星平台上的高光谱红外大气垂直探测仪,能为对流尺度区域模式预报提供所需的高时空和高光谱分辨率的大气状态信息。本文利用高分辨率区域模式WRF及其同化系统WRFDA对GIIRS观测的偏差（观测亮温减去模拟亮温,记为O?B）分布特征进行了全景分析,结果表明:长波通道O?B偏差和标准差普遍小于中波通道,且都存在一段受污染的通道。O?B偏差的日变化和偏差与卫星天顶角的关系相对较弱,而所有筛选通道的偏差都与亮温值及卫星的扫描阵列位置有关,偏差的水平分布主要表现出“阵列偏差”的特征。2020年重新定标后,GIIRS观测数据质量比2019年有明显提高。在此基础上进一步进行了偏差订正试验,试验发现选取扫描阵列作为偏差订正的主要因子,都能有效地改进2019年和2020年筛选出的GIIRS通道的偏差,订正后O?B和O?A的系统性误差（偏差）都变小。该研究结果可为全球/区域模式中同化GIIRS长波及中波通道的辐射资料提供参考。     

Estimating the retrievability of temperature profiles from satellite infrared measurements

A method is developed to assess retrievability, namely the retrieval potential for atmospheric temperature profiles, from satellite infrared measurements in clear-sky conditions. This technique is based upon generalized linear inverse theory and empirical orthogonal function analysis. Utilizing the NCEP global temperature reanalysis data in January and July from 1999 to 2003, the retrievabilities obtained with the Atmospheric Infrared Sounder （AIRS） and the High Resolution Infrared Radiation Sounder/3 （HIRS/3） sounding channel data are derived respectively for each standard pressure level on a global scale. As an incidental result of this study, the optimum truncation number in the method of generalized linear inverse is deduced too. The results show that the retrievabilities of temperature obtained with the two datasets are similar in spatial distribution and seasonal change characteristics. As for the vertical distribution, the retrievabilities are low in the upper and lower atmosphere, and high between 400 hPa and 850 hPa. For the geographical distribution, the retrievabilities are low in the low-latitude oceanic regions and in some regions in Antarctica, and relatively high in mid-high latitudes and continental regions. Compared with the HIRS/3 data, the retrievability obtained with the AIRS data can be improved by an amount between 0.15 and 0.40.     

METOP星载干涉式超高光谱分辨率红外大气探测仪(IASI)及其产品

为了深入了解国际卫星红外大气探测仪器的新特点,以欧洲METOP卫星装载的超高光谱红外大气探测仪(IASI:Infrared Atmospheric Sounding Interferometer)为例,介绍卫星红外探测仪研发背景、仪器特征和数据产品处理,并与第1个高光谱大气红外探测仪(AIRS:Atmosphenic Inftared Sounder)作了比较分析.采用干涉分光技术的IASI为地球大气遥感提供了丰富的研究资料,可用于反演大气、海洋、云和大气成分,对我国风云卫星的高光谱红外大气垂直探测仪器的研发具有重要参考价值.     

基于ERA5的快速辐射传输模式与FY-4A成像仪观测结果的偏差分析

2016年12月发射升空的FY-4A是中国第二代静止轨道气象卫星,该星上搭载了可提供东半球近实时高分辨率卫星观测数据的扫描辐射成像仪——AGRI（Advanced Geostationary Radiation Imager）。在其观测数据应用于大气参数反演或同化前,数据偏差的定量化分析是一个必要环节。采用快速辐射传输模式RTTOV（Radiative Transfer for the TIROS Operational Vertical Sounder）,基于欧洲中期天气预报中心（ECMWF）第5个全球再分析数据产品（ERA5）对AGRI的7个红外通道（通道08—14）进行了模拟,并利用MODIS云检测产品对模拟结果进行了晴空筛选,以期得到一些对AGRI的定量应用有价值的偏差分析结果。观测-模拟（O-B）的偏差分析结果显示:海洋和陆地上,通道10（7.1 μm）存在明显大于其他红外通道的系统性偏差,这很可能来源于ERA5在对流层中层对水汽的高估。通道08为近红外短波通道,地表反射作用影响强烈,陆地上存在较大的平均偏差,而海洋上平均偏差小于0.4 K。通道14在ERA5近地层气温偏高及定标偏差的影响下,海洋存在接近1 K的平均偏差;陆地上存在2 K左右的平均偏差。其余各红外通道在海洋和陆地上的平均偏差分别在0.6和1.3 K以下。偏差影响因子分析结果显示:地表海拔高度、观测天顶角对偏差也存在一定程度的影响;海洋上偏差分布存在季节变化可能来源于再分析资料中海表温度估算的季节性误差。      

Striping Noise Analysis and Mitigation for Microwave Temperature Sounder-2 Observations

The Microwave Temperature Sounder(MWTS)-2 has a total of 13 temperature-sounding channels with the capability of observing radiance emissions from near the surface to the stratosphere. Similar to the Advanced Technology Microwave Sounder(ATMS), striping pattern noise, primarily in the cross-track direction, exists in MWTS-2 radiance observations. In this study, an algorithm based on principal component analysis(PCA) combined with ensemble empirical mode decomposition(EEMD) is described and applied to MWTS-2 brightness temperature observations. It is arguably necessary to smooth the first three principal component(PC) coefficients by removing the first four intrinsic mode functions(IMFs) using the EEMD method(denoted as PC3/IMF4). After the PC3/IMF4 noise mitigation, the striping pattern noise is effectively removed from the brightness temperature observations. The noise level in MWTS-2 observations is significantly higher than that detected in ATMS observations. In May 2014, the scanning profile of MWTS-2 was adjusted from varying-speed scanning to constantspeed scanning. The impact on striping noise levels brought on by this scan profile change is also analyzed here. The striping noise in brightness temperature observations worsened after the profile change. Regardless of the scan profile change, the striping noise mitigation method reported in this study can more or less suppress the noise levels in MWTS-2 observations.     

A New Infrared Atmospheric Sounding Interferometer Channel Selection and Assessment of Its Impact on Met Office NWP Forecasts

A new set of Infrared Atmospheric Sounding Interferometer(IASI) channels was re-selected from 314 EUMETSAT channels.In selecting channels,we calculated the impact of the individually added channel on the improvement in the analysis outputs from a one-dimensional variational analysis(1D-Var) for the Unified Model(UM) data assimilation system at the Met Office,using the channel score index(CSI) as a figure of merit.Then,200 channels were selected in order by counting each individual channel's CSI contribution.Compared with the operationally used 183 channels for the UM at the Met Office,the new set shares 149 channels,while the other 51 channels are new.Also examined is the selection from the entropy reduction method with the same 1D-Var approach.Results suggest that channel selection can be made in a more objective fashion using the proposed CSI method.This is because the most important channels can be selected across the whole IASI observation spectrum.In the experimental trial runs using the UM global assimilation system,the new channels had an overall neutral impact in terms of improvement in forecasts,as compared with results from the operational channels.However,upper-tropospheric moist biases shown in the control run with operational channels were significantly reduced in the experimental trial with the newly selected channels.The reduction of moist biases was mainly due to the additional water vapor channels,which are sensitive to the upper-tropospheric water vapor.     

An Algorithm for Detecting Ice Cloud at Different Altitudes by Combining Dual CrIS Full Spectrum Resolution CO2 Channels

Using infrared sensors to detect ice clouds in different atmospheric layers is still a challenge. The different scattering and absorption properties of longwave and shortwave infrared channels can be utilized to fulfill this purpose.In this study, the release of Suomi-NPP Cross-track Infrared Sounder(Cr IS) full spectrum resolution is used to select and pair channels from longwave(～ 15 μm) and shortwave(～4.3 μm) CO_2 absorption bands under stricter conditions, so as to better detect ice clouds. Besides, the differences of the weighting function peaks and cloud insensitive level altitudes of the paired channels are both within 50 h Pa so that the variances due to atmospheric conditions can be minimized. The training data of clear sky are determined by Visible Infrared Imaging Radiometer Suite(VIIRS) cloud mask product and used to find the linear relationship between the paired longwave and shortwave CO_2 absorption channels. From the linear relationship, the so-called cloud emission and scattering index(CESI) is derived to detect ice clouds. CESI clearly captures the center and the ice cloud features of the Super Typhoon Hato located above 415 h Pa. Moreover, the CESI distributions agree with cloud top pressure from the VIIRS in both daytime and nighttime in different atmospheric layers.     

同化FY3B/IRAS资料对GRAPES分析场影响诊断

GRAPES(Globe and Regional Assimilation and Prediction System)变分系统能够同化常规资料和非常规卫星资料,这些被同化的资料究竟对同化系统得到的分析场有何影响,目前国内外尚未见相关的研究文献。为此,首次采用基于信息熵信号自由度思想,诊断风云三号B星(FY3B)红外分光计(Infrared Atmospheric Sounder,IRAS)资料对GRAPES分析场的影响。诊断过程中,采用数值逼近方法,统计2012年12月24日18时到2013年1月22日00时共114个时次IRAS资料对GRAPES分析场影响,结果表明,IRAS中高层通道亮温资料对GRAPES分析场影响比地表通道20观测亮温的影响大,地表通道8和9观测亮温对分析场影响较大。前24个GRAPES变分同化时次每个时次IRAS通道亮温对分析场影响的贡献率分析结果显示,高层通道和H_2O通道贡献率较大。个例分析结果表明,在同化探空资料基础上加入IRAS资料后,温度和湿度增量场变化幅度较大,表明IRAS资料对分析场有降温和增湿作用。     

One-Dimensional Variational Retrieval of Temperature and Humidity Profiles from the FY4A GIIRS

A physical retrieval approach based on the one-dimensional variational(1 D-Var) algorithm is applied in this paper to simultaneously retrieve atmospheric temperature and humidity profiles under both clear-sky and partly cloudy conditions from FY-4 A GIIRS(geostationary interferometric infrared sounder) observations. Radiosonde observations from upper-air stations in China and level-2 operational products from the Chinese National Satellite Meteorological Center(NSMC)during the periods from December 2019 to January 2020(winter) and from July 2020 to August 2020(summer) are used to validate the accuracies of the retrieved temperature and humidity profiles. Comparing the 1 D-Var-retrieved profiles to radiosonde data, the accuracy of the temperature retrievals at each vertical level of the troposphere is characterized by a root mean square error(RMSE) within 2 K, except for at the bottom level of the atmosphere under clear conditions. The RMSE increases slightly for the higher atmospheric layers, owing to the lack of temperature sounding channels there.Under partly cloudy conditions, the temperature at each vertical level can be obtained, while the level-2 operational products obtain values only at altitudes above the cloud top. In addition, the accuracy of the retrieved temperature profiles is greatly improved compared with the accuracies of the operational products. For the humidity retrievals, the mean RMSEs in the troposphere in winter and summer are both within 2 g kg^–1. Moreover, the retrievals performed better compared with the ERA5 reanalysis data between 800 h Pa and 300 h Pa both in summer and winter in terms of RMSE.     

A Limb Correction Method for the Microwave Temperature Sounder 2 and Its Applications

The Microwave Temperature Sounder 2 (MWTS-2) is a cross-track radiometer that has 13 channels of sampling radiances emitted from different vertical levels of the atmosphere. Because of the varying scan angles of each field of view within a scan line, observations from the MWTS-2 are subject to strong scan-position-dependent features, i.e., the limb effect. When examining brightness temperatures (TBs), weather signals observed at every temperature-sounding channel are often concealed by scan-dependent patterns. This study, therefore, proposes a limb correction method to remove scan-dependent features so that the underlying weather signals can be uncovered. Limb-corrected TBs can be used to monitor large-scale patterns over the globe as well as extreme weather events such as typhoons. Limb-corrected TBs are also more correlated with atmospheric physical variables such as temperature and liquid water path.     

气象卫星亚像元火点面积和亮温估算方法

该文给出了利用极轨气象卫星AVHRR多个红外通道资料在不同火点强度条件下估算亚像元火点面积和温度的方法，包括利用牛顿迭代法求解不同红外波段（如:中红外3.7 μm通道和热红外11 μm通道，或短波红外1.6 μm通道和热红外11 μm通道）混合像元表达式建立的二元非线性方程组, 得到亚像元火点面积和温度；通过建立亚像元火点面积和温度与11 μm混合像元亮温增量和11 μm与12 μm混合像元亮温增量差异的关系，以查表方式获得对于相同红外波段不同通道（如11 μm和12 μm通道）亚像元火点面积和温度估算值。     

毫米/亚毫米波探测大气温度和湿度的通道选择

利用热带大气温湿廓线计算了热带地区毫米/亚毫米波段微波大气透过率权重函数。对权重函数峰值高度的分析结果显示:对流层低层的大气温度可以选择118 GHz通道的远翼频率来探测,而对高层大气温度进行探测时,选择425 GHz通道的远翼频率较为合适;在大气湿度探测方面,183 GHz通道组合适合探测对流层中层大气的湿度,高层大气湿度探测应该首先考虑380 GHz通道组合来实现。根据大气温度探测通道和大气湿度探测通道的权重函数分布,鉴于国内现有遥感仪器的制造水平,建议选择118 GHz 3个通道与425 GHz 8个通道共11个大气温度探测通道和183 GHz 3个通道与380 GHz 5个通道共8个大气湿度探测通道作为未来静止轨道微波探测的候选通道。     

A practical approach for the assimilation of cloudy infrared radiances and its evaluation using airs simulated observations

Abstract

A variational estimation procedure for the simultaneous retrieval of cloud parameters and thermodynamic profiles from infrared radiances is proposed. The method is based on a cloud emissivity model which accounts for the frequency dependence of cloud absorption and scattering and possible mixed phase situations. An effective cloud top height and emissivity are assumed. Monte Carlo experiments performed in a 1D‐var assimilation context using simulated Atmospheric Infrared Radiance Sounder (AIRS) observations from 100 channels demonstrate the substantial added value, in theory, of cloudy radiance assimilation as opposed to clear‐channel assimilation. Improved temperature and humidity retrievals are obtained for a broad layer above the cloud as well as below cloud level under partial cloud cover conditions. The impact is most pronounced in broken to overcast situations involving mid‐level clouds. In these situations, the effective cloud top height and emissivity are retrieved with estimated rms errors typically lower than 30 hPa and 3%, respectively. Expected relative errors on the retrieved effective particle size are of the order of 30–50%. The methodology is directly applicable to real hyperspectral infrared data upon inclusion, for local estimation, of the cloud parameters in the Canadian 4D‐var assimilation system.     

The study of in-orbit calibration accuracy of NOAA satellite infrared sounder and its effect on temperature profile retrievals

The calibration accuracy of High Resolution Infrared Radiation Sounder Mod. 2 (H1RS / 2) on NOAA-10 sat-ellite is analyzed in this paper. The non-linear effect in the linear calibration curve induces a deviation of 1.5 degrees (k) of brightness temperature in the tenth channel (8.3 urn, water vapor absorption) of the HIRS / 2 and the non-linear effect affects the other channels to a different extent. Based on analyzing non-linearity in two-point calibration curve, a tri-point calibration equation is given. A numerical test of effects of the linear and non-linear calibration models on the accuracy of atmospheric temperature retrievals is carried out.     

HIRS/2 亮温直接应用于数值预报的试验研究

将极轨卫星HIRS/2探测仪的高分辨率亮温资料，经过定标、定位和其它一些预处理后，不经过反演，直接和其它观测资料一起用最优内插(OI)方法进行分析，给出数值预报初始场，并用国家气象中心的准业务有限区预报模式做预报试验。对比试验结果表明，在数值预报中直接使用亮温比使用反演资料更好。