肖弘毅 韓 威 白一泓

1.中國(guó)氣象局地球系統(tǒng)數(shù)值預(yù)報(bào)中心，北京，100081

2.災(zāi)害天氣國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京，100081

3.中國(guó)氣象局中國(guó)遙感衛(wèi)星輻射測(cè)量和定標(biāo)重點(diǎn)開(kāi)放實(shí)驗(yàn)室/國(guó)家衛(wèi)星氣象中心（國(guó)家空間天氣監(jiān)測(cè)預(yù)警中心），北京，100081

4.許健民氣象衛(wèi)星創(chuàng)新中心，北京，100081

1 引言

20世紀(jì)70年代以來(lái)，已有若干代星載微波成像儀投入使用，并且在氣象、遙感、水文、地質(zhì)等多個(gè)領(lǐng)域得到廣泛的關(guān)注和研究（Moncet，et al，2011）。1978年雨云七號(hào)（Nimbus-7）衛(wèi)星和海洋衛(wèi)星（SeaSat）上分別搭載了掃描式多通道微波輻射計(jì)SMMR（Scanning Multichannel Microwave Radiometer）（Madrid，1978；Gloersen，et al，1984）；美國(guó)國(guó)防氣象衛(wèi)星系列DMSP（Defense Meteorological Satellite Program）自1987年起首次搭載了微波成像儀特別傳感器SSM/I（Special Sensor Microwave/Imager）（Hollinger，et al，1990），并自2003年起首次搭載了微波成像儀/探測(cè)器特別傳感器SSMIS（Special Sensor Microwave Imager/Sounder）（Goodrum，et al，2000）；1997年 熱 帶 測(cè) 雨 衛(wèi)星TRMM（Tropical Rainfall Measuring Mission）上搭載了微波成像儀TMI（TRMM Microwave Imager）（Kummerow，et al，2000）；2002年美國(guó)地球觀測(cè)系統(tǒng)EOS（Earth Observation System）Aqua衛(wèi)星上搭載了先進(jìn)微波掃描輻射計(jì)AMSR-E（Advanced Microwave Scanning Radiometer–Earth Observation System）（Kawanishi，et al，2003）；2003年科里奧利衛(wèi)星（Coriolis）上搭載了全極化微波輻射計(jì)WindSat（Wind Sat Polarimeteric Radiometer）（Gaiser，et al，2004）；2008年中國(guó)風(fēng)云三號(hào)系列極軌衛(wèi)星FY-3（FengYun-3）上首次搭載微波成像儀MWRI（Micro-Wave Radiation Imager）（Dong，et al，2009）；2011年中國(guó)海洋二號(hào)氣象衛(wèi)星系列HY-2（Hai Yang-2）上首次搭載微波成像儀MWRI（Microwave Radiometer Imager）（Gao，et al，2019）；2012年日本宇宙航空研究開(kāi)發(fā)機(jī)構(gòu)JAXA（Japan Aerospace Exploration Agency）發(fā)射的第一期全球水循環(huán)變化 觀 測(cè) 衛(wèi) 星GCOM-W1（Global Change Observa tion Mission–Water）上搭載了第2代先進(jìn)微波掃描輻射計(jì)AMSR2（Advanced Microwave Scanning Radiometer 2）（Oki，et al，2010；JAXA，2013）；2014年發(fā)射的全球降水測(cè)量衛(wèi)星GPM（Global Precipitation Measurement mission）上搭載了微波成 像 儀GMI（GPM Microwave Imager）（Lean，et al，2017）。迄今為止，微波成像儀已成為中外極軌氣象衛(wèi)星上搭載的主流微波探測(cè)儀器之一，并且在大氣遙感中扮演著越來(lái)越重要的角色（Kazumori，et al，2014）。

為了對(duì)近年日益受到關(guān)注的全球環(huán)境變化提供監(jiān)測(cè)，JAXA建立了“地球環(huán)境變化觀測(cè)任務(wù)（Global Change Observation Mission，GCOM）”。GCOM旨在建立可在全球范圍內(nèi)對(duì)有效地球物理參數(shù)實(shí)施10—15 a尺度連續(xù)觀測(cè)的系統(tǒng)，以闡明全球氣候變化、天氣預(yù)報(bào)、水循環(huán)等機(jī)制（JAXA，2013）。GCOM-W系列衛(wèi)星負(fù)責(zé)對(duì)水循環(huán)進(jìn)行觀測(cè)，GCOM-W1是 其 第1代 衛(wèi) 星，于2012年5月18日從種子島航天中心發(fā)射，進(jìn)入下午星軌道業(yè)務(wù)運(yùn)行。其上搭載的AMSR2為AMSR-E的繼承儀器，較后者新增了頻率為7.3 GHz的兩個(gè)通道，以緩解C波段無(wú)線電頻率間的干擾（鄒曉蕾等，2015）。作為衛(wèi)星探測(cè)器，AMSR2的觀測(cè)可以覆蓋全球，彌補(bǔ)洋面、沙漠、高原、極地等區(qū)域常規(guī)觀測(cè)資料稀少的問(wèn)題（束艾青等，2019）；同時(shí)，AMSR2和其他星載微波成像儀一樣可以發(fā)揮微波的穿透優(yōu)勢(shì)，對(duì)云雨區(qū)域地表或海表的物理信息進(jìn)行監(jiān)測(cè)（錢玲等，2019）。AMSR2可以提供有關(guān)降水、云、大氣濕度、水汽、溫度、土壤濕度、積雪分布、海面溫濕場(chǎng)、海面風(fēng)速和海表鹽度（束艾青等，2019；錢玲等，2019；俞兆文等，2017a，2017b，2018；郭黎等，2017）等大氣和地表信息，對(duì)其同化有望改善數(shù)值天氣預(yù)報(bào)的分析和預(yù)報(bào)技巧。

