第6章物理過(guò)程參數(shù)化和資料同化問(wèn)題

1、陳耀登0陳耀登0陳耀登參數(shù)化、初邊界條件數(shù)值天氣預(yù)報(bào) 大氣科學(xué)學(xué)院大氣科學(xué)學(xué)院 陳耀登陳耀登 QQ：3545881陳耀登1陳耀登1陳耀登參數(shù)化、初邊界條件課程任務(wù)課程任務(wù)數(shù)值天氣預(yù)報(bào)的基礎(chǔ)知識(shí)、基本原理和基本方法數(shù)值天氣預(yù)報(bào)的基礎(chǔ)知識(shí)、基本原理和基本方法課程目標(biāo)課程目標(biāo)理解數(shù)值天氣預(yù)報(bào)的基本原理；理解數(shù)值天氣預(yù)報(bào)的基本原理；掌握數(shù)值預(yù)報(bào)基本操作步驟；掌握數(shù)值預(yù)報(bào)基本操作步驟；為從事數(shù)值天氣預(yù)報(bào)業(yè)務(wù)、數(shù)值預(yù)報(bào)產(chǎn)品釋用、為從事數(shù)值天氣預(yù)報(bào)業(yè)務(wù)、數(shù)值預(yù)報(bào)產(chǎn)品釋用、數(shù)值模擬研究等打下基礎(chǔ)數(shù)值模擬研究等打下基礎(chǔ)；數(shù)值天氣預(yù)報(bào)數(shù)值天氣預(yù)報(bào)理論和方法理論和方法陳耀登2陳耀登2陳耀登參數(shù)化、初邊界條件天氣預(yù)報(bào)

2、的三次躍進(jìn)（1）地面和高空天氣圖的繪制和應(yīng)用； 天氣預(yù)報(bào)從經(jīng)驗(yàn)走向定性（2）數(shù)值天氣預(yù)報(bào)的實(shí)施； 天氣預(yù)報(bào)從定性走向定量（3）大氣混沌性：集合預(yù)報(bào) 單一的確定論走向多值概率陳耀登3陳耀登3陳耀登參數(shù)化、初邊界條件陳耀登4陳耀登4陳耀登參數(shù)化、初邊界條件20112011年年7 7月月2525號(hào)夜晚，威海乳山市局部突降特大暴雨號(hào)夜晚，威海乳山市局部突降特大暴雨，民房倒塌民房倒塌16001600多間，死亡和失蹤人數(shù)逾百人多間，死亡和失蹤人數(shù)逾百人1.11.1中尺度天氣及重要性中尺度天氣及重要性陳耀登5陳耀登5陳耀登參數(shù)化、初邊界條件衛(wèi)星云圖衛(wèi)星云圖江蘇江蘇北京北京中尺度降水系統(tǒng)中尺度降水系統(tǒng)(200

3、km-2000km)(200km-2000km)雷達(dá)回波雷達(dá)回波20km20112011年年7 7月月2525日晚日晚0808時(shí)時(shí)中尺度對(duì)流系統(tǒng)中尺度對(duì)流系統(tǒng)（20km-200km20km-200km）中尺度對(duì)流單體中尺度對(duì)流單體(2km20km)(2km20km)乳山站乳山站3 3小時(shí)小時(shí)325mm325mm降水降水1.11.1中尺度天氣及重要性中尺度天氣及重要性陳耀登6陳耀登6陳耀登參數(shù)化、初邊界條件陳耀登基本問(wèn)題、歷史及展望三份重要文件三份重要文件1.中國(guó)氣象局氣象科技創(chuàng)新體系建設(shè)指導(dǎo)氣象科技創(chuàng)新體系建設(shè)指導(dǎo)意見(jiàn)（意見(jiàn)（2014-2020年）年）2014.10.302.中國(guó)氣象局國(guó)家氣象

4、科技創(chuàng)新工程（國(guó)家氣象科技創(chuàng)新工程（20142020年）實(shí)施方案年）實(shí)施方案2014.10.303.教育部、中國(guó)氣象局教育部、中國(guó)氣象局關(guān)于加強(qiáng)氣象人才關(guān)于加強(qiáng)氣象人才培養(yǎng)工作的指導(dǎo)意見(jiàn)培養(yǎng)工作的指導(dǎo)意見(jiàn)2015.2.3陳耀登基本問(wèn)題、歷史及展望 中國(guó)氣象局：2014.10.30 氣象科技創(chuàng)新體系建設(shè)指導(dǎo)意見(jiàn)（氣象科技創(chuàng)新體系建設(shè)指導(dǎo)意見(jiàn)（2014-2020年）年）國(guó)家氣象科技創(chuàng)新工程（國(guó)家氣象科技創(chuàng)新工程（20142020年）實(shí)施方案年）實(shí)施方案陳耀登基本問(wèn)題、歷史及展望 一、優(yōu)化氣象相關(guān)專業(yè)和人才結(jié)構(gòu)。一、優(yōu)化氣象相關(guān)專業(yè)和人才結(jié)構(gòu)。辦好大氣科學(xué)大氣科學(xué)、應(yīng)用氣象學(xué)、大氣物理與大氣環(huán)境等學(xué)

