韋曉澄 , 孫逢林
(1. 中國(guó)氣象科學(xué)研究院, 北京 100081;2. 中國(guó)氣象局中國(guó)遙感衛(wèi)星輻射測(cè)量和定標(biāo)重點(diǎn)開(kāi)放實(shí)驗(yàn)室/國(guó)家衛(wèi)星氣象中心(國(guó)家空間天氣監(jiān)測(cè)預(yù)警中心), 北京 100081;3. 許健民氣象衛(wèi)星創(chuàng)新中心,北京 100081)
準(zhǔn)確、及時(shí)地獲取云中的溫度、濕度結(jié)構(gòu)、降水參數(shù)信息及云微物理參數(shù)(如液水路徑、冰晶含量與有效粒子半徑)是天氣預(yù)報(bào)分析工作中極為重要的一環(huán),能為強(qiáng)對(duì)流、暴雨和強(qiáng)臺(tái)風(fēng)等極端天氣的預(yù)警提供重要的參考信息。除此之外,云微物理參數(shù)和光學(xué)參數(shù)對(duì)地表和地球的輻射收支、大氣的加熱率和冷卻率的垂直分布等造成直接影響[1],從而影響整個(gè)地球大氣系統(tǒng)的氣候態(tài)[2-3]。
目前星載儀器中,對(duì)云進(jìn)行觀測(cè)的手段主要分為紅外遙感和微波遙感兩種。通過(guò)對(duì)紅外成像通道數(shù)據(jù)的分析,能獲取云頂?shù)奈锢硖卣?,且其能長(zhǎng)時(shí)間對(duì)某一地區(qū)進(jìn)行連續(xù)觀測(cè)[4-5]。極軌衛(wèi)星微波載荷因其對(duì)云的穿透能力強(qiáng),能獲取云中溫度、濕度等氣象要素與各項(xiàng)云微物理參數(shù),對(duì)數(shù)值預(yù)報(bào)效果改善的貢獻(xiàn)居全部觀測(cè)的首位,在強(qiáng)對(duì)流、暴雨和強(qiáng)臺(tái)風(fēng)等極端天氣預(yù)警預(yù)報(bào)業(yè)務(wù)工作中發(fā)揮著重要的作用[6]。但是,傳統(tǒng)的極軌衛(wèi)星微波載荷具有再訪時(shí)間長(zhǎng),觀測(cè)覆蓋范圍小的特點(diǎn),當(dāng)前的極軌雙星觀測(cè)體系只能保障6 小時(shí)的觀測(cè)周期[7],這無(wú)疑會(huì)大大降低極端天氣預(yù)警的時(shí)效性。因此,實(shí)現(xiàn)靜止軌道微波觀測(cè)一直都是近年來(lái)的主要研究目標(biāo)之一,如美國(guó)的靜止軌道微波觀測(cè)平 臺(tái)(Geostationary Microwave Observatory,GEM)[8]以及歐洲的微波大氣探測(cè)靜止觀測(cè)平臺(tái)(Geostationary Observatory for Microwave Atmospheric Soundings,GOMAS)。這些未來(lái)的靜止軌道微波探測(cè)器是直接瞄準(zhǔn)在短臨預(yù)報(bào)和臺(tái)風(fēng)等災(zāi)害性天氣系統(tǒng)的監(jiān)測(cè)的應(yīng)用。
然而,在靜止軌道星載儀器的對(duì)地觀測(cè)過(guò)程中,衛(wèi)星觀測(cè)到的云位置在不同視角下會(huì)有所偏差,即會(huì)產(chǎn)生視差問(wèn)題(下文記為視差)。影響星載探測(cè)儀器視差的因子有很多,包括衛(wèi)星分別到星下點(diǎn)與到探測(cè)點(diǎn)之間連線形成的張角、儀器的空間分辨率以及觀測(cè)點(diǎn)處的云高等。考慮到微波類儀器空間分辨率相對(duì)較大,國(guó)內(nèi)外視場(chǎng)偏差的研究主要針對(duì)紅外探測(cè)儀器,對(duì)微波儀器的研究相對(duì)偏少,因此,為了后續(xù)更好地應(yīng)用靜止微波探測(cè)資料,我們針對(duì)微波儀器的視差問(wèn)題開(kāi)展相關(guān)研究。
現(xiàn)在雖無(wú)投入業(yè)務(wù)使用的靜止軌道微波儀器,但是對(duì)于未來(lái)靜止軌道微波探測(cè)器在天氣方面的應(yīng)用,視差問(wèn)題是否會(huì)造成資料的應(yīng)用影響是本文重點(diǎn)探討的科學(xué)問(wèn)題。但由于沒(méi)有實(shí)際觀測(cè)數(shù)據(jù),本研究將基于我國(guó)風(fēng)云四號(hào)氣象衛(wèi)星(FY-4A)的相關(guān)數(shù)據(jù)和仿真模擬資料來(lái)展開(kāi)探討,主要從視差訂正方法、衛(wèi)星空間分辨率對(duì)觀測(cè)視差的影響、衛(wèi)星天頂角對(duì)視差的影響進(jìn)行理論研究,并結(jié)合實(shí)際觀測(cè)個(gè)例來(lái)分析靜止軌道微波探測(cè)器的視差特征。
基于中國(guó)風(fēng)云四號(hào)(FY-4A)A 星觀測(cè)的經(jīng)緯度數(shù)據(jù)進(jìn)行未來(lái)攜帶微波儀器的靜止軌道衛(wèi)星觀測(cè)位置的模擬,并且基于當(dāng)前極軌衛(wèi)星上微波探測(cè)儀器空間分辨率參數(shù),針對(duì)未來(lái)靜軌微波儀器設(shè)計(jì)了幾種空間分辨率進(jìn)行模擬試驗(yàn)。其中,風(fēng)云四號(hào)氣象衛(wèi)星A 星是由上海航天技術(shù)研究院研發(fā),于2016年12月11日在我國(guó)西昌衛(wèi)星發(fā)射中心發(fā)射的新一代地球靜止軌道定量遙感氣象衛(wèi)星。它的軌道高度為36 000 km,星下點(diǎn)經(jīng)度為104.7 °E,上面搭載了多通道掃描成像輻射計(jì)、干涉式大氣垂直探測(cè)儀、閃電成像儀和空間天氣檢測(cè)儀器包[9]。比起我國(guó)第一代靜止氣象衛(wèi)星風(fēng)云二號(hào)(FY-2),F(xiàn)Y-4A 不僅由自旋穩(wěn)定式觀測(cè)模式改成了三軸穩(wěn)定的控制方案,通道數(shù)也由5個(gè)通道增加至14 個(gè)通道。其中,包括6 個(gè)空間分辨率為0.5~1 km 的可見(jiàn)光通道[10],空間分辨率為2 km 的中短波紅外通道以及8個(gè)空間分辨率為4 km 的長(zhǎng)波紅外通道[11]。此外,最明顯的進(jìn)步是FY-4A 可以每15 分鐘獲得一張全圓盤(pán)圖(FY-2 號(hào)是1 小時(shí)一張圖),每3 分鐘就對(duì)臺(tái)風(fēng)尺度的區(qū)域?qū)崿F(xiàn)一次觀測(cè),為極端天氣的預(yù)警和監(jiān)測(cè)保駕護(hù)航。
現(xiàn)有的極軌氣象衛(wèi)星星載微波探測(cè)器的空間分辨率各不相同。我國(guó)風(fēng)云三號(hào)系列衛(wèi)星上搭載的微波溫度探測(cè)儀(MWTS)最低點(diǎn)分辨率約為50 km,微波濕度探測(cè)儀(MWHS)最低點(diǎn)分辨率為15 km[12]。美國(guó)NOAA-KLM 系列極軌衛(wèi)星的先進(jìn)微波探測(cè)器(AMSU)中,AMSU-A 星下點(diǎn)空間分辨率約為48 km,AMSU-B 星下點(diǎn)空間分辨率為16 km[13-14]。2011 年10 月28 日美國(guó)發(fā)射成功的NPP 衛(wèi)星上搭載的新一代微波探測(cè)儀(ATMS)不同通道的星下點(diǎn)分辨率為16~75 km 不等。又因?yàn)殪o止軌道微波探測(cè)器的觀測(cè)亮溫模擬過(guò)程中所使用的ERA-5 再分析場(chǎng)數(shù)據(jù)的時(shí)間分辨率為1 小時(shí),空間分辨率為0.25 °×0.25 °,所以我們理論敏感性試驗(yàn)所設(shè)定的衛(wèi)星空間分辨率選取近似的25 km(約0.25 °)、50 km和75 km。
傳統(tǒng)的極軌輻射計(jì)采用的是逐行機(jī)械掃描(靜止軌道微波輻射計(jì)的二維波束掃描問(wèn)題仍亟待解決),為了保障分辨率能最大限度地統(tǒng)一,在計(jì)算過(guò)程中,我們以星下點(diǎn)所在的位置進(jìn)行緯向檢索,當(dāng)所要求的的空間分辨率為a km 時(shí),盡量保證星下點(diǎn)所在的列上的點(diǎn)與其緯向相鄰點(diǎn)間的距離盡量接近a,及每一行上的點(diǎn)與其經(jīng)向鄰點(diǎn)間的距離約為a,最大不超過(guò)2a。
因?yàn)槟壳安](méi)有在軌的已投入業(yè)務(wù)使用的微波輻射計(jì)。本文基于FY-4A 的紅外亮溫?cái)?shù)據(jù),使用美國(guó)衛(wèi)星資料同化聯(lián)合中心(The US Joint Center for Satellite Data Assimilation, JCSDA)開(kāi)發(fā)的快速輻射傳輸模式(Community Radiative Transfer Mode, CRTM)[15-16]模擬得到靜止軌道微波儀器觀測(cè)亮溫,其中使用歐洲中期天氣預(yù)報(bào)再分析中心(European Centre for Medium-Range Weather Forecasts Reanalysis, ERA-5)的數(shù)據(jù)作為正演大氣環(huán)境場(chǎng)。數(shù)值模擬的靜止軌道微波探測(cè)器通道頻率與美國(guó)S-NPP 衛(wèi)星上搭載的先進(jìn)技術(shù)微波探測(cè)儀(Advanced Microwave Technology Sensor, ATMS)一致,各個(gè)通道的頻率與空間分辨率如表1所示。在模擬研究中,使用了FY-4A成像儀的業(yè)務(wù)云產(chǎn)品數(shù)據(jù)(基于發(fā)展成熟的風(fēng)云四號(hào)成像儀科學(xué)算法[17-18]計(jì)算而來(lái)),包括云掩膜(Cloud Mask)數(shù)據(jù)以及云頂高度數(shù)據(jù)(Cloud Top Height)。
表1 ATMS 儀器通道設(shè)置情況介紹
云高數(shù)據(jù)則是來(lái)自FY-4A 的二級(jí)云產(chǎn)品資料[18],并由此匹配得到微波儀器每個(gè)視場(chǎng)的云高。圖1 中為基于再分析資料和CRTM 模擬得到的2018 年8 月17 日00 點(diǎn)(UTC)的靜止軌道微波探測(cè)器觀測(cè)亮溫結(jié)果,分別展示了通道頻率為31.40 GHZ(低頻窗區(qū)通道)、54.40 GHz(溫度探測(cè)通道)、88.20 GHz(高頻窗區(qū)通道)和183.31±7.00 GHZ(濕度探測(cè)通道)的亮溫模擬結(jié)果。
圖1 基于CRTM模式和ERA-5資料模擬得到的2018年8月17日00時(shí)(UTC)的靜止軌道微波探測(cè)器觀測(cè)亮溫
靜止衛(wèi)星定點(diǎn)在赤道上空,除了對(duì)星下點(diǎn)的觀測(cè),衛(wèi)星都是沿斜路徑(有一定的衛(wèi)星天頂角)觀測(cè)地球大氣系統(tǒng),定位到地球表面(圖2,定位到B 點(diǎn))。但云漂浮在大氣中,B 位置并不是云在地表投影的真實(shí)位置,B點(diǎn)和D點(diǎn)的距離認(rèn)為是視場(chǎng)偏差。一般來(lái)說(shuō),云頂高度越高,距離星下點(diǎn)越遠(yuǎn),衛(wèi)星觀測(cè)云的位置偏差也就越大[19]。在研究臺(tái)風(fēng)時(shí),尤其是多源資料時(shí)空匹配時(shí),進(jìn)行云位置的偏差修正是數(shù)據(jù)處理過(guò)程中的重要步驟之一。
研究中使用的每個(gè)視場(chǎng)的位置訂正方法與靜止軌道紅外云圖的位置訂正方法類似。如圖2 所示,云頂高度為h,地球半徑為r,衛(wèi)星與地球球心的距離為R,星下點(diǎn)位置為A,觀測(cè)所得云的位置為B,訂正后云的位置,即真云云頂位置為C,真云云下點(diǎn)位置為D。
圖2 視差訂正前后云位置發(fā)生變化的示意圖 B為衛(wèi)星原本測(cè)量的點(diǎn),但因?yàn)锽處云的干擾,測(cè)得的值為C處的值。
則由文獻(xiàn)[19]可得:
公式(4)中(當(dāng)lonC≥lonB時(shí)取“-”,反之取“+”)。下文中所用的云位置偏差即為該像素點(diǎn)上的云在訂正前后的位置之差。在本次研究中,取風(fēng)云四號(hào)A星的星下點(diǎn)位置,即A點(diǎn)坐標(biāo)(104.7,0)。
據(jù)趙敏等[20]的統(tǒng)計(jì),2000—2018 年?yáng)|亞地區(qū)云頂高度主要分布在1~11 km 之間,因此,本研究中分別選用1 km、5 km、10 km 和15 km 四類云頂高度進(jìn)行理論分析,來(lái)探討微波類儀器是否需要考慮視差、并且進(jìn)一步分析衛(wèi)星空間分辨率和云頂高度對(duì)微波視差的影響。圖3~圖5中以2018年8 月17 日00 時(shí)(UTC)的亮溫為基礎(chǔ),分別展示了通 道 頻率 為57.29 GHz 、88.20 GHz 和183.31±7.00 GHz 中視差訂正前后每個(gè)視場(chǎng)內(nèi)的亮溫差。結(jié)果表明:在不同通道中,云頂高度較高時(shí),視差問(wèn)題越發(fā)明顯,隨著儀器空間分辨率的降低,視差問(wèn)題產(chǎn)生的影響也逐漸減小。這是因?yàn)楫?dāng)觀測(cè)點(diǎn)密度降低時(shí),無(wú)法準(zhǔn)確定位每個(gè)觀測(cè)點(diǎn)處的變化情況,無(wú)法準(zhǔn)確判斷云高產(chǎn)生的視差,從而模糊該處的觀測(cè)結(jié)果。在實(shí)際的例子中,視差問(wèn)題產(chǎn)生的影響可能因不同地方的云高而產(chǎn)生差異。
圖3 不同云頂高度及不同分辨率條件下,通道10(57.29)中每個(gè)像元在視差訂正前后的亮溫變化第一列為2018年8月17日00時(shí)的模擬亮溫?cái)?shù)據(jù)在FY-4A所反演得到的云高條件下,進(jìn)行視差訂正前后的亮溫差。其余四列分別是同時(shí)刻的模擬亮溫?cái)?shù)據(jù)在模擬云高的條件下,進(jìn)行視差訂正前后的亮溫差,模擬云高分別為1 km、5 km、10 km和15 km。
圖5 同圖3,但為通道18
另外,我們還進(jìn)行了在不同通道及不同空間分辨率的條件下,反演得到的模擬亮溫的偏差分析。圖6~圖8 分別展示了空間分辨率為25 km時(shí),通道10(57.29 GHz)、通道16(88.20 GHz)和通道18(183.31±7.00 GHz)的模擬亮溫結(jié)果,濕度通道和高頻窗區(qū)通道中,視差訂正前后的亮溫偏差可達(dá)4 K,再次證明微波視場(chǎng)內(nèi)的視差問(wèn)題不容忽視。相比起其他通道,溫度探測(cè)通道的視差問(wèn)題尤為明顯,云高變高時(shí),視差訂正前后亮溫變化也更為劇烈,且訂正后的亮溫幾乎都低于訂正前的亮溫。
圖6 空間分辨率為25 km時(shí),通道10(57.29 GHz)中每個(gè)像元在視差訂正前后亮溫變化的分布
圖8 空間分辨率為25 km時(shí),通道18(183.31±7.00 GHz)中每個(gè)像元在視差訂正前后亮溫變化的分布
圖4 同圖3,但為通道16
圖7 空間分辨率為25 km時(shí),通道16(88.20 GHz)中每個(gè)像元在視差訂正前后亮溫變化的分布
總的來(lái)說(shuō),盡管微波儀器的空間分辨率較低,但仍然存在視差問(wèn)題,且在云高較高的情況下,視差的存在導(dǎo)致的各像元處的亮溫偏差不小,無(wú)法忽略,這也與以往的認(rèn)知不同。請(qǐng)注意不同通道對(duì)云層的穿透能力不同,視差問(wèn)題的大小也不同。
靜止軌道微波探測(cè)器通過(guò)天線等部分勻速轉(zhuǎn)動(dòng),形成天線波束的圓錐掃描,由此獲取地球表面和大氣的輻射數(shù)據(jù)。因此,在對(duì)目標(biāo)云團(tuán)進(jìn)行觀測(cè)時(shí),成像儀的天頂角也會(huì)造成云的位置偏差。本研究中也對(duì)此進(jìn)行了探討。成像儀的天頂角的定義成像儀和觀測(cè)點(diǎn)間的連線與成像儀和星下點(diǎn)間的連線間的夾角,即若星下點(diǎn)的位置為(lonA,latA),真云位置(觀測(cè)點(diǎn))的位置為(lonB,latB),則觀測(cè)點(diǎn)的天頂角為如在圖2中,當(dāng)角OSB 增大時(shí),B 點(diǎn)和C 點(diǎn)間的距離也會(huì)增大,當(dāng)微波儀器空間分辨率較高時(shí),B 點(diǎn)和C 點(diǎn)間可能間隔了數(shù)個(gè)像元,即產(chǎn)生了視差問(wèn)題。圖9(見(jiàn)下頁(yè))為基于FY-4A 的各個(gè)像元的經(jīng)緯度位置的模擬結(jié)果,模擬云高分別為1 km、5 km、10 km和15 km,空間分辨率設(shè)定為20 km、40 km、60 km和80 km。在這樣的模擬條件下,每個(gè)像元處的天頂角不同,產(chǎn)生的云位置偏差也不同,分別以天頂角為x軸,產(chǎn)生的云位置偏差為y軸,以散點(diǎn)圖的形式生成如圖結(jié)果。結(jié)合圖6~圖9,我們發(fā)現(xiàn)云頂高度越高,產(chǎn)生的視差問(wèn)題越大;在同等云高條件下,天頂角越大,產(chǎn)生的視差越大;且分辨率并沒(méi)有直接對(duì)視差產(chǎn)生影響。
圖9 天頂角、云頂高度及分辨率與對(duì)應(yīng)像素點(diǎn)處云位置偏差的聯(lián)系
為了能更為準(zhǔn)確直觀地展示靜止微波探測(cè)器上可能出現(xiàn)視場(chǎng)偏差,我們重點(diǎn)針對(duì)臺(tái)風(fēng)和高云個(gè)例進(jìn)行分析,包括2018年7月10日的臺(tái)風(fēng)“瑪莉亞”和同天在中國(guó)陸面的某個(gè)云頂高度大于15 km的對(duì)流云區(qū)。
圖10 展示的是2018 年第8 號(hào)超強(qiáng)臺(tái)風(fēng)“瑪莉亞”?!艾斃騺啞庇?018年7月4日(北京時(shí)間,下同)在美國(guó)關(guān)島以東的洋面生成,于2018 年7 月10 日11 時(shí)左右到達(dá)我國(guó)的臺(tái)風(fēng)預(yù)警24 小時(shí)警戒線,最后于2018 年7 月11 日09 時(shí)左右在福建連江黃岐半島登陸,登陸時(shí)中心最大風(fēng)速高達(dá)42 m/s,造成的直接經(jīng)濟(jì)損失接近30 億元,超百萬(wàn)人受災(zāi)。圖11 是將圖10 中116~123 °E,20~30 °N 的區(qū)域放大的結(jié)果。結(jié)合兩幅圖可以看出,隨著空間分辨率的降低,靜止衛(wèi)星的云掩膜產(chǎn)品會(huì)產(chǎn)生誤判,F(xiàn)Y-4A 的云產(chǎn)品中判別為無(wú)云的地區(qū)也會(huì)因?yàn)樵频奈恢闷疃鴱膱D上顯示為有云,F(xiàn)Y-4A 的云產(chǎn)品中判別為有云的地區(qū)也會(huì)出現(xiàn)誤判為無(wú)云。這個(gè)誤差也會(huì)影響衛(wèi)星對(duì)臺(tái)風(fēng)的定位、追蹤與強(qiáng)度估計(jì)。
相比起極軌微波探測(cè)器,靜止軌道微波探測(cè)器的時(shí)間分辨率更高,在臺(tái)風(fēng)和強(qiáng)對(duì)流的定位、追蹤等方面的應(yīng)用上有著明顯優(yōu)勢(shì),但視差問(wèn)題會(huì)對(duì)其應(yīng)用產(chǎn)生影響。圖10、圖11和圖12中也展示了同一天的4個(gè)時(shí)刻,不同空間分辨率下視差訂正前后的云的位置,盡管高時(shí)間分辨率的微波探測(cè)器能對(duì)臺(tái)風(fēng)進(jìn)行連續(xù)追蹤,但視差問(wèn)題依舊會(huì)影響臺(tái)風(fēng)定位的準(zhǔn)確性。
圖10 不同分辨率下臺(tái)風(fēng)“瑪莉亞”的位置與分辨率為4 km時(shí)觀測(cè)得到的臺(tái)風(fēng)位置的對(duì)比(視差影響)紅線為FY-4A觀測(cè)得到的云的輪廓。
圖11 不同分辨率下臺(tái)風(fēng)“瑪莉亞”的位置與分辨率為4 km時(shí)觀測(cè)得到的臺(tái)風(fēng)位置的對(duì)比
視差的影響在分析高云的例子時(shí)體現(xiàn)得更為明顯。圖12 中是2018 年7 月10 日中國(guó)陸面上的云頂高度大于15 km 的對(duì)流云個(gè)例其視差訂正前后差異對(duì)比。從圖12 中可以清楚地看到,與靠近赤道的臺(tái)風(fēng)相比,當(dāng)分辨率降低后,緯度稍高的高云例子在進(jìn)行視差訂正前后云的位置偏差變得更大,當(dāng)分辨率達(dá)40 km 及以上時(shí),甚至?xí)霈F(xiàn)觀測(cè)時(shí)缺漏了某些云團(tuán)的情況,這無(wú)疑會(huì)對(duì)陸面上對(duì)流系統(tǒng)的觀測(cè)造成極大的影響。
圖12 不同分辨率下高云的位置與分辨率為4 km時(shí)觀測(cè)得到的高云位置的對(duì)比
本次研究通過(guò)利用數(shù)值模擬和實(shí)際風(fēng)云四號(hào)靜止氣象衛(wèi)星A 星觀測(cè)結(jié)果對(duì)未來(lái)靜止軌道微波探測(cè)器的視差特征進(jìn)行了詳細(xì)分析。
(1) 靜止軌道微波探測(cè)器在探測(cè)過(guò)程中會(huì)受到云高的影響,觀測(cè)時(shí)會(huì)造成一定的視差,較大的視差會(huì)明顯影響預(yù)報(bào)員對(duì)天氣過(guò)程的實(shí)時(shí)監(jiān)測(cè)分析。隨著云高和觀測(cè)視角的增加,造成的觀測(cè)視差越大,必要時(shí)需對(duì)云的位置進(jìn)行有效訂正。
(2) 靜止微波探測(cè)器較粗的空間分辨率不會(huì)對(duì)視差的大小直接造成影響,但會(huì)對(duì)視差訂正前后的亮溫變化產(chǎn)生影響,空間分辨率越高,視差問(wèn)題影響越大。
(3) 盡管靜止微波探測(cè)器具有較高的時(shí)間分辨率,但是在臺(tái)風(fēng)或強(qiáng)對(duì)流的預(yù)警過(guò)程中,視差的存在會(huì)影響預(yù)警的準(zhǔn)確性。
除此之外,因微波具有云雨穿透能力,紅外探測(cè)所得的輻射云頂高度不一定等同于微波探測(cè)得到的真實(shí)云頂高度,當(dāng)對(duì)靜止微波探測(cè)進(jìn)行視差訂正時(shí)需加以考慮。