自2007年至今，中國(guó)氣象局?jǐn)?shù)值預(yù)報(bào)中心致力于新一代數(shù)值天氣預(yù)報(bào)系統(tǒng)“全球/區(qū)域同化與預(yù)報(bào)系統(tǒng)-全球預(yù)報(bào)系統(tǒng)（Global/Regional Assimilation and Prediction System–Global Forecast System，GRAPES_GFS）”的研發(fā)（Chen，et al，2008；薛 紀(jì) 善 等，2008；Ma，et al，2018）。為 了 匹配GRAPES全球預(yù)報(bào)模式的框架，2008年搭建了三維變分同化框架（Xue，et al，2008）。在此基礎(chǔ)上，使用三維變分同化分析方案的四維變分同化平臺(tái)于2016年在GRAPES中得以建立（王金成等，2017）。2018年，GRAPES全球四維變分同化平臺(tái)實(shí)現(xiàn)業(yè)務(wù)化（Zhang，et al，2019）。2021年，GRAPES_GFS數(shù)值天氣預(yù)報(bào)系統(tǒng)更名為CMA_GFS（Wang，et al，2021；Zhuang，et al，2021）。

CMA_GFS于2020年更新到3.0版，已有先進(jìn)微 波 探 測(cè) 單 元A/B（Advanced Microwave Sounding Unit–A/B，AMSU-A/B）（Qin，et al，2020；王金成等，2016；李剛等，2016b）、微波濕度計(jì)（Micro-Wave Humidity Sounder，MWHS）（朱利劍，2019）、微 波 溫 度 計(jì)（Micro-Wave Temperature Sounder，MWTS）（Li，et al，2014）等多臺(tái)微波探測(cè)器和高分辨率紅外探測(cè)器（High-resolution Infra-Red Sounder，HIRS）、紅外大氣探測(cè)干涉儀（Infrared Atmospheric Sounding Interferometer，IASI）（李 剛 等，2016a；鄧松等，2017）、高光譜大氣紅外探測(cè)器AIRS（Atmospheric Infra-Red Sounder）（朱 文 剛等，2013；王根等，2017）、靜止干涉紅外探測(cè)儀（Geostationary Interferometric Infrared Sounder，GIIRS）（尹若瑩等，2019；Yin，et al，2020）、紅外高光譜大氣 探 測(cè) 儀（Hyperspectral Infrared Atmospheric Sounder，HIRAS）（樊浩等，2022）等多臺(tái)紅外探測(cè)器在全球四維變分平臺(tái)中實(shí)現(xiàn)了輻射亮溫直接同化。2019年才首次實(shí)現(xiàn)風(fēng)云三號(hào)C、D兩星微波成像儀在CMA_GFS中的同化應(yīng)用（Xiao，et al，2020）。到目前為止，微波成像儀資料仍然是CMA_GFS同化較少的資料之一。相比之下，國(guó)際上，SSMIS、AMSR2、GMI等多臺(tái)微波成像儀已經(jīng)在歐洲中期天氣預(yù)報(bào)中心、日本氣象局全球模式、美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心、英國(guó)氣象局全球模式、法國(guó)氣象局全球/區(qū)域模式中實(shí)現(xiàn)了晴空區(qū)同化（Geer，et al，2018）；其中歐洲中心已經(jīng)將其全天空同化投入業(yè)務(wù)應(yīng)用，日本和美國(guó)也已在開(kāi)發(fā)相應(yīng)的全天空同化技術(shù)（Geer，et al，2018）。中國(guó)方面，Yang等（2016）和楊春等（2017）、俞兆文等（2017a，2017b，2018）、錢玲等（2019）、束艾青等（2019）先后在美國(guó)國(guó)家大 氣 研 究 中 心（National Center for Atmospheric Research，NCAR）開(kāi)發(fā)的WRF（Weather Research and Forecasting）模式的資料同化系統(tǒng)（WRFDA）（Barker，et al，2012）中實(shí)現(xiàn)了AMSR2的晴空區(qū)和云雨區(qū)資料同化，并在多個(gè)臺(tái)風(fēng)個(gè)例的分析與預(yù)報(bào)中驗(yàn)證了效果。由此可見(jiàn)，實(shí)現(xiàn)GCOM-W AMSR2的資料同化可以有效增加CMA_GFS對(duì)微波成像儀資料的使用，并有望對(duì)強(qiáng)對(duì)流和降水的分析與預(yù)報(bào)產(chǎn)生一定的正效果。