5、科專業(yè)，加強(qiáng)數(shù)值天氣預(yù)報(bào)數(shù)值天氣預(yù)報(bào)、大氣探測(cè)、公共氣象服務(wù)、氣象防災(zāi)減災(zāi)等方向的人才培養(yǎng)。 支持行業(yè)特色高校行業(yè)特色高校統(tǒng)籌招生計(jì)劃增量與存量，適度擴(kuò)大大氣科學(xué)類專業(yè)本科生和研究適度擴(kuò)大大氣科學(xué)類專業(yè)本科生和研究生招生規(guī)模生招生規(guī)模。教育部、中國(guó)氣象局教育部、中國(guó)氣象局關(guān)于加強(qiáng)氣象人關(guān)于加強(qiáng)氣象人才培養(yǎng)工作的指導(dǎo)意見(jiàn)才培養(yǎng)工作的指導(dǎo)意見(jiàn)2015.2.3陳耀登基本問(wèn)題、歷史及展望一個(gè)重要事件一個(gè)重要事件 2015.4.20：教育部、中國(guó)氣象局成立“中國(guó)氣象人才培養(yǎng)聯(lián)盟” 氣象領(lǐng)域高等院校（氣象類學(xué)院）、科研院所、企業(yè)以及國(guó)家級(jí)氣象業(yè)務(wù)單位、省（自治區(qū)、直轄市）氣象局均可申請(qǐng)加入聯(lián)盟。陳耀登基本

6、問(wèn)題、歷史及展望宇如聰：解讀（宇如聰：解讀（2015.4.23） 氣象教育歷久彌新 學(xué)科建設(shè)與人才培養(yǎng)面臨挑戰(zhàn)： 天氣預(yù)報(bào)、數(shù)值模擬、大氣化學(xué)、地球觀測(cè)、人工影響天氣等方面仍有不足 陳耀登12陳耀登12陳耀登參數(shù)化、初邊界條件主要內(nèi)容主要內(nèi)容第一章第一章 大氣運(yùn)動(dòng)的基本方程組大氣運(yùn)動(dòng)的基本方程組 模式的基本方程組(球坐標(biāo)、局地直角坐標(biāo))；方程的垂直坐標(biāo)變換； 模式的分類。 第二章第二章 地圖投影坐標(biāo)系中的基本方程組地圖投影坐標(biāo)系中的基本方程組 地圖投影的基本概念； 三種常用的地圖投影方式的基本概念、特點(diǎn)；正交曲線坐標(biāo)系的基本知識(shí)；地圖投影坐標(biāo)系下的方程：水平坐標(biāo)變換。第四章 準(zhǔn)地轉(zhuǎn)正壓模式第五

7、章 正壓原始方程模式第六章 斜壓原始方程模式第七章 物理過(guò)程參數(shù)化和初邊界條件第八章 數(shù)值天氣預(yù)報(bào)基本問(wèn)題及展望陳耀登13陳耀登13陳耀登參數(shù)化、初邊界條件主要內(nèi)容主要內(nèi)容第三章 數(shù)值計(jì)算方案差分格式(前差格式、后差格式、中央差格式)；差分格式的基本性質(zhì)（相容性、穩(wěn)定性，收斂性）；時(shí)間積分格式（時(shí)間層，隱式，非隱式）；差分格式的誤差問(wèn)題（縮小格距或采用較高階的差分格式均能提高計(jì)算的精確度）非線性計(jì)算穩(wěn)定性問(wèn)題（混淆誤差）陳耀登14陳耀登14陳耀登參數(shù)化、初邊界條件主要內(nèi)容主要內(nèi)容第四章第四章 正壓原始方程模式正壓原始方程模式正壓原始方程(基本假定)的導(dǎo)出及其積分性質(zhì)；正壓模式計(jì)算穩(wěn)定性及其判據(jù)

8、（空間步長(zhǎng)和時(shí)間步長(zhǎng)的關(guān)系）；正壓模式的空間差分格式、時(shí)間積分格式；空間濾波和平滑（抑制非線性不穩(wěn)定）；求解過(guò)程（上機(jī)實(shí)習(xí)）。陳耀登15陳耀登15陳耀登參數(shù)化、初邊界條件主要內(nèi)容主要內(nèi)容第五章第五章 斜壓原始方程模式斜壓原始方程模式斜壓引入的原因；斜壓原始方程的垂直坐標(biāo)。第六章 物理過(guò)程參數(shù)化和資料同化第七章 數(shù)值天氣預(yù)報(bào)基本問(wèn)題及展望陳耀登16陳耀登16陳耀登參數(shù)化、初邊界條件個(gè)人簡(jiǎn)介（陳耀登）個(gè)人簡(jiǎn)介（陳耀登） 福建沙縣人，博士，副教授。 大氣科學(xué)學(xué)院數(shù)值模擬系主任。 主要從事氣象數(shù)值模擬與資料同化相關(guān)工作。 2007年畢業(yè)于河海大學(xué)，獲博士學(xué)位；同年到我校大氣科學(xué)學(xué)院參加工作。陳耀登17

