曹 倩 , 章毅之 , 易艷紅 , 劉雅楠 , 劉良玉
1.江西省氣象科學(xué)研究所, 江西 南昌 3300962.宜春市氣象局, 江西 宜春 336000
颮線(xiàn)是由一系列活躍的對(duì)流云帶側(cè)向排列而成的線(xiàn)狀或帶狀的中尺度對(duì)流系統(tǒng),過(guò)境時(shí)易產(chǎn)生雷暴、暴雨、大風(fēng)、冰雹和龍卷等天氣。其中,弓形颮線(xiàn)常常與下?lián)舯┝骰驈?qiáng)烈的地面直線(xiàn)大風(fēng)相伴(Nolen,1959;Wakimoto et al,2006;Wheatley et al,2006),比普通的直線(xiàn)形颮線(xiàn)更易產(chǎn)生強(qiáng)災(zāi)害性天氣(Przybylinski,1995;Atkins et al,2005;梁建宇和孫建華,2012;劉香娥和郭學(xué)良,2012),因此一直是國(guó)內(nèi)外同行研究的重點(diǎn)和難點(diǎn)。
近年來(lái),對(duì)弓形颮線(xiàn)的研究工作較多。例如,Weisman(2001)通過(guò)數(shù)值模擬研究發(fā)現(xiàn)后側(cè)入流急流是形成弓形回波的重要原因。金龍等(2013)利用多種觀(guān)測(cè)資料和雙多普勒雷達(dá)反演的三維風(fēng)場(chǎng),分析一次弓形回波的三維結(jié)構(gòu)和演變機(jī)制發(fā)現(xiàn),降水的拖曳及蒸發(fā)冷卻在地面形成強(qiáng)冷池,冷池觸發(fā)的干冷后側(cè)入流在弓形回波前側(cè)下沉,促使回波演變成弓形。公衍鐸等(2019)和羅琪等(2019)綜合多種觀(guān)測(cè)資料和NECP分析資料,分別分析了發(fā)生在我國(guó)北方的強(qiáng)弓狀颮線(xiàn)的形成和維持機(jī)制,均發(fā)現(xiàn)對(duì)流層中層強(qiáng)后側(cè)入流和大的溫度露點(diǎn)差導(dǎo)致強(qiáng)下層輻散氣流形成,是颮線(xiàn)演變成弓狀結(jié)構(gòu)的主要原因。
上述研究對(duì)弓形颮線(xiàn)形成的中尺度系統(tǒng)有重要的指示作用,而對(duì)影響其演變機(jī)制的揭示還不夠,特別是地形對(duì)颮線(xiàn)發(fā)展演變的影響涉及較少,使得颮線(xiàn)上山、下山過(guò)程中如何演變很難把握,進(jìn)而難以預(yù)測(cè)其引發(fā)災(zāi)害的強(qiáng)度及位置。江西省坐落于長(zhǎng)江中下游,中部丘陵起伏,東西南三面環(huán)山,在地形的影響下強(qiáng)對(duì)流天氣多發(fā)(鄧詩(shī)茹和鄒海波,2015;劉亞楠等,2021;張娟娟等,2021)。以往造成江西地區(qū)大范圍雷暴大風(fēng)的颮線(xiàn)大多產(chǎn)生于臨近省份,并逐漸移入江西。為了提高江西地區(qū)颮線(xiàn)天氣的預(yù)報(bào)預(yù)警能力,嘗試對(duì)移入江西的颮線(xiàn)做地形敏感性試驗(yàn),以探究地形對(duì)颮線(xiàn)演變特征的影響是非常有必要的。因此,文中首先利用GSI-3DVar同化系統(tǒng)同化多部多普勒雷達(dá)的徑向風(fēng)和反射率因子資料獲得一個(gè)相對(duì)準(zhǔn)確的初始場(chǎng)后,采用WRF模式對(duì)2016年4月16日對(duì)由湖南移入江西的一次弓形颮線(xiàn)天氣過(guò)程進(jìn)行數(shù)值模擬研究和地形敏感性試驗(yàn),以揭示地形對(duì)弓形颮線(xiàn)形成和演變的影響。
文中研究的強(qiáng)對(duì)流天氣于2016年4月16日01時(shí)(北京時(shí),下同)在湖南生成,此后逐步發(fā)展并快速向東移動(dòng),于04時(shí)演變成尺度較小的直線(xiàn)形颮線(xiàn)并進(jìn)入江西,之后直線(xiàn)形颮線(xiàn)經(jīng)過(guò)位于江西、湖南兩省交界處的九嶺山后演變成標(biāo)準(zhǔn)的弓形颮線(xiàn),08時(shí)進(jìn)入消散階段。此次颮線(xiàn)過(guò)程伴隨雷暴大風(fēng)和短時(shí)強(qiáng)降水天氣,造成27個(gè)縣(市)出現(xiàn)短時(shí)強(qiáng)降水(≥20 mm/h),以都昌縣46.2 mm/h為最大,以及8個(gè)區(qū)域自動(dòng)氣象站出現(xiàn)短時(shí)10級(jí)以上雷雨大風(fēng),其中以金溪縣合市站33.8 m/s為最大。
利用美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心(NCEP)提供的每6 h一次、分辨率為1°×1°的再分析資料,分析此次強(qiáng)對(duì)流天氣過(guò)程的大尺度環(huán)境背景。2016年4月15日20時(shí)500 hPa環(huán)流形勢(shì)(圖略)顯示,中高緯地區(qū),巴爾喀什湖一帶存在一個(gè)強(qiáng)大的高壓脊,脊前不斷有冷空氣沿新疆西部進(jìn)入我國(guó)。中低緯地區(qū),從河套地區(qū)到華北地區(qū)低槽加深東移,低槽后部不斷有冷空氣擴(kuò)散南下,青藏高原地區(qū)有一低槽位于四川省西部,江西省中北部處于槽前西南急流中,并且溫度槽落后于高度槽,強(qiáng)烈的斜壓效應(yīng)使得高空低槽不斷發(fā)展東移,到16日02時(shí)低槽已經(jīng)東移至重慶和貴州東部。槽后擴(kuò)散南下的冷空氣與西南暖濕氣流在湖南地區(qū)上空交匯是觸發(fā)此次強(qiáng)對(duì)流天氣的重要因素。分析850 hPa和700 hPa形勢(shì)場(chǎng)(圖1)發(fā)現(xiàn),湖南省西部和江西省中北部均受西南暖濕氣流控制,其中700 hPa層的相對(duì)濕度為60%—80%,而850 hPa層上大部分區(qū)域相對(duì)濕度均大于80%。西南暖濕氣流為颮線(xiàn)的發(fā)生提供了充足的水汽,且低層水汽比高層更加充沛。850 hPa層上湖南地區(qū)出現(xiàn)了西北風(fēng)與西南風(fēng)切變,這也有利于觸發(fā)此次颮線(xiàn)天氣。
九嶺山為東北—西南走向狹長(zhǎng)的條狀山峰,最大高度超過(guò)1 500 m(圖2)。圖3給出了直線(xiàn)形颮線(xiàn)經(jīng)過(guò)九嶺山逐漸演變成弓形颮線(xiàn)的過(guò)程實(shí)況。16日04時(shí)左右線(xiàn)狀強(qiáng)回波帶從湖南移入江西宜春市(圖3a),此時(shí)其水平尺度還不大,層云降水區(qū)也不寬。之后,線(xiàn)狀強(qiáng)回波帶逐漸東移爬山,其東北側(cè)有分散的對(duì)流單體生成并與主體合并(圖3b、c)。05時(shí),颮線(xiàn)結(jié)構(gòu)清晰,宜春市存在一條東北—西南向帶狀強(qiáng)對(duì)流回波區(qū),最大回波強(qiáng)度超過(guò)55 dBz,此時(shí)颮線(xiàn)前1 h所經(jīng)之處均變成層狀云降水區(qū)。此后,颮線(xiàn)主體東移下山,其形狀由直線(xiàn)形逐漸演變成弓形,中心回波強(qiáng)度均維持在55dBz以上(圖3d、e)。
圖1 2016年4月16日02時(shí)700 hPa(a)、850 hPa(b)風(fēng)場(chǎng)(矢線(xiàn),單位:m/s)、溫度場(chǎng)(紅實(shí)線(xiàn),單位:℃)和相對(duì)濕度場(chǎng)(填色,單位:%)
圖2 九嶺山海拔分布(單位:m;白色虛線(xiàn)內(nèi)為九嶺山)
06:30,颮線(xiàn)主體發(fā)展成熟,強(qiáng)對(duì)流區(qū)高度組織化,成為一個(gè)標(biāo)準(zhǔn)的弓形颮線(xiàn),并伴有寬廣的層云降水區(qū),以及在強(qiáng)對(duì)流區(qū)和層云區(qū)中間存在有界弱回波區(qū)(圖3f)。
文中,采用Liang(2007)和Chen等(2017)改進(jìn)后的基于IVAP( integrating velocity-azimuth process)方法的雷達(dá)徑向風(fēng)觀(guān)測(cè)算子,在同化雷達(dá)觀(guān)測(cè)時(shí)觀(guān)測(cè)參數(shù)由雷達(dá)徑向風(fēng)及其空間分布特征計(jì)算而來(lái),同時(shí)同化徑向風(fēng)和切向風(fēng)信息。IVAP觀(guān)測(cè)算子的觀(guān)測(cè)空間(Y1、Y2)和分析空間(H1、H2)分別為:
(1)
(2)
GSI(Gridpoint Statistical Interpolation)云分析模塊可以結(jié)合地面觀(guān)測(cè)、衛(wèi)星觀(guān)測(cè)、雷達(dá)觀(guān)測(cè)等多種觀(guān)測(cè)資料計(jì)算云量,反演計(jì)算云冰、云水、雪、雨、雹等混合比。文中采用GSI中的云分析方法同化雷達(dá)反射率因子。使用RUC層云方案(Weygandt et al,2006)計(jì)算云水和冰的混合比,使用Thomps-on方案(Thompson et al,2004)計(jì)算雪、雨和雹的混合比。采用APPS方案(Hu et al,2006)調(diào)整云內(nèi)的溫度(假定云內(nèi)溫度垂直變化為濕絕熱過(guò)程),使云中溫度與云場(chǎng)一致。采用RUC層云方案(Weygandt et al,2006)調(diào)整云中的水汽。
圖3 觀(guān)測(cè)的2016年4月16日04:00(a)、04:30(b)、05:00(c)、05:30(d)、06:00(e)、06:30(f)雷達(dá)組合反射率(單位:dBz)
研究使用WRF-ARW V3.9.1和GSI-3DVar V3.4分別作為預(yù)報(bào)模式和資料同化系統(tǒng)。采用三重嵌套方案,模擬區(qū)域中心為(30°N,114°E)(圖略)。外層網(wǎng)格數(shù)為301×301,水平格距為9 km;第二層網(wǎng)格數(shù)為391×397,水平格距為3 km;最內(nèi)層網(wǎng)格數(shù)為601×466,水平格距為1 km。垂直方向有50層,模式頂高為50 hPa。由每6 h一次、分辨率為1°×1°的NCEP再分析資料,通過(guò)WPS初始化模塊處理得到2016年4月16日02時(shí)的插值場(chǎng),以之作為初始背景場(chǎng)和側(cè)邊界場(chǎng)。同時(shí),采用GSI自帶的云分析方法(Hu et al,2006)同化江西省南昌、宜春、景德鎮(zhèn)、上饒和撫州5部多普勒雷達(dá)的反射率因子資料,又采用IVAP方法(Liang,2007;Chen et al,2017)同化上述5部雷達(dá)的徑向風(fēng)資料。之后,以02時(shí)同化后的分析場(chǎng)作為初始場(chǎng),向前預(yù)報(bào)6 h至08時(shí)結(jié)束,積分時(shí)間步長(zhǎng)為30 s,每30 min輸出一次結(jié)果。模式主要參數(shù)化方案為WDM6云微物理方案、RRTM長(zhǎng)波輻射方案、Dudhia短波輻射方案、Noah陸面過(guò)程、YSU邊界層方案、Monin-Obukhov地表方案,在第一重嵌套中采用Kain-Fritsch積云對(duì)流方案,第二重和第三重嵌套中關(guān)閉積云對(duì)流方案。文中主要對(duì)第三層網(wǎng)格的模擬結(jié)果進(jìn)行分析。
為了探討九嶺山地形對(duì)直線(xiàn)形颮線(xiàn)演變成弓形颮線(xiàn)過(guò)程的影響,設(shè)計(jì)了兩組試驗(yàn):1)控制試驗(yàn)(CTLT),模式使用真實(shí)的地形,模擬此次颮線(xiàn)經(jīng)過(guò)九嶺山后的發(fā)展演變過(guò)程;2) 模擬過(guò)程中將九嶺山地形高度減半(HALF),以分析地形對(duì)弓形颮線(xiàn)形成的影響。
將控制試驗(yàn)?zāi)M的雷達(dá)組合反射率(圖4)與實(shí)況觀(guān)測(cè)(圖3)對(duì)比分析發(fā)現(xiàn),16日04:00(圖4a),線(xiàn)狀強(qiáng)回波帶從宜春市西部進(jìn)入江西,其北側(cè)出現(xiàn)了與實(shí)況類(lèi)似的新生對(duì)流單體;隨著線(xiàn)狀強(qiáng)回波帶東移爬山,北側(cè)有分散的對(duì)流單體與其合并,強(qiáng)回波帶水平尺度逐漸變大(圖4b、c);此后,系統(tǒng)東移下山的過(guò)程中,強(qiáng)回波帶逐漸由直線(xiàn)形演變成弓形,中心強(qiáng)度大于55 dBz(圖4d—f)。對(duì)比整個(gè)模擬過(guò)程,控制試驗(yàn)?zāi)M的強(qiáng)對(duì)流區(qū)中心強(qiáng)度、位置和走向與實(shí)況有所差別,且模擬的層狀云降水區(qū)過(guò)??;但是控制試驗(yàn)?zāi)軌蚧驹佻F(xiàn)了直線(xiàn)形颮線(xiàn)經(jīng)過(guò)九嶺山后其尺度逐漸增大并演變成弓形颮線(xiàn)的過(guò)程。因此,CTLT試驗(yàn)的模擬結(jié)果可以用來(lái)探討地形對(duì)弓形颮線(xiàn)形成的影響。
HALF試驗(yàn)?zāi)M結(jié)果(圖5)顯示,將九嶺山地形高度減半后,模擬的颮線(xiàn)發(fā)展演變過(guò)程與CTLT試驗(yàn)?zāi)M結(jié)果有明顯區(qū)別??傮w來(lái)看,HALF試驗(yàn)和CTLT試驗(yàn)?zāi)M的颮線(xiàn)移動(dòng)速度均與實(shí)況一致。04:00—05:00時(shí)段(圖5a—c),HALF試驗(yàn)和CTLT試驗(yàn)?zāi)M的颮線(xiàn)略相似均表現(xiàn)為直線(xiàn)形颮線(xiàn)水平尺度逐漸增大并向東移動(dòng),但CTLT試驗(yàn)?zāi)M的強(qiáng)對(duì)流區(qū)較HALF試驗(yàn)略寬。颮線(xiàn)下山后(圖5d—f),HALF試驗(yàn)?zāi)M的颮線(xiàn)仍然表現(xiàn)為直線(xiàn)形并繼續(xù)向東移動(dòng),未能演變成弓形颮線(xiàn),模擬的颮線(xiàn)發(fā)展演變過(guò)程與實(shí)況有很大差距。由此可知,九嶺山地形直接影響颮線(xiàn)能否由直線(xiàn)形演變成弓形。
圖4 CTLT試驗(yàn)?zāi)M的2016年4月16日04:00(a)、04:30(b)、05:00(c)、05:30(d)、06:00(e)、06:30(f)雷達(dá)組合反射率(單位:dBz)
圖5 同圖4,但為HALF試驗(yàn)
后部入流急流和冷池是影響颮線(xiàn)演變成弓形的重要因素(Weisman,2001;金龍等,2013)。文中將通過(guò)對(duì)比分析兩組試驗(yàn)的三維風(fēng)場(chǎng)結(jié)構(gòu)和近地面冷池,探討地形對(duì)弓狀颮線(xiàn)形成影響的內(nèi)在原因。
圖6顯示了CTLT試驗(yàn)?zāi)M的颮線(xiàn)由直線(xiàn)形演變成弓形過(guò)程中不同高度的風(fēng)暴相對(duì)風(fēng)場(chǎng)(即水平風(fēng)場(chǎng)剔除颮線(xiàn)系統(tǒng)的水平移動(dòng)速度)。04:30(圖6a1—c1),直線(xiàn)形颮線(xiàn)南段經(jīng)過(guò)九嶺山下山后,在1 km高度上出現(xiàn)了明顯的后側(cè)入流,其平均速度約為16.5 m/s,2 km高度上還未出現(xiàn)后側(cè)入流,在3 km 高度上出現(xiàn)風(fēng)速輻合。05:00—05:30時(shí)段(圖6a2—c2、6a3—c3),直線(xiàn)形颮線(xiàn)東移至完全下山過(guò)程中,在1 km和2 km高度上均出現(xiàn)了后側(cè)入流,其平均速度約分別為23.2 m/s和17.8 m/s,1 km高度上的氣流強(qiáng)于2 km高度,3 km高度上仍然存在風(fēng)速的輻合。直線(xiàn)形颮線(xiàn)東移下山后(圖6a4—c4、6a5—c5),颮線(xiàn)由直線(xiàn)形逐漸演變成弓形,其后側(cè)1 km和2 km高度上的后側(cè)入流急流仍然非常強(qiáng),其平均速度分別維持在24.5 m/s和19 m/s左右。在3 km高度上弓形颮線(xiàn)的北端出現(xiàn)了氣旋性渦旋,颮線(xiàn)后側(cè)也出現(xiàn)了一致的西風(fēng),但后向入流的強(qiáng)度明顯弱于低層。中低層后側(cè)入流急流和氣旋渦度的出現(xiàn)均有利于颮線(xiàn)由直線(xiàn)形演變成弓形。
HALF試驗(yàn)?zāi)M結(jié)果(圖7)顯示,降低九嶺山地形高度后,所有時(shí)次1 km和2 km高度上的后側(cè)入流氣流比CTLT試驗(yàn)均顯著偏小,尤其是氣流下山后,1 km和2 km高度上的后側(cè)入流氣流速度比CTLT試驗(yàn)分別偏小6 m/s和10 m/s左右,因此HALF試驗(yàn)中2 km高度上的后側(cè)入流氣流非常弱。HALF試驗(yàn)中,在3 km高度上存在與1 km和2 km高度同位置、同強(qiáng)度的后側(cè)入流,颮線(xiàn)前側(cè)未出現(xiàn)明顯的入流,這些也與CTLT試驗(yàn)結(jié)果存在明顯差異。HALF試驗(yàn)與CTLT試驗(yàn)相似之處為05:30以后在3 km高度上颮線(xiàn)的北端也出現(xiàn)了氣旋性渦旋。由此可知,氣流在經(jīng)過(guò)九嶺山地形后下山增強(qiáng),使得颮線(xiàn)后側(cè)入流急流明顯增強(qiáng),進(jìn)而導(dǎo)致颮線(xiàn)由直線(xiàn)形演變成弓形;當(dāng)?shù)匦尾淮嬖诤?,氣流下山則不會(huì)明顯增強(qiáng),颮線(xiàn)繼續(xù)維持直線(xiàn)形并繼續(xù)向前傳播。
圖8和圖9分別給出了CTLT試驗(yàn)和HALF試驗(yàn)?zāi)M的近地面小時(shí)變溫和950 hPa地面風(fēng)場(chǎng)。CTLT試驗(yàn)結(jié)果顯示,在04:00—05:00時(shí)段(圖8a—c)颮線(xiàn)東移下山過(guò)程中,冷池強(qiáng)度很弱,其中心值僅為-2 ℃,但在此過(guò)程中冷池面積逐漸增大,且伴隨冷池出現(xiàn)的冷池出流也逐漸增強(qiáng)。此后,繼續(xù)東移過(guò)程中(圖8d—e),冷池強(qiáng)度逐漸增大,冷池出流的強(qiáng)度也逐漸增大。當(dāng)颮線(xiàn)演變成標(biāo)準(zhǔn)弓狀回波時(shí)(圖8f),冷池中心強(qiáng)度低于-4 ℃。HALF試驗(yàn)結(jié)果(圖9)顯示,在04:00—05:00時(shí)段冷池及伴隨的冷池出流的強(qiáng)度僅略弱于CTLT試驗(yàn),但是颮線(xiàn)下山后的05:30—06:30時(shí)段,冷池及伴隨的冷池出流的強(qiáng)度與CTRL試驗(yàn)差距更大。
對(duì)比分析實(shí)況(圖10)和CTLT、HALF試驗(yàn)?zāi)M的04時(shí)、05時(shí)、06時(shí)地面小時(shí)變溫場(chǎng)發(fā)現(xiàn),相較于HALF試驗(yàn),CTRL試驗(yàn)?zāi)M的冷池位置、范圍和強(qiáng)度與實(shí)況更接近。由此可知,由于九嶺山地形的存在,颮線(xiàn)在下山東移過(guò)程中冷池強(qiáng)度增強(qiáng),導(dǎo)致颮線(xiàn)系統(tǒng)內(nèi)的溫度梯度增大,使得颮線(xiàn)后側(cè)的入流也增強(qiáng),進(jìn)而有利于颮線(xiàn)向弓形演變;當(dāng)?shù)匦尾淮嬖跁r(shí),颮線(xiàn)東移過(guò)程中其后側(cè)的冷池強(qiáng)度減弱,伴隨的冷池出流減弱,颮線(xiàn)系統(tǒng)內(nèi)的溫度梯度同時(shí)也減小,颮線(xiàn)后側(cè)的入流也減弱,此時(shí)颮線(xiàn)在東移過(guò)程中未能由直線(xiàn)形演變成弓形。
利用GSI-3DVar同化系統(tǒng)中IVAP雷達(dá)徑向風(fēng)觀(guān)測(cè)算子同化多部多普勒雷達(dá)徑向風(fēng)資料,采用云分析方法同化多部多普勒雷達(dá)的反射率資料后,將同化結(jié)果作為初始場(chǎng),再利用WRF模式對(duì)2016年4月16日發(fā)生在江西省西北部的一次颮線(xiàn)天氣過(guò)程進(jìn)行數(shù)值模擬,探討地形對(duì)颮線(xiàn)形成的影響及其原因,得到以下主要結(jié)論:
1) 此次颮線(xiàn)在高空低槽東移、西南急流、低層風(fēng)切變和地面弱冷空氣共同影響的背景下形成,水汽條件充沛。真實(shí)地形情景下模擬的雷達(dá)回波演變特征與實(shí)況相似,颮線(xiàn)在爬山的過(guò)程中尺度逐漸增大,下山后其形狀由直線(xiàn)形逐漸演變成弓形。
圖6 CTLT試驗(yàn)?zāi)M的2016年4月16日1 km(左)、2 km(中)、3 km(右)高度風(fēng)暴相對(duì)風(fēng)場(chǎng)(矢線(xiàn),單位:m/s)和雷達(dá)反射率因子(填色,單位:dBz)(藍(lán)色等值線(xiàn)表示海拔大于500 m地形)(a.04:30,b.05:00,c.05:30,d.06:00,e.06:30)
圖7 同圖6,但為HALF試驗(yàn)
2) 颮線(xiàn)從九嶺山東側(cè)下山過(guò)程中氣流增強(qiáng),形成較強(qiáng)的后側(cè)入流急流和近地面冷池,其形狀也由直線(xiàn)形逐漸演變成弓形;當(dāng)降低九嶺山地形高度后,颮線(xiàn)東移過(guò)程中后側(cè)入流氣流明顯減弱,且冷池范圍和強(qiáng)度均減弱,其形狀未能由直線(xiàn)形演變成弓形。
圖8 CTLT模擬的2016年4月16日04:00(a)、04:30(b)、05:00(c)、05:30(d)、06:00(e)、06:30(f)950 hPa風(fēng)場(chǎng)(矢線(xiàn),單位:m/s)和小時(shí)變溫場(chǎng)(填色,單位:℃)
圖9 同圖8,但為HALF試驗(yàn)
綜上所述,九嶺山地形直接影響颮線(xiàn)能否由直線(xiàn)形演變成弓形。然而,以上結(jié)論僅僅針對(duì)一次颮線(xiàn)過(guò)程,得到的結(jié)論不一定具有普適性,還需要綜合多個(gè)颮線(xiàn)個(gè)例進(jìn)行分析。并且,對(duì)于機(jī)理的分析還僅僅是定性分析,在今后的研究中需要通過(guò)定量計(jì)算來(lái)深入探究地形影響的內(nèi)在機(jī)制。
圖10 觀(guān)測(cè)的2016年4月16日04:00(a)、05:00(b)、06:00(c)小時(shí)變溫場(chǎng)(單位:℃)