制約微波成像儀及其他微波資料同化效果的關(guān)鍵問(wèn)題包含三方面：首先，質(zhì)量控制較為粗糙，往往采用簡(jiǎn)易的殘差檢驗(yàn)等有限的幾種方案對(duì)不合理像元進(jìn)行剔除，不能從物理意義出發(fā)對(duì)污染像元進(jìn)行精準(zhǔn)屏蔽，造成保留了一些不可用的觀測(cè)噪聲，浪費(fèi)了一些包含有益信息的衛(wèi)星觀測(cè)信息（Qin，et al，2022）；其次，未能對(duì)儀器中存在的系統(tǒng)性偏差進(jìn)行合理、可靠的訂正，或者隨著儀器的老化，偏差特性已經(jīng)出現(xiàn)了漂移，但未能及時(shí)更新偏差訂正方案（Chattopadhyay，et al，2021）；第三，由于真值未知，由觀測(cè)和真值之差決定的觀測(cè)誤差難以準(zhǔn)確估計(jì)，以往的估算方法不能很好地剝離觀測(cè)誤差和模式誤差，導(dǎo)致觀測(cè)不能為同化系統(tǒng)帶來(lái)準(zhǔn)確的信息增量（Desroziers，et al，2005）。針對(duì)上述問(wèn)題，在充分調(diào)研的基礎(chǔ)上，構(gòu)造有針對(duì)性的多步驟質(zhì)量控制方案，對(duì)AMSR2由于自身儀器特性和各通道特性而受到的來(lái)自不同污染源的干擾進(jìn)行分類檢測(cè)和剔除；遴選預(yù)報(bào)因子組合對(duì)AMSR2的系統(tǒng)性偏差進(jìn)行訂正；采用基于變分同化后驗(yàn)估計(jì)統(tǒng)計(jì)觀測(cè)誤差系數(shù)，避免以往的觀測(cè)誤差取值相較實(shí)際觀測(cè)誤差過(guò)高或過(guò)低，使得求解分析場(chǎng)時(shí)對(duì)觀測(cè)場(chǎng)和背景場(chǎng)的權(quán)重分配失衡。在此基礎(chǔ)上，實(shí)現(xiàn)GCOM-W AMSR2輻射率資料在CMA全球四維變分同化系統(tǒng)中的應(yīng)用。

2 資料概述

AMSR2的通道設(shè)計(jì)與已經(jīng)在CMA_GFS 4DVar平臺(tái)上實(shí)現(xiàn)同化的風(fēng)云3號(hào)D星微波成像儀相似，但多出C波段的4個(gè)通道，共有覆蓋6.925—89.0 GHz的7個(gè)頻率，每個(gè)頻率各有垂直和水平兩種極化方式，共計(jì)14個(gè)通道（Kazumori，et al，2014）。AMSR2以55°掃描角對(duì)地球進(jìn)行圓錐掃描，每分鐘掃描40次，每次掃描提供243個(gè)掃描點(diǎn)，每個(gè)掃描點(diǎn)觀測(cè)10 km范圍的像元（89 GHz的兩個(gè)通道均設(shè)計(jì)成一對(duì)平行的觀測(cè)副端，因此該通道像元為5 km）。因此，AMSR2每通道每6 h同化時(shí)間窗共可提供約350萬(wàn)個(gè)觀測(cè)。GCOM-W AMSR2與FY-3D MWRI在通道設(shè)計(jì)上的對(duì)比見(jiàn)表1。

表1 GCOM-W AMSR2和FY-3D MWRI的通道設(shè)計(jì)Table1 Comparison of channelsdesign in GCOM-W AMSR2 and FY3D MWRI

G-Portal（Global Portal System，https://gportal.jaxa.jp/gpr/）是從JAXA地球觀測(cè)衛(wèi)星探測(cè)器獲取的產(chǎn)品。該網(wǎng)站提供GCOM-W AMSR2的兩種Level 1級(jí)數(shù)據(jù)、8種Level 2級(jí)數(shù)據(jù)及若干種Level 3級(jí)數(shù)據(jù)。文中使用的數(shù)據(jù)包括AMSR2 Level 1R級(jí)亮溫資料和Level 2級(jí)云液態(tài)水（CLW）資料，時(shí)段為2018年7月13日—8月25日。

3 研究方法

3.1 稀疏化方案

GCOM-W AMSR2大約每12 h覆蓋全球一次，產(chǎn)生近700萬(wàn)個(gè)觀測(cè)數(shù)據(jù)。為了避免資料同化中的超定問(wèn)題以及觀測(cè)場(chǎng)比重過(guò)大，需對(duì)原始資料進(jìn)行稀疏化處理。CMA全球四維變分同化業(yè)務(wù)系統(tǒng)對(duì)所有衛(wèi)星原始資料按照200 km的稀疏化半徑進(jìn)行均勻重采樣，GCOM-W AMSR2與其他儀器保持一致的稀疏化半徑設(shè)置。每6 h（1個(gè)四維變分同化時(shí)間窗）GCOM-W AMSR2有2000—3000個(gè)觀測(cè)數(shù)據(jù)進(jìn)入同化系統(tǒng)。

3.2 質(zhì)量控制方案

采取以下質(zhì)量控制方案對(duì)不合理的觀測(cè)資料進(jìn)行剔除：

（1）觀測(cè)殘差檢驗(yàn)：剔除偏差訂正后觀測(cè)亮溫偏差大于3 K的像元（俞兆文等，2017a，2018；Yang，et al，2016；楊春等，2017；劉志權(quán)等，2007）。

（2）觀測(cè)誤差檢驗(yàn)：剔除偏差訂正后觀測(cè)亮溫偏差大于3倍觀測(cè)誤差標(biāo)準(zhǔn)差的像元（楊春等，2017）。

（3）異常值檢測(cè)：剔除觀測(cè)亮溫低于70 K或高于320 K的像元（黃薇等，2013）。

（4）下墊面類型檢測(cè)：剔除下墊面為陸地、海岸線等復(fù)雜地型的像元（俞兆文等，2017a，2018；楊 春等，2017；黃薇 等，2013；Liu，et al，2012）。AMSR2的L1R級(jí)原始數(shù)據(jù)包含對(duì)應(yīng)不同通道的多個(gè)海陸掩碼數(shù)據(jù)集，其分辨率隨通道頻率的增大而增大。最后選取了6.925 GHz的海陸掩碼數(shù)據(jù)集，亦即分辨率最低的數(shù)據(jù)集進(jìn)行下墊面類型檢測(cè)，因?yàn)樵摂?shù)據(jù)集可以確保各通道中的非洋面像元全部被剔除。