9、陳耀登17陳耀登參數(shù)化、初邊界條件模式物理過(guò)程參數(shù)化模式物理過(guò)程參數(shù)化 陳耀登18陳耀登18陳耀登參數(shù)化、初邊界條件1 1分鐘分鐘1 1小時(shí)小時(shí)熱對(duì)流龍卷等雷暴對(duì)流單體城市效應(yīng)中中尺度尺度梅雨鋒臺(tái)風(fēng)等斜壓波1 1天天1 1月月20m200m2km20km200km2000km10000km龍卷龍卷熱對(duì)流熱對(duì)流微尺度中氣旋中氣旋大龍卷大龍卷短重力短重力波波微尺度雷暴雷暴對(duì)流單對(duì)流單體體城市效城市效應(yīng)應(yīng)中尺度中尺度中尺度對(duì)流系對(duì)流系統(tǒng)統(tǒng)山地湖山地湖泊天氣泊天氣中尺度組織化雷暴雷暴對(duì)流對(duì)流單體單體對(duì)流對(duì)流系統(tǒng)系統(tǒng)大大尺度尺度 中中尺度尺度2 20 km30min30min 6h6h時(shí)時(shí)間間空間空間小

10、尺度小尺度雷暴雷暴3h3h 1day1day中中尺度尺度20 200 km大氣運(yùn)動(dòng)是多尺度的大氣運(yùn)動(dòng)是多尺度的陳耀登19陳耀登19陳耀登參數(shù)化、初邊界條件陳耀登20陳耀登20陳耀登參數(shù)化、初邊界條件數(shù)值天氣預(yù)報(bào)中的主要物理過(guò)程數(shù)值天氣預(yù)報(bào)中的主要物理過(guò)程1、濕熱力過(guò)程（凝結(jié)過(guò)程）、濕熱力過(guò)程（凝結(jié)過(guò)程） 1）積云對(duì)流過(guò)程）積云對(duì)流過(guò)程 2）大尺度凝結(jié)過(guò)程）大尺度凝結(jié)過(guò)程陳耀登21陳耀登21陳耀登參數(shù)化、初邊界條件數(shù)值天氣預(yù)報(bào)中的主要物理過(guò)程數(shù)值天氣預(yù)報(bào)中的主要物理過(guò)程2、邊界層過(guò)程（近地面和行星邊界）、邊界層過(guò)程（近地面和行星邊界）陳耀登22陳耀登22陳耀登參數(shù)化、初邊界條件數(shù)值天氣預(yù)報(bào)中的主

11、要物理過(guò)程數(shù)值天氣預(yù)報(bào)中的主要物理過(guò)程3、輻射過(guò)程、輻射過(guò)程4、微物理過(guò)程、微物理過(guò)程陳耀登23陳耀登23陳耀登參數(shù)化、初邊界條件WRF Physics Radiation 輻射方案 Longwave (ra_lw_physics) Shortwave (ra_sw_physics) Surface Surface layer (sf_sfclay_physics) 近地面層方案 Land/water surface (sf_surface_physics) 陸面方案 PBL (bl_pbl_physics) 行星邊界層方案 Cumulus (cu_physics) 積云對(duì)流方案 Microp

12、hysics (mp_physics) 微物理方案陳耀登24陳耀登24陳耀登參數(shù)化、初邊界條件 物理過(guò)程參數(shù)化的必要性物理過(guò)程參數(shù)化的必要性 次網(wǎng)格過(guò)程：次網(wǎng)格過(guò)程： 由于大氣是一種連續(xù)運(yùn)動(dòng)的由于大氣是一種連續(xù)運(yùn)動(dòng)的連續(xù)介質(zhì)連續(xù)介質(zhì)，不管模式的分辨率如何高，總有一些，不管模式的分辨率如何高，總有一些接接近于或小于網(wǎng)格距尺度的運(yùn)動(dòng)或過(guò)程近于或小于網(wǎng)格距尺度的運(yùn)動(dòng)或過(guò)程，無(wú)，無(wú)法在模式中確切地反映出來(lái)，這種運(yùn)動(dòng)過(guò)法在模式中確切地反映出來(lái)，這種運(yùn)動(dòng)過(guò)程稱為程稱為次網(wǎng)格過(guò)程次網(wǎng)格過(guò)程。陳耀登25陳耀登25陳耀登參數(shù)化、初邊界條件陳耀登26陳耀登26陳耀登參數(shù)化、初邊界條件 次網(wǎng)格尺度物理過(guò)程：次網(wǎng)格尺