（5）海冰檢測(cè)：用北極輻射與湍流相互作用研究（Arctic Radiation and Turbulence Interaction Study，ARTIST）海冰算法（ASI）（Spreen，et al，2008；蘇潔等，2013）識(shí)別并剔除下墊面為海冰的像元。

（6）降雨檢測(cè)：剔除滿足下述4個(gè)條件中任意一 個(gè) 的 像 元（郭 黎 等，2017；Bettenhausen，et al，2006；Zhao，2013；趙屹立等，2013）

式中，T為亮溫，下標(biāo)V代表垂直極化，H代表水平極化，右上角數(shù)字代表近似頻率。

（7）云檢測(cè)：同時(shí)使用兩個(gè)方案進(jìn)行云檢測(cè)，任一方案識(shí)別為云區(qū)的像元即予以剔除。

B．對(duì)每個(gè)通道，剔除CLW值大于該通道在表2中相應(yīng)閾值的像元（俞兆文等，2017a，2018；Yang，et al，2016；楊春等，2017；Kazumori，et al，2008）。

表2 AMSR2各通道的CLW檢驗(yàn)閾值Table 2 CLW check thresholds corresponding to the AMSR2 channels

（8）無(wú)線電頻率干擾檢測(cè)：AMSR2最低4個(gè)頻段的通道受到無(wú)線電頻率干擾的影響（Kazumori，et al，2016；Tian，et al，2016a；Newman，et al，2016；Lawrence，et al，2017；Wu，et al，2019）。一 般來(lái)說(shuō)，AMSR2各通道的觀測(cè)亮溫隨著通道頻率的增大而上升（Li，et al，2004；Zou，et al，2012；鄒曉蕾等，2013）。例如，23.8 GHz通道的亮溫在大部分陸地上理應(yīng)比18.7 GHz通道的亮溫要高，亦即滿足然而，當(dāng)18.7 GHz通道存在無(wú)線電頻率干擾時(shí)，低頻通道的觀測(cè)亮溫將會(huì)上升，出現(xiàn)（鄒曉蕾等，2015；Li，et al，2004）。因此，采用如下判據(jù)篩查低頻通道中的無(wú)線電頻率干擾（de Nijs，et al，2015；Zabolotskikh，et al，2015；Tian，et al，2016b；馮呈呈，2015；Wu，et al，2011）

式中，i為AMSR2的通道編號(hào)（表1），依次取1—8，TV(i)代 表第i個(gè)垂直極化通道的觀測(cè)亮溫。如果檢測(cè)到無(wú)線電頻率干擾，就屏蔽該像元的亮溫?cái)?shù)據(jù)。

（9）太陽(yáng)耀光角（太陽(yáng)耀光（sunglint）現(xiàn)象是指在微波儀器低頻通道中，當(dāng)太陽(yáng)光在地球表面上的反射光與微波探測(cè)器的入射方向相近時(shí)，地球表面將呈現(xiàn)銀鏡色，干擾正常觀測(cè)。因此，太陽(yáng)耀光角被定義為反射太陽(yáng)光與微波儀器入射方向的夾角，當(dāng)夾角接近于0，就有可能發(fā)生太陽(yáng)耀光現(xiàn)象）檢驗(yàn)：剔除6.925、7.3、10.65 GHz的6個(gè)低頻通道中太陽(yáng)耀光角小于25°的像元（Kazumori，et al，2008，2014；俞兆文等，2017a，2018；Yang，et al，2016；楊春等，2017）。

AMSR2 L1R級(jí)數(shù)據(jù)中提供了太陽(yáng)方位角、衛(wèi)星天頂角和太陽(yáng)天頂角的數(shù)據(jù)，經(jīng)過(guò)推導(dǎo)，太陽(yáng)耀光角可以寫成

式中，φ為太陽(yáng)方位角，θ1為 衛(wèi)星天頂角，θ2為太陽(yáng)天頂角。

采用（3）—（9）項(xiàng)質(zhì)量控制算法，以受云雨天氣影響最極端的89 GHz垂直極化通道為例，效果如圖1所示。為了體現(xiàn)觀測(cè)數(shù)據(jù)的分布特性，暫未采用（1）、（2）兩項(xiàng)質(zhì)量控制（后續(xù)試驗(yàn)資料均做了全部質(zhì)量控制）。從OMB（Observation-Minus-Background）的概率密度分布可以看出，質(zhì)量控制之前的觀測(cè)資料遠(yuǎn)偏離正態(tài)分布，有大量的異常值，且均值遠(yuǎn)大于0 K；質(zhì)量控制之后，OMB較大的數(shù)據(jù)已經(jīng)悉數(shù)被剔除，剩余的數(shù)據(jù)基本滿足正態(tài)分布，OMB的均值和標(biāo)準(zhǔn)差都已大幅度減小。從OMB的水平分布（圖略）可以看出，質(zhì)量控制之前的數(shù)據(jù)在陸地、海岸線、海冰等區(qū)域由于地表輻射率的影響，數(shù)據(jù)不可用，在洋面云雨區(qū)也與背景場(chǎng)有較大的偏離；經(jīng)過(guò)質(zhì)量控制，這些不合理的數(shù)據(jù)已經(jīng)絕大部分得以屏蔽，保留的數(shù)據(jù)大體上對(duì)應(yīng)于海面晴空區(qū)，有較小的OMB。由此可見(jiàn)，采用的上述質(zhì)量控制方案有針對(duì)性地遴選了合理、可用的觀測(cè)資料。

圖1 GCOM-W AMSR2（a）6.925 GHz水平極化通道和（b）89 GHz垂直極化通道在質(zhì)量控制前 （陰影柱）、后（實(shí)心柱）的OMB概率密度分布Fig.1 OMB probability density profiles（PDF）of the GCOM-W AMSR2（a）06H and（b）89V channels before（shadow bar）and after （black bar）quality control（QC）

為了驗(yàn)證太陽(yáng)耀光現(xiàn)象對(duì)觀測(cè)的影響，從6.925 GHz水平極化通道在僅采用（3）—（8）項(xiàng)質(zhì)量控制后的OMB水平分布（圖2a）可以看到，在印度洋地區(qū)存在明顯的深紅色（約3 K）的條帶未能被質(zhì)量控制方案識(shí)別，這些條帶即是受到太陽(yáng)耀光現(xiàn)象干擾產(chǎn)生的噪聲。在添加第（9）項(xiàng)算法后，質(zhì)量控制方案有效地剔除了該區(qū)域的偏差（圖2b）。

圖2 在已加入第（3）—（8）項(xiàng)質(zhì)量控制的基礎(chǔ)上，添加第（9）項(xiàng)質(zhì)量控制前（a）、后（b），6.925 GHz水平極化通道的OMB水平分布Fig.2 OMB distribution profiles of 06H channel（a）before and（b）after the ninth quality control procedure being activated after all the（3）—（8）quality control procedures have been applied

3.3 偏差訂正方案

盡管在SSMI/S（Bell，et al，2008）、TMI（Geer，et al，2010）、FY-3C MWRI（Lawrence，et al，2017；張淼等，2019；Xie，et al，2019）等多個(gè)同類儀器中觀測(cè)到升降軌偏差，但AMSR2被證實(shí)并不存在該方面的問(wèn)題（Lawrence，et al，2017；張淼等，2019）。因此，經(jīng)典的預(yù)報(bào)因子組合被用于GCOM-W AMSR2的偏差訂正。該偏差訂正方案包括氣團(tuán)屬性和掃描位置，其中，氣團(tuán)偏差訂正采用4種預(yù)報(bào)因子的組合：模式初始場(chǎng)厚度（1000—300 hPa）、模式初始場(chǎng)厚度（200—50 hPa）、模式初始場(chǎng)表面溫度，以及模式初始場(chǎng)水汽含量（Yang，et al，2016；劉志權(quán)等，2007）。2018年7月13—25日的GCOM-W AMSR2資料被用于統(tǒng)計(jì)偏差訂正系數(shù)。圖3為采用上述方案計(jì)算的偏差訂正系數(shù)對(duì)GCOM-W AMSR2觀測(cè)資料進(jìn)行偏差訂正試驗(yàn)的結(jié)果?？梢钥闯?，經(jīng)過(guò)偏差訂正，所有通道OMB的平均值和標(biāo)準(zhǔn)差都有了明顯的減小。此外，OMB的概率密度分布曲線（圖略）在偏差訂正后明顯更趨于正態(tài)分布。因此，上述偏差訂正方案能夠有效改善原始觀測(cè)資料的系統(tǒng)性偏差。

圖3 GCOM-W AMSR2 14個(gè)通道偏差訂正前（空心）、后（實(shí)心）OMB平均值（紅方框）與標(biāo)準(zhǔn)差（藍(lán)三角）Fig.3 Mean value（red square）and standard deviation（blue triangle）of OMB before（empty）and after （solild）the bias correction （BC）is applied to the 14 channels on GCOM-W AMSR2

3.4 通道選擇

經(jīng)過(guò)上述對(duì)GCOM-W AMSR2 14個(gè)通道資料質(zhì)量的初步診斷分析，確定偏差小于0.2 K且標(biāo)準(zhǔn)差小于3 K為通道選擇判據(jù)，共有10個(gè)通道滿足該判據(jù)，分別是06V、06H、07V、07H、10V、10H、19V、23V、37V、89V。

成都市經(jīng)過(guò)多年的探索與實(shí)踐，逐漸摸索出兒童課外閱讀推廣的三種模式，一是區(qū)域整體推進(jìn)，二是學(xué)校特色推進(jìn)，三是教師自發(fā)推進(jìn)。

3.5 觀測(cè)誤差統(tǒng)計(jì)

利用Desroziers提出的變分同化后驗(yàn)估計(jì)（Desroziers，et al，2005）對(duì)觀測(cè)誤差進(jìn)行估算。觀測(cè)誤差可以寫成OMB和OMA（Observation-Minus-Analysis）的幾何平均數(shù)。利用上述方法對(duì)2018年7月13日—8月25日的GCOM-W AMSR2觀 測(cè) 資料統(tǒng)計(jì)觀測(cè)誤差，所得結(jié)果如表3所示。

表3 GCOM-W AMSR2參加同化的10個(gè)通道的觀測(cè)誤差Table 3 Observation errors corresponding to the ten channels of GCOM-W AMSR2 assimilated in the present study

3.6 試驗(yàn)設(shè)置