13、度物理過(guò)程： 湍流、對(duì)流、凝結(jié)和輻射過(guò)程湍流、對(duì)流、凝結(jié)和輻射過(guò)程都包含都包含有次網(wǎng)格過(guò)程。在數(shù)值預(yù)報(bào)中通過(guò)采用有次網(wǎng)格過(guò)程。在數(shù)值預(yù)報(bào)中通過(guò)采用參參數(shù)化的方法數(shù)化的方法來(lái)考慮這些過(guò)程，即用大尺度來(lái)考慮這些過(guò)程，即用大尺度變量來(lái)描述次網(wǎng)格過(guò)程對(duì)大尺度運(yùn)動(dòng)的變量來(lái)描述次網(wǎng)格過(guò)程對(duì)大尺度運(yùn)動(dòng)的統(tǒng)統(tǒng)計(jì)效應(yīng)計(jì)效應(yīng)。陳耀登27陳耀登27陳耀登參數(shù)化、初邊界條件（也是潛熱釋放過(guò)程）分為：（也是潛熱釋放過(guò)程）分為：大尺度上升運(yùn)動(dòng)大尺度上升運(yùn)動(dòng)造成造成:上升慢，范圍廣，鋒面降水上升慢，范圍廣，鋒面降水積云對(duì)流積云對(duì)流產(chǎn)生：上升快，范圍小，熱帶降水產(chǎn)生：上升快，范圍小，熱帶降水陳耀登28陳耀登28陳耀登參數(shù)化、

14、初邊界條件熱力學(xué)方程 （差分）絕熱運(yùn)動(dòng)及溫度的擴(kuò)散造成的溫度變化國(guó)家氣象中心原先使用的LFM模式中的方案水汽凝結(jié)產(chǎn)生的非絕熱加熱水汽凝結(jié)產(chǎn)生的非絕熱加熱L為凝結(jié)潛熱為凝結(jié)潛熱陳耀登29陳耀登29陳耀登參數(shù)化、初邊界條件大尺度水平運(yùn)動(dòng)和垂直運(yùn)動(dòng)所引起的濕度變化水汽方程水汽收支變化水汽收支變化陳耀登30陳耀登30陳耀登參數(shù)化、初邊界條件水汽收支變化水汽收支變化水汽凝結(jié)產(chǎn)生的非絕熱加熱水汽凝結(jié)產(chǎn)生的非絕熱加熱大尺度凝結(jié)；大尺度凝結(jié)； 對(duì)流凝結(jié)；對(duì)流凝結(jié)；對(duì)流垂直輸送引起的濕度變化對(duì)流垂直輸送引起的濕度變化tqCLCCCFL為凝結(jié)潛熱為凝結(jié)潛熱陳耀登31陳耀登31陳耀登參數(shù)化、初邊界條件對(duì)流凝結(jié)對(duì)流凝

15、結(jié) 小尺度的積云對(duì)流小尺度的積云對(duì)流 大尺度運(yùn)動(dòng)總體影響的大尺度運(yùn)動(dòng)總體影響的關(guān)系及相互作用關(guān)系及相互作用 1 1）直接法：）直接法：求解描述兩類不同尺度運(yùn)動(dòng)的耦合方求解描述兩類不同尺度運(yùn)動(dòng)的耦合方程組，直接求解積云尺度運(yùn)動(dòng)程組，直接求解積云尺度運(yùn)動(dòng) 2 2）間接法：）間接法：用參數(shù)方法考慮小尺度運(yùn)動(dòng)對(duì)大尺度用參數(shù)方法考慮小尺度運(yùn)動(dòng)對(duì)大尺度運(yùn)動(dòng)的總體影響，而不考慮小尺度運(yùn)動(dòng)的細(xì)微結(jié)構(gòu)運(yùn)動(dòng)的總體影響，而不考慮小尺度運(yùn)動(dòng)的細(xì)微結(jié)構(gòu)陳耀登32陳耀登32陳耀登參數(shù)化、初邊界條件積云對(duì)流參數(shù)化將模式不能顯式分辨的對(duì)流凝結(jié)和對(duì)流引起的熱量、水分和動(dòng)量的輸送與模式的預(yù)報(bào)變量聯(lián)系起來(lái)。實(shí)際做法：用網(wǎng)格點(diǎn)上的變