為了檢驗(yàn)GCOM-W AMSR2輻射亮溫資料同化對(duì)全球模式分析和預(yù)報(bào)的影響，在CMA_GFS 3.0版中開(kāi)展了試驗(yàn)。CMA_GFS 3.0版的基本設(shè)置為：水平分辨率0.5°，垂直87層非均勻分層，模式層頂0.1 hPa，時(shí)間步長(zhǎng)450 s，三維參考大氣，預(yù)估-修正的SISL時(shí)間積分算法，Garcia非地形重力波，Goff Gratch飽和水汽壓計(jì)算公式，半隱式系數(shù)設(shè)置為0.55。試驗(yàn)所選取的主要物理過(guò)程參數(shù)化方案（Ma，et al，2018）包括：CMA雙參數(shù)微物理方案（劉奇俊等，2003）、預(yù)報(bào)云方案（Tiedtke，1993）、RRTMG長(zhǎng) 短 波 輻 射 方 案（Pincus，et al，2003；Morcrette，et al，2008）、CoLM（Common Land Model）陸面過(guò)程方案（Dai，et al，2003）、NMRF邊界層參數(shù)化方案（Hong，et al，1996）、NSAS（New Simplified Arakawa-Schubert）積云對(duì)流參數(shù)化方案（劉奇俊等，2003；Arakawa，et al，1974；Pan，et al，1995）、重力波拖曳方案。

本研究共設(shè)計(jì)2組批量試驗(yàn)：第1組為對(duì)照試驗(yàn)（縮寫為CTRL），同化所有常規(guī)資料（包括云導(dǎo)風(fēng)資料、海面/地面站觀測(cè)資料、探空資料、掩星資料等）和業(yè)務(wù)上默認(rèn)開(kāi)啟的衛(wèi)星資料；第2組為敏感性試驗(yàn)（縮寫為AMSR2），在對(duì)照試驗(yàn)的基礎(chǔ)上，新增同化GCOM-W AMSR2輻射亮溫資料，并對(duì)其采用上文的通道選擇方案、稀疏化方案、質(zhì)量控制方案、偏差訂正方案、觀測(cè)誤差方案。2組試驗(yàn)均開(kāi)啟所有衛(wèi)星儀器的動(dòng)態(tài)偏差訂正。試驗(yàn)時(shí)段均為2018年7月25日至8月24日（1個(gè)月）。

每組批量試驗(yàn)均包含兩個(gè)模塊。在第1個(gè)模塊中，試驗(yàn)會(huì)在每天的03、09、15和21時(shí)（世界時(shí)，下同）分別進(jìn)行4次同化，并分別得到當(dāng)日06、12、18時(shí)及次日00時(shí)4個(gè)時(shí)次的分析場(chǎng)。第1次同化（2018年7月25日03時(shí)）采用的背景場(chǎng)由當(dāng)日00時(shí)的全球再分析數(shù)據(jù)通過(guò)3 h預(yù)報(bào)得到。其他每次同化均循環(huán)采用前一次同化經(jīng)過(guò)6 h預(yù)報(bào)之后生成的背景場(chǎng)。每個(gè)分析時(shí)次均對(duì)前、后各3 h的資料進(jìn)行同化。在第2個(gè)模塊中，從每天的00和12時(shí)兩個(gè)時(shí)間點(diǎn)開(kāi)始，分別進(jìn)行為期240 h的預(yù)報(bào)。

4 結(jié)果分析

由于AMSR2對(duì)水汽較為敏感，它的同化有望為濕度場(chǎng)增加有用信息。因此，以歐洲中期天氣預(yù)報(bào)中心（ECMWF）的第5代大氣再分析資料（ERA5）為參考，以比濕為濕度場(chǎng)檢驗(yàn)要素，對(duì)濕度場(chǎng)的均方根誤差（RMSE）進(jìn)行計(jì)算，說(shuō)明GCOM-W AMSR2輻射亮溫資料同化對(duì)分析的影響，結(jié)果如圖4所示。

圖4 批量試驗(yàn)濕度場(chǎng)均方根誤差的改進(jìn)百分比（黑線：北半球；藍(lán)線：南半球；紅線：赤道地區(qū)。黑色細(xì)橫線陰影區(qū)：北半球濕度分析改進(jìn)廓線的誤差范圍；藍(lán)色細(xì)豎線陰影區(qū)：南半球濕度分析改進(jìn)廓線的誤差范圍；紅色圓點(diǎn)陰影區(qū)：赤道地區(qū)濕度分析改進(jìn)廓線的誤差范圍）Fig.4 The root-mean-square-error （RMSE）reduction rate of humidity analysis results of batch experiments（the black，blue and red lines correspond to the Northern Hemisphere，the Southern Hemisphere and the equatorial region，respectively；the shadow area filled by black horizontal thin lines，blue vertical thin lines and red dots correspond to error ranges in the Northern Hemisphere，the Southern Hemisphere and the equatorial region，respectively）

總體而言，同化AMSR2資料的效果南半球優(yōu)于北半球，北半球優(yōu)于赤道地區(qū)，這可以歸因于南半球常規(guī)觀測(cè)資料稀少，AMSR2衛(wèi)星資料的同化可以有效地填補(bǔ)該區(qū)域的觀測(cè)資料，而赤道地區(qū)由于存在長(zhǎng)期的大面積降水云，在當(dāng)前的質(zhì)量控制方案下可以進(jìn)入同化系統(tǒng)的AMSR2資料比其他兩個(gè)區(qū)域少。對(duì)于南半球，同化AMSR2資料后，900—700 hPa和400 hPa以上高度的濕度分析場(chǎng)有顯著改進(jìn)，最大改進(jìn)在150 hPa處，達(dá)到22%，而對(duì)于低層的改進(jìn)與此前同化FY-3D MWRI的結(jié)果一致（Xiao，et al，2020），這與微波成像儀的權(quán)重函數(shù)主要對(duì)低層大氣較為敏感有關(guān)。對(duì)于赤道地區(qū)，同化AMSR2資料對(duì)850 hPa高度以下濕度分析場(chǎng)的改進(jìn)要明顯優(yōu)于北半球和南半球，這可以從AMSR2對(duì)低層大氣水汽較為敏感得到解釋。

為了進(jìn)一步揭示同化AMSR2資料對(duì)分析場(chǎng)質(zhì)量的改進(jìn)，分別繪制了對(duì)照試驗(yàn)和敏感試驗(yàn)的濕度分析場(chǎng)偏差的垂直廓線（圖5）?？梢钥吹剑珻MA_GFS對(duì)照試驗(yàn)的濕度分析場(chǎng)在低層大氣（低于750 hPa）存在普遍的負(fù)偏差，而在中高層大氣（750—300 hPa）存在一定的正偏差。同化AMSR2輻射率資料后，不同地區(qū)各層大氣的分析場(chǎng)偏差從整體上得到了一致的改進(jìn)，無(wú)論是北半球、南半球或是熱帶地區(qū)，偏差的絕對(duì)值都在減小。只有南半球和赤道地區(qū)的近地面層（925 hPa附近）偏差略有增大。

圖5 （a）北半球、（b）南半球和（c）赤道地區(qū)的濕度分析場(chǎng)偏差垂直廓線（黑線：CTRL，紅線：AMSR2）Fig.5 Vertical profiles of biases of humidity analysis fields in（a）the Northern Hemisphere，（b）the Southern Hemisphere and （c）the equatorial region （the black and red lines correspond to CTRL and AMSR2，respectively）

AMSR2的水汽敏感特性有望在降水預(yù)報(bào)中發(fā)揮優(yōu)勢(shì)。同化GCOM-W AMSR2對(duì)降水預(yù)報(bào)的影響如圖6所示。站點(diǎn)降水檢驗(yàn)結(jié)果表明，各種程度的72 h降水預(yù)報(bào)，同化GCOM-W AMSR2輻射亮溫資料后的TS均明顯高于對(duì)照試驗(yàn)。這說(shuō)明GCOM-W AMSR2資料的應(yīng)用對(duì)降水的中期預(yù)報(bào)技巧有較好的提升，這與已有報(bào)道的AMSR2同化對(duì)諸如臺(tái)風(fēng)（束艾青等，2019；錢玲等，2019；俞兆文等，2017a，2018；Yang，et al，2016；楊春等，2017）、暴雨（俞兆文等，2017b）等強(qiáng)對(duì)流天氣事件的預(yù)報(bào)有正效果的結(jié)論是一致的。為了更好地突出AMSR2同化在強(qiáng)降水預(yù)報(bào)方面的表現(xiàn)，選取2018年8月19日受臺(tái)風(fēng)“溫比亞”影響，在山東省全境發(fā)生，直接導(dǎo)致壽光洪災(zāi)的極端暴雨個(gè)例（鄭怡等，2019）進(jìn)行TS評(píng)分檢驗(yàn)。結(jié)果表明，敏感試驗(yàn)在此次個(gè)例中對(duì)各種程度的72 h降水預(yù)報(bào)評(píng)分均較對(duì)照試驗(yàn)有大幅度提升，特別是對(duì)50 mm以上強(qiáng)降水的預(yù)報(bào)水平有顯著提升。兩組試驗(yàn)對(duì)100 mm以上暴雨的預(yù)報(bào)技巧均仍然有提升空間，然而由于受到空間分辨率的制約，這在全球模式中是一個(gè)普遍存在的問(wèn)題（宮宇等，2018），有待進(jìn)一步的研究和考察。

圖6 （a）全球批量試驗(yàn)和（b）2018年“8.19”山東大范圍暴雨個(gè)例試驗(yàn)的72 h站點(diǎn)降水檢驗(yàn)TS評(píng)分Fig.6 TSof point precipitation verification at 72 h from（a）global forecast batch experiments and（b）single case experiment of "8.19" Shandong rainstorm in 2018

從敏感試驗(yàn)相對(duì)對(duì)照試驗(yàn)的預(yù)報(bào)評(píng)分卡結(jié)果（圖7a）來(lái)看，同化GCOM-W AMSR2輻射亮溫資料對(duì)全球數(shù)值天氣預(yù)報(bào)存在總體上的正效果。北半球的效果大體為中性，而南半球的正效果最為明顯，這可以從北半球的觀測(cè)資料已經(jīng)非常充分，而南半球以海洋為主常規(guī)觀測(cè)資料匱乏中得到解釋。赤道地區(qū)同樣表現(xiàn)出明顯的正效果，特別是u風(fēng)場(chǎng)存在明顯改進(jìn)，與已有關(guān)于AMSR2同化可改進(jìn)海面風(fēng)速預(yù)報(bào)的研究結(jié)果一致（錢玲等，2019；楊春等，2017），只有高度場(chǎng)短期預(yù)報(bào)的RMS評(píng)分和850 hPa溫度場(chǎng)的中短期預(yù)報(bào)表現(xiàn)出一些負(fù)效果。此外，東亞地區(qū)也存在一些中性偏正的效果，并且仍然集中在u風(fēng)場(chǎng)（錢玲等，2019；楊春等，2017）。敏感試驗(yàn)在全球范圍內(nèi)各個(gè)高度多個(gè)物理量的距平相關(guān)、平均偏差、均方根誤差和標(biāo)準(zhǔn)差（圖7b、c）也相對(duì)對(duì)照試驗(yàn)有所優(yōu)化，與已有研究注意到的AMSR2同化對(duì)海平面氣壓場(chǎng)的改進(jìn)一致（錢玲等，2019；楊春等，2017）。綜上所述，GCOM-W AMSR2輻射亮溫資料的同化，可以有效地提高CMA全球預(yù)報(bào)系統(tǒng)的分析和預(yù)報(bào)技巧。