16、量值來(lái)表示次網(wǎng)格尺度積云的凝結(jié)加熱以及垂直輸送效應(yīng)的作用。Cumulus Parameterization陳耀登33陳耀登33陳耀登參數(shù)化、初邊界條件（1）對(duì)流什么時(shí)候發(fā)生？（2）對(duì)流發(fā)生后產(chǎn)生什么影響？（3）參數(shù)的估計(jì)（1）對(duì)流產(chǎn)生的判別條件：a.條件不穩(wěn)定：濕靜力能模式條件l 積云對(duì)流積云對(duì)流陳耀登34陳耀登34陳耀登參數(shù)化、初邊界條件b.水汽輻合：c.低層風(fēng)速輻合： 滿足以上三個(gè)條件的格點(diǎn)可根據(jù)郭氏方法計(jì)算降水量。如果不考慮降落過(guò)程液態(tài)水的蒸發(fā)，氣柱中凝結(jié)的水汽就是地面降水量。（2）對(duì)流影響的估計(jì)：空氣拄的大尺空氣拄的大尺度輻合度輻合經(jīng)驗(yàn)常數(shù)經(jīng)驗(yàn)常數(shù)產(chǎn)生對(duì)流陳耀登35陳耀登35陳耀登參數(shù)

17、化、初邊界條件假如水汽輻合量中有一部分凝結(jié)下降，則降水量為：根據(jù)霍維邁爾經(jīng)驗(yàn)公式未凝結(jié)的水汽未凝結(jié)的水汽環(huán)境空氣環(huán)境空氣的比濕的比濕環(huán)境空氣環(huán)境空氣的溫度的溫度環(huán)境空氣的飽和比濕（2）對(duì)流影響的估計(jì)：陳耀登36陳耀登36陳耀登參數(shù)化、初邊界條件(a)對(duì)流凝結(jié)加熱率：（3）參數(shù)的估計(jì)：陳耀登37陳耀登37陳耀登參數(shù)化、初邊界條件(b)對(duì)流引起水汽垂直輸送，使空氣柱內(nèi)水汽重新分布，在模式中是直接計(jì)算大尺度水汽輻合和對(duì)流作用的總效應(yīng)。令對(duì)水汽方程則水汽方程可寫為則水汽方程可寫為（3）參數(shù)的估計(jì)：陳耀登38陳耀登38陳耀登參數(shù)化、初邊界條件最后，水汽方程可寫成：直接估計(jì)其中：Anthes（1977）未

18、凝結(jié)的水汽與各層的飽和差成正比陳耀登39陳耀登39陳耀登參數(shù)化、初邊界條件Cumulus schemeRecommendations about use For dx 10 km: probably need cumulus scheme For dx 3 km: probably do not need scheme However, there are cases where the earlier triggering of convection by cumulus schemes help For dx=3-10 km, scale separation is a question

19、No schemes are specifically designed with this range of scales in mind陳耀登40陳耀登40陳耀登參數(shù)化、初邊界條件三）、邊界層過(guò)程（動(dòng)量、熱量、水汽）三）、邊界層過(guò)程（動(dòng)量、熱量、水汽）l邊界層的處理：模式把最低層作為邊界層，采用經(jīng)驗(yàn)輸送公式計(jì)算湍流動(dòng)量輸送：湍流粘性應(yīng)力：邊界層頂：湍流應(yīng)力邊界層頂：湍流應(yīng)力0地面地面陳耀登41陳耀登41陳耀登參數(shù)化、初邊界條件四）、輻射過(guò)程四）、輻射過(guò)程輻射傳輸?shù)幕疚锢磉^(guò)程： 太陽(yáng)的輻射； 地-氣系統(tǒng)輻射； 不同地表狀態(tài)云、氣溶膠、水汽、臭氧、二氧化碳等對(duì)輻射傳輸?shù)挠绊?。陳耀?2陳耀登

20、42陳耀登參數(shù)化、初邊界條件四）、輻射過(guò)程四）、輻射過(guò)程輻射傳輸方程：在地球大氣條件下，求解非常復(fù)雜，只能在一些假定下求得解析解，因此輻射傳輸方程的求解，一直是大氣輻射學(xué)研究的重要內(nèi)容。一般通過(guò)參數(shù)化形式進(jìn)行表達(dá)并求解。陳耀登43陳耀登43陳耀登參數(shù)化、初邊界條件三）次網(wǎng)格過(guò)程參數(shù)化存在問(wèn)題：三）次網(wǎng)格過(guò)程參數(shù)化存在問(wèn)題： 盡管用參數(shù)化方案來(lái)描述次網(wǎng)格物理盡管用參數(shù)化方案來(lái)描述次網(wǎng)格物理過(guò)程的方法已取得了相當(dāng)好的效果，但仍過(guò)程的方法已取得了相當(dāng)好的效果，但仍有許多未解決的有許多未解決的問(wèn)題問(wèn)題:1）參數(shù)化不能考慮大尺度對(duì)小尺度的影響及參數(shù)化不能考慮大尺度對(duì)小尺度的影響及其反饋?zhàn)饔?；其反饋?zhàn)饔茫?/p>