圖7 批量試驗(yàn)的預(yù)報(bào)結(jié)果（a.綜合評(píng)分卡，b.1000 hPa濕度場(chǎng)的平均偏差，c.1000 hPa高度場(chǎng)的ACC）Fig.7 Forecasting results of batch experiments（a.forecast score card，b.mean bias of specific humidity at 1000 hPa，c.anomaly correlation（ACC）of geopotential field at 1000 hPa）

5 結(jié)論與討論

針 對(duì)GCOM-W AMSR2在CMA_GFS 4DVar平臺(tái)中的直接同化應(yīng)用開(kāi)展了系統(tǒng)的研究和試驗(yàn)：

（1）在已有研究的基礎(chǔ)上，整合了一套針對(duì)GCOM-W AMSR2的質(zhì)量控制算法，可對(duì)異常觀測(cè)（觀測(cè)誤差檢驗(yàn)、觀測(cè)殘差檢驗(yàn)、異常值檢驗(yàn)）、地面污染（陸地檢驗(yàn)、海冰檢驗(yàn)）、天氣影響（降雨檢驗(yàn)、云檢測(cè)）、人為因素（無(wú)線電干擾）和AMSR2特有的光學(xué)問(wèn)題（太陽(yáng)耀光檢驗(yàn)）進(jìn)行剔除，并用OMB對(duì)質(zhì)量控制效果進(jìn)行診斷。結(jié)果表明，質(zhì)量控制算法可以有效排除存在問(wèn)題的像元，并保留合理亮溫進(jìn)入同化系統(tǒng)。

（2）根據(jù)GCOM-W AMSR2的通道特性，選擇經(jīng)典預(yù)報(bào)因子組合對(duì)其進(jìn)行偏差訂正，并用OMB對(duì)偏差訂正效果進(jìn)行驗(yàn)證。結(jié)果表明，偏差訂正后AMSR2各通道觀測(cè)的平均偏差都顯著減小至接近于0，同時(shí)偏差的標(biāo)準(zhǔn)差也明顯下降。此外，各通道觀測(cè)的OMB分布也更接近理論預(yù)期的正態(tài)分布。因此，偏差訂正算法可以有效地校正觀測(cè)中存在的、無(wú)法通過(guò)質(zhì)量控制方案移除的系統(tǒng)性誤差。

（3）在質(zhì)量控制和偏差訂正對(duì)AMSR2各通道性能進(jìn)行診斷的基礎(chǔ)上，從總共14個(gè)通道中選取10個(gè)通道進(jìn)入后續(xù)研究，對(duì)其設(shè)計(jì)稀疏化前處理方案，用變分同化后驗(yàn)估計(jì)統(tǒng)計(jì)其觀測(cè)誤差，并在CMA_GFS 4DVar平臺(tái)上搭建GCOM-W AMSR2的同化框架。實(shí)現(xiàn)了AMSR2輻射亮溫的直接同化。

（4）開(kāi)展批量試驗(yàn)，對(duì)CMA_GFS 4DVar同化GCOM-W AMSR2對(duì)分析和預(yù)報(bào)的影響進(jìn)行檢驗(yàn)。結(jié)果表明，GCOM-W AMSR2的水汽敏感特性可以既體現(xiàn)在CMA_GFS對(duì)濕度分析場(chǎng)刻畫的改進(jìn)上，也反映在各種程度降水的中期預(yù)報(bào)技巧的提升中，同時(shí)也反映在富含水汽的赤道地區(qū)預(yù)報(bào)評(píng)分的提高上。此外，同化GCOM-W AMSR2發(fā)揮了衛(wèi)星觀測(cè)的優(yōu)勢(shì)，有效改善了南半球地區(qū)的預(yù)報(bào)水平。

鑒于CMA加載RTTOV輻射傳輸模式時(shí)尚未接入RTTOV_SCATT模塊，無(wú)法正確模擬受散射影響的觀測(cè)亮溫，因此暫時(shí)未能實(shí)現(xiàn)云雨區(qū)的衛(wèi)星微波觀測(cè)資料同化。另外，受限于地表發(fā)射率的計(jì)算精度和地表溫度的數(shù)據(jù)誤差，陸面地區(qū)特別是沙漠、海冰等復(fù)雜下墊面區(qū)域，對(duì)低層大氣敏感的微波通道由于受到地面信息的污染，其資料同化在學(xué)界仍然是國(guó)際性的挑戰(zhàn)。此外，由于微波成像儀通道數(shù)量較少，且預(yù)處理步驟對(duì)觀測(cè)資料進(jìn)行了稀疏化，同時(shí)質(zhì)量控制方案大量剔除了最能發(fā)揮其水汽敏感特性的熱帶地區(qū)數(shù)據(jù)，種種因素都導(dǎo)致能夠進(jìn)入同化系統(tǒng)的AMSR2觀測(cè)數(shù)據(jù)相對(duì)有限。隨著CMA經(jīng)、緯度網(wǎng)格分辨率的不斷提升，稀疏化半徑有望在將來(lái)相應(yīng)地減小。目前，復(fù)雜下墊面上低層微波通道的資料同化研究已經(jīng)在CMA_GFS中開(kāi)展，未來(lái)也計(jì)劃嘗試云雨區(qū)的微波通道同化應(yīng)用，從而發(fā)揮微波探測(cè)高穿透性的特點(diǎn)。在GCOM-W AMSR2等星載微波儀器實(shí)現(xiàn)全地表、全天候同化的基礎(chǔ)上，有望進(jìn)一步提升衛(wèi)星資料的使用量與全球覆蓋率，更充分地挖掘和驗(yàn)證GCOM-W AMSR2在改進(jìn)分析和預(yù)報(bào)方面的應(yīng)用潛力。