21、2）參數(shù)的數(shù)值缺乏客觀的確定方法；參數(shù)的數(shù)值缺乏客觀的確定方法；3）模式對(duì)參數(shù)化的差異過(guò)于敏感等。模式對(duì)參數(shù)化的差異過(guò)于敏感等。 陳耀登44陳耀登44陳耀登參數(shù)化、初邊界條件中國(guó)科學(xué)：地球科學(xué)，2012，42（12）：1966-1978陳耀登45陳耀登45陳耀登參數(shù)化、初邊界條件陳耀登46陳耀登46陳耀登參數(shù)化、初邊界條件桑美黑格比BMJKF不同的cu scheme模擬的路徑與CMA的最佳路徑之間的路徑誤差隨時(shí)間增大，而且不同的cu scheme模擬的路徑差異也較大。陳耀登47陳耀登47陳耀登參數(shù)化、初邊界條件不同的cu scheme模擬的臺(tái)風(fēng)強(qiáng)度相差較大，在24h左右開(kāi)始出現(xiàn)較大差異，且差異

22、隨時(shí)間增大。陳耀登48陳耀登48陳耀登參數(shù)化、初邊界條件臺(tái)風(fēng)路徑和強(qiáng)度的模擬對(duì)積云對(duì)流參數(shù)化方案的選擇很敏感；不同個(gè)例的最佳方案不同；原因是：不同方案采用了不同的假定，以及對(duì)積云降水的處理。陳耀登49陳耀登49陳耀登參數(shù)化、初邊界條件中國(guó)科學(xué):地球科學(xué), 2014, 44(3): 548-562.陳耀登50陳耀登50陳耀登參數(shù)化、初邊界條件 邊界層內(nèi)地表和自由大氣的交換傳輸作用是十分復(fù)雜的物理過(guò)程 具體問(wèn)題需具體分析 個(gè)例分析的結(jié)果不具普適性 適當(dāng)選擇邊界層參數(shù)化方案對(duì)不同研究對(duì)象的模擬有重要意義（需要做敏感性實(shí)驗(yàn)分析其原因）陳耀登51陳耀登51陳耀登參數(shù)化、初邊界條件高原山地氣象研究，201

23、1，31（2）：18-25陳耀登52陳耀登52陳耀登參數(shù)化、初邊界條件熱帶氣象學(xué)報(bào)，2009，25（4）：435-441陳耀登53陳耀登53陳耀登參數(shù)化、初邊界條件熱帶氣象學(xué)報(bào)，2012，28（4）：461-470陳耀登54陳耀登54陳耀登參數(shù)化、初邊界條件一、數(shù)值預(yù)報(bào)中的初值問(wèn)題：一、數(shù)值預(yù)報(bào)中的初值問(wèn)題：初、邊界條件初、邊界條件 陳耀登55陳耀登55陳耀登參數(shù)化、初邊界條件一、數(shù)值預(yù)報(bào)中的初始化問(wèn)題：一、數(shù)值預(yù)報(bào)中的初始化問(wèn)題：采用原始方程模式作預(yù)報(bào)時(shí)，直接采用未經(jīng)任何處理的觀測(cè)值或分析值作為模式的初始場(chǎng)，容易導(dǎo)致高頻振蕩，產(chǎn)生計(jì)算不穩(wěn)定：第一、觀測(cè)或分析資料的誤差導(dǎo)致風(fēng)場(chǎng)和氣壓場(chǎng)的不平衡

24、；第一、觀測(cè)或分析資料的誤差導(dǎo)致風(fēng)場(chǎng)和氣壓場(chǎng)的不平衡；第二、初始資料和數(shù)值模式之間的不平衡。第二、初始資料和數(shù)值模式之間的不平衡。因此，需要對(duì)模式的初值進(jìn)行處理，即模式資料的模式資料的“初始初始化化”。問(wèn)題的提出：?jiǎn)栴}的提出：陳耀登56陳耀登56陳耀登參數(shù)化、初邊界條件模式初始化的常用方法：模式初始化的常用方法：靜力初始化動(dòng)力初始化變分初始化靜力初始化又稱靜處理，是指用一些已知的風(fēng)壓場(chǎng)平衡關(guān)系，或用運(yùn)動(dòng)方程等求得的診斷方程來(lái)處理初值，使風(fēng)場(chǎng)同氣壓場(chǎng)平衡或近似平衡的方法。地轉(zhuǎn)風(fēng)初值地轉(zhuǎn)風(fēng)初值平衡風(fēng)初值平衡風(fēng)初值考慮輻散的初始風(fēng)場(chǎng)考慮輻散的初始風(fēng)場(chǎng)陳耀登57陳耀登57陳耀登參數(shù)化、初邊界條件動(dòng)力初

25、始化又稱動(dòng)處理，這種方法是借助于原始方程模式本身所具有的動(dòng)力特性，經(jīng)過(guò)一些合理的步驟，使重力慣性波阻尼或被濾去，而得到接近平衡的初值。動(dòng)處理的具體方法包括動(dòng)力恢復(fù)法和正規(guī)波法等。變分初始化它既可以用于靜處理，也可以用于動(dòng)處理。 該方法是通過(guò)變分原理，使初始資料在一定動(dòng)力約束下調(diào)整，達(dá)到各種初始場(chǎng)之間協(xié)調(diào)一致的方法。陳耀登58陳耀登58陳耀登參數(shù)化、初邊界條件靜力初始化：簡(jiǎn)單、計(jì)算量小，但效果一般；動(dòng)力初始化：自動(dòng)調(diào)整、有效濾除重力波高頻部分、 計(jì)算量較大，非線性的處理不一定收斂變分初始化：物理意義清楚，結(jié)果好壞依賴于動(dòng)力 約束優(yōu)劣。簡(jiǎn)單約束：求解方便；復(fù)雜 約束：數(shù)學(xué)解存在較多問(wèn)題。常用方法評(píng)

26、價(jià)：常用方法評(píng)價(jià)：陳耀登59陳耀登59陳耀登參數(shù)化、初邊界條件一般初始化方法存在的問(wèn)題一般初始化方法存在的問(wèn)題消除虛假的重力波對(duì)計(jì)算穩(wěn)定性有好處，但可能也把重力波真實(shí)部分消除了數(shù)值預(yù)報(bào)的系統(tǒng)強(qiáng)度往往偏弱。上述方法都是讓風(fēng)壓場(chǎng)達(dá)到平衡或近似平衡，以消除虛假重力波，從而消除其帶來(lái)的計(jì)算不穩(wěn)定。天氣的發(fā)展演變和環(huán)流形式的劇烈變化，都是和風(fēng)壓場(chǎng)的不平衡激發(fā)出的重力慣性波有關(guān)，由于重力慣性波的存在，過(guò)分集中的不平衡能量能得到迅速釋放，使風(fēng)壓場(chǎng)重新調(diào)整適應(yīng)。如何使用不平衡初值以提高數(shù)值預(yù)報(bào)的準(zhǔn)確率，是研究和業(yè)如何使用不平衡初值以提高數(shù)值預(yù)報(bào)的準(zhǔn)確率，是研究和業(yè)務(wù)中重點(diǎn)關(guān)注的問(wèn)題之一。務(wù)中重點(diǎn)關(guān)注的問(wèn)題之一

27、。陳耀登60陳耀登60陳耀登參數(shù)化、初邊界條件二、數(shù)值預(yù)報(bào)中的邊界條件問(wèn)題：二、數(shù)值預(yù)報(bào)中的邊界條件問(wèn)題：數(shù)值模式需要給定邊界條件：水平側(cè)邊界條件垂直邊界條件2、法向速度為零的邊界條件（剛體邊界）3、海綿邊界條件（過(guò)渡帶：外側(cè)固定，里側(cè)為計(jì)算值）4、外推邊界條件常用的水平側(cè)邊界條件包括：1、固定的邊界條件(預(yù)報(bào)量不隨時(shí)間變化)大氣科學(xué)學(xué)院閔錦忠大氣科學(xué)學(xué)院閔錦忠陳耀登61陳耀登61陳耀登參數(shù)化、初邊界條件三、嵌套網(wǎng)格邊界處理問(wèn)題三、嵌套網(wǎng)格邊界處理問(wèn)題單向嵌套：?jiǎn)蜗蚯短祝?先計(jì)算粗網(wǎng)格，粗網(wǎng)格為細(xì)網(wǎng)格提供邊界； 然后再計(jì)算細(xì)網(wǎng)格。雙向嵌套：雙向嵌套： 先計(jì)算粗網(wǎng)格，粗網(wǎng)格為細(xì)網(wǎng)格提供邊界； 計(jì)

28、算細(xì)網(wǎng)格；在邊界交界處，將細(xì)網(wǎng)格結(jié)果反饋給粗網(wǎng)格。陳耀登62陳耀登62陳耀登參數(shù)化、初邊界條件二、數(shù)值預(yù)報(bào)中的資料同化問(wèn)題：二、數(shù)值預(yù)報(bào)中的資料同化問(wèn)題：u英國(guó)數(shù)學(xué)家英國(guó)數(shù)學(xué)家Richardson(1922)-數(shù)值預(yù)報(bào)嘗試數(shù)值預(yù)報(bào)嘗試陳耀登63陳耀登63陳耀登參數(shù)化、初邊界條件二、數(shù)值預(yù)報(bào)中的資料同化問(wèn)題：二、數(shù)值預(yù)報(bào)中的資料同化問(wèn)題：The major reasons for the failure of Richardsons forecast: 觀測(cè)站網(wǎng)的密度和資料的精確度不夠； 用的是完全的原始方程組； 所取的時(shí)間和空間的間隔不合適，使計(jì)算出現(xiàn)不穩(wěn)定陳耀登64陳耀登64陳耀登參數(shù)化、初

29、邊界條件數(shù)值預(yù)報(bào)的不確定性數(shù)值預(yù)報(bào)的不確定性 數(shù)值預(yù)報(bào)的不確定性大氣本身的混沌特性初值誤差模式誤差陳耀登65陳耀登65陳耀登參數(shù)化、初邊界條件二、數(shù)值預(yù)報(bào)中的資料同化問(wèn)題：二、數(shù)值預(yù)報(bào)中的資料同化問(wèn)題：如何最有效地利用各種氣象觀測(cè)資料，為數(shù)值天氣預(yù)報(bào)提供一個(gè)更好的初始場(chǎng)，是提高數(shù)值天氣預(yù)報(bào)準(zhǔn)確性的重要方法。這就是氣象資料同化問(wèn)題。（初始資料的同化）（初始資料的同化） 在進(jìn)行數(shù)值天氣預(yù)報(bào)的過(guò)程中，我們還能不斷地接收到氣象觀測(cè)資料。它們相對(duì)初始場(chǎng)而言，代表了某一地點(diǎn)、某一時(shí)段內(nèi)氣象場(chǎng)的實(shí)際演變。使用它們來(lái)不斷地修正原來(lái)的預(yù)報(bào)結(jié)果，對(duì)于提高數(shù)值天氣預(yù)報(bào)的質(zhì)量非常有益。（預(yù)報(bào)過(guò)程中的資料同化）（預(yù)報(bào)過(guò)

30、程中的資料同化）兩類資料同化的方法兩類資料同化的方法（初始資料的同化）（初始資料的同化）：預(yù)報(bào)過(guò)程中的資料同化預(yù)報(bào)過(guò)程中的資料同化：陳耀登67陳耀登67陳耀登參數(shù)化、初邊界條件陳耀登68陳耀登68陳耀登參數(shù)化、初邊界條件陳耀登69陳耀登69陳耀登參數(shù)化、初邊界條件陳耀登70陳耀登70陳耀登參數(shù)化、初邊界條件陳耀登71陳耀登71陳耀登參數(shù)化、初邊界條件Remote Sensing Observations陳耀登72陳耀登72陳耀登參數(shù)化、初邊界條件我國(guó)高空探測(cè)網(wǎng)我國(guó)高空探測(cè)網(wǎng) 全國(guó)120個(gè)探空站，每天進(jìn)行08、20時(shí)（北京時(shí)）兩次從地面到高空25000米以上高空探測(cè)。獲取高空大氣溫度、濕度、氣壓

31、和風(fēng)向風(fēng)速資料。 其中87個(gè)站參加全球資料交換。 共有7個(gè)探空站為全球GCOS探空站。陳耀登73陳耀登73陳耀登參數(shù)化、初邊界條件我國(guó)新一代天氣雷達(dá)布局我國(guó)新一代天氣雷達(dá)布局計(jì)劃布設(shè)計(jì)劃布設(shè)158158部部新一代天氣雷達(dá)新一代天氣雷達(dá); ;已經(jīng)完成了已經(jīng)完成了130130部建設(shè)部建設(shè), ,并投入并投入業(yè)務(wù)使用。業(yè)務(wù)使用。（不包括臺(tái)灣省、（不包括臺(tái)灣省、香港和澳門特區(qū)）香港和澳門特區(qū)）陳耀登74陳耀登74陳耀登參數(shù)化、初邊界條件四維同化是把不同時(shí)刻、不同地區(qū)、不同性質(zhì)的氣象資料，四維同化是把不同時(shí)刻、不同地區(qū)、不同性質(zhì)的氣象資料，通過(guò)統(tǒng)計(jì)與動(dòng)力關(guān)系使之在動(dòng)力與熱力上協(xié)調(diào)起來(lái)，以求得通過(guò)統(tǒng)計(jì)與動(dòng)力

32、關(guān)系使之在動(dòng)力與熱力上協(xié)調(diào)起來(lái)，以求得更為合理的初始場(chǎng)；或者用來(lái)更新預(yù)報(bào)值，以得到質(zhì)量更好更為合理的初始場(chǎng)；或者用來(lái)更新預(yù)報(bào)值，以得到質(zhì)量更好的連續(xù)天氣分析及預(yù)報(bào)結(jié)果的問(wèn)題。的連續(xù)天氣分析及預(yù)報(bào)結(jié)果的問(wèn)題。資料同化（資料同化（data assimilation)： 大氣資料同化；海洋資料同化；陸面資料同化 多項(xiàng)式法；逐步訂正法；多元最優(yōu)插值法；變分同化法；集合卡爾曼濾波法陳耀登75陳耀登75陳耀登參數(shù)化、初邊界條件客觀分析-起源 分析：在天氣預(yù)報(bào)中專指“實(shí)況天氣圖的繪畫”畫 天氣圖。分析是靠手工進(jìn)行的, 不客觀。兩個(gè)預(yù)報(bào)員畫同一張?zhí)鞖鈭D,由于學(xué)歷不同、經(jīng)驗(yàn)不同、動(dòng)手能力和敬業(yè)精神不同,他們會(huì)得到很不相同的結(jié)果。陳耀登76陳耀登76陳耀登參數(shù)化、初邊界條件客觀分析-定義 定義：使用滿足某一