• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    建立我國數(shù)值預(yù)報同化診斷系統(tǒng)的行動方案頂層設(shè)計

    2023-05-11 02:50:52朱國富
    關(guān)鍵詞:數(shù)值觀測資料

    朱國富

    (中國氣象局地球系統(tǒng)數(shù)值預(yù)報中心,北京 100081)

    0 引言

    數(shù)值天氣預(yù)報已成為公認(rèn)的現(xiàn)代氣象事業(yè)及其發(fā)展的最核心的科技。資料同化是數(shù)值天氣預(yù)報的一個主要組成部件,它基于統(tǒng)計最優(yōu)估計方法產(chǎn)生大氣狀態(tài)的最優(yōu)估計,來提供預(yù)報模式的盡可能準(zhǔn)確的初始條件。更精確的資料同化方法的使用,諸如“四維變分同化方法”和“同化系統(tǒng)誤差參數(shù)設(shè)定的優(yōu)化改進(jìn)”這些主要科學(xué)進(jìn)展,連同“可利用觀測資料的大量增加”,使得預(yù)報模式的初始條件得以改善,真正帶來了數(shù)值預(yù)報性能的顯著提高[1]。

    歷經(jīng)多年,數(shù)值天氣預(yù)報中資料同化系統(tǒng)已經(jīng)發(fā)展成非常復(fù)雜的系統(tǒng)?,F(xiàn)在的同化系統(tǒng)要處理非常廣泛的空基、地基多種來源的氣象觀測數(shù)據(jù),將這些數(shù)據(jù)與大氣狀態(tài)的先驗(yàn)信息結(jié)合起來,并使用綜合復(fù)雜的預(yù)報模式線性化模型進(jìn)行四維資料同化以確保觀測信息得到與模式動力一致的傳播和吸納。對如此復(fù)雜的系統(tǒng),進(jìn)行有效和協(xié)調(diào)的診斷是必要的,例如,僅就被同化的觀測數(shù)量,近幾十年來呈“指數(shù)級”增長,未來還將繼續(xù)增長;這需要基本的診斷工具進(jìn)行一整套日常例行的觀測信息監(jiān)控來保障正常的信息輸入,并需要復(fù)雜先進(jìn)的診斷工具評估數(shù)據(jù)質(zhì)量和信息冗余,以便對最有希望的觀測進(jìn)行優(yōu)先排序處理。同化診斷中觀測信息監(jiān)控的一個經(jīng)典例子是,利用1981年月平均資料,歐洲中期天氣預(yù)報中心(ECMWF)同化系統(tǒng)發(fā)現(xiàn),在距離開普敦東南2300 km的一個南印度洋孤島(南非Marion Island)上,背景風(fēng)和觀測報告風(fēng)之間平均有10°~12°的系統(tǒng)性差異。對此進(jìn)行的一項調(diào)查顯示,多年前安裝計算風(fēng)向的軟件時,指定的風(fēng)向參考是磁北極,而不是地理北極!之后ECMWF在1982年開啟了一個業(yè)務(wù)和研究部門間的重點(diǎn)合作研究項目,發(fā)現(xiàn)一些觀測平臺在持續(xù)地發(fā)送不正確的天氣報告,甚至有些簡直不按照約定標(biāo)準(zhǔn)編寫數(shù)據(jù)代碼!1985年,ECMWF被世界氣象組織(WMO)指定為監(jiān)測全球高空資料的領(lǐng)導(dǎo)中心[2]。這是ECMWF之所以成為領(lǐng)導(dǎo)中心地位的一個具體經(jīng)歷??梢?,也正是源自觀測信息監(jiān)控,可以理解對分析場產(chǎn)生有害影響的觀測黑名單是作為基礎(chǔ)一環(huán)的質(zhì)量控制方式。

    診斷一直是ECMWF研究和開發(fā)的一個活躍領(lǐng)域。而且,數(shù)值預(yù)報系統(tǒng)的診斷是ECMWF總體戰(zhàn)略所必需的組成部分:面向ECMWF發(fā)展“從此處向何處去”,一個診斷戰(zhàn)略是必不可少的[3]。

    診斷工作不僅僅是獲得一個總體了解,因?yàn)樵\斷是面向數(shù)值預(yù)報系統(tǒng)的所有環(huán)節(jié)及其組成部分,所有的研究人員都需要在自己的特定領(lǐng)域中制作相應(yīng)的詳盡診斷。另一方面,隨著數(shù)值預(yù)報的發(fā)展,不可避免地,制作出預(yù)報的所有環(huán)節(jié)任務(wù)更加細(xì)分化和專門化。正是如此,提出了一個問題:未來我們是否需要加強(qiáng)跨傳統(tǒng)部門界線的溝通、協(xié)調(diào)和協(xié)作。面對這一問題,2010年ECMWF成立了“診斷工作組”(WGD:Working Group on Diagnostics),并作為WGD的第一個主要成果撰寫出技術(shù)報告“診斷研究進(jìn)展”(Developments in diagnostics research)[3]。這份報告為考慮一般意義的診斷策略提供了借鑒,并為從無到有地開展常態(tài)化同化系統(tǒng)診斷工作提供了參考,包括診斷的要義、診斷工作定位的角色作用及其實(shí)現(xiàn)途徑或方式、診斷工作內(nèi)容的框架等基本方面,例如數(shù)值預(yù)報的診斷框架(圖1)。

    從圖1可以一般地描述數(shù)值預(yù)報意義上診斷工作的目的:深入動力框架、物理過程、資料同化各個組成部分,通過識別系統(tǒng)的缺陷、減小各環(huán)節(jié)的偏差、改進(jìn)物理過程不確定性和觀測不確定性的描述等工作,來識別、調(diào)查和理解數(shù)值模式系統(tǒng)的預(yù)報性能問題和模式氣候狀態(tài),以助力數(shù)值預(yù)報系統(tǒng)發(fā)展戰(zhàn)略的實(shí)現(xiàn);并在所涉及的工作中將研究和業(yè)務(wù)的各個部門聯(lián)合起來。

    后來ECMWF的診斷工作[4],其目的不變。診斷工作的發(fā)展體現(xiàn)在以下方面:診斷的施展空間面向預(yù)報系統(tǒng)在“所有預(yù)報時效”敏銳度和可靠性的改進(jìn);診斷以“極端天氣預(yù)報”為目標(biāo);與診斷關(guān)聯(lián)的互動超越ECMWF內(nèi)部的研究和業(yè)務(wù)各部門而擴(kuò)大到ECMWF成員國之間。而且,預(yù)報系統(tǒng)(以及觀測)的準(zhǔn)確性不斷提高,意味著誤差比以往任何時候都要小,需要發(fā)展更精確的診斷工具來量化這些誤差。

    本文將面向中國氣象局(CMA)四維變分同化系統(tǒng)自主研發(fā)的現(xiàn)實(shí)狀況和發(fā)展需求,思考國際上氣象數(shù)值預(yù)報先進(jìn)業(yè)務(wù)中心的相關(guān)工作,結(jié)合在CMA同化系統(tǒng)自主研發(fā)實(shí)踐中的理解,形成建立CMA常態(tài)化同化診斷系統(tǒng)的頂層設(shè)計,旨在初步形成一套日常例行的同化診斷能力。

    1 ECMWF診斷工作的啟示

    1.1 診斷的要義(是客觀要求的“問題導(dǎo)向”,不是主觀的“找事情”;是“理解”,不是“較真”)

    對于不同工作和人群的視角,“診斷”看似有不同的解釋,但所有的解釋都體現(xiàn)了“理解”這一概念。像ECMWF這樣的業(yè)務(wù)預(yù)報中心的診斷工作屬于“對(預(yù)報系統(tǒng))誤差的診斷”(diagnosis of error),是調(diào)查、發(fā)現(xiàn)和理解數(shù)值預(yù)報系統(tǒng)中觀測、分析或預(yù)報各環(huán)節(jié)的誤差,進(jìn)而解決問題來推動數(shù)值預(yù)報的發(fā)展。ECMWF診斷工作組從成立之初的2010年到2019年,盡管具體的診斷工作和診斷工具有了發(fā)展更新,但不變的是診斷要義:理解[2-3]。這對開展診斷工作的啟示是:數(shù)值預(yù)報系統(tǒng)的診斷是通過“問題導(dǎo)向”來調(diào)查研究制作出預(yù)報的各環(huán)節(jié)的誤差,和通過“理解”來推動數(shù)值預(yù)報的發(fā)展,也就是它“基于問題、歸于理解和終于協(xié)作”來解決問題。

    事實(shí)上,這對于我國CMA數(shù)值模式的自主研發(fā)更有著現(xiàn)實(shí)意義。自主研發(fā)基于理解,理解是科技自主自強(qiáng)的首要要求,是氣象核心科技數(shù)值預(yù)報高質(zhì)量發(fā)展的前提條件,不追求理解的自主研發(fā)是不可思議的。注重理解也是科學(xué)家精神的基本內(nèi)涵。

    1.2 認(rèn)識診斷工作的出發(fā)點(diǎn)(首要在“前端”和“基本”)

    在ECMWF最初成立的診斷工作組的診斷研究進(jìn)展[3]中,觀測信息監(jiān)控(Observation monitoring)是ECMWF診斷系統(tǒng)的一個關(guān)鍵組成部分;其中的許多診斷是基于觀測空間下背景場(B)和分析場(A)與觀測(O)的離差,即所謂的(O-B)和(O-A)。對于再分析的診斷,重點(diǎn)是數(shù)據(jù)的時間一致性。ECMWF集合預(yù)報系統(tǒng)的診斷集中關(guān)注“集合離散度和集合平均誤差之間的不一致性”如何與“初始條件和預(yù)報模式中不確定性原因的描述缺陷”關(guān)聯(lián)起來。從中我們能聯(lián)系到其中內(nèi)在的原因和事實(shí):

    1)對于觀測信息監(jiān)控,許多診斷基于(O-B)和(O-A)。這是因?yàn)閷Γ∣-B)和(O-A)進(jìn)行的數(shù)據(jù)統(tǒng)計為觀測、同化或預(yù)報等方面的特征描述提供了強(qiáng)有力的基礎(chǔ),而且其中的(O-B)(稱作“新息”(innovation))是同化系統(tǒng)的基本輸入,是用來更新背景場而得到分析場的、用作訂正的新信息;如果基本輸入錯誤,則直接對分析場產(chǎn)生有害影響,如引言中的觀測信息監(jiān)控的一個經(jīng)典例子[2]。

    2)對于再分析,它是基于觀測、同化、預(yù)報三個串聯(lián)方面的先進(jìn)性和一致性來得到先進(jìn)的穩(wěn)定一致的長時間序列格點(diǎn)資料。其中觀測數(shù)據(jù)的時間一致性不僅處在這個串聯(lián)鏈條的最前端,而且是再分析的關(guān)鍵性重點(diǎn);因?yàn)樵俜治鍪腔谑褂卯?dāng)下最先進(jìn)的業(yè)務(wù)同化系統(tǒng)和預(yù)報模式系統(tǒng)而制作出來的,這樣再分析中同化和預(yù)報這兩個方面的先進(jìn)性和一致性而自然得到了保障。

    3)對于集合預(yù)報,評估集合預(yù)報可靠性的一個通用診斷是集合離散度和集合平均誤差之間的關(guān)系[5]?!案怕暑A(yù)報崩潰”被定義為當(dāng)集合離散度很小、但集合平均誤差很大時發(fā)生的情況。因此集合離散度和集合平均誤差是集合預(yù)報診斷中必需首要考慮的“基本所在”。

    這些對我們開展診斷工作的啟示是:診斷工作首要和首先從“鏈條前端”和“基本所在”開始。為此,需要在認(rèn)知上理解清楚“哪些是基本的”和在實(shí)施上考慮哪些是“現(xiàn)在能做的”和“客觀上必須做的”。

    1.3 診斷工作的架構(gòu)和范圍(所有環(huán)節(jié)及其組成部分)

    數(shù)值預(yù)報的步驟過程表現(xiàn)為:通過資料同化方法得到分析場,由它作為預(yù)報模式的初值;然后積分預(yù)報模式到未來某一時刻,便得到該時刻的預(yù)報場。相應(yīng)地,數(shù)值預(yù)報系統(tǒng)有兩個主要組成部分:初值形成的資料同化系統(tǒng)和時間積分的預(yù)報模式系統(tǒng)。資料同化系統(tǒng)包含引入各類觀測的若干觀測算子和作用于每一個觀測來傳播觀測信息的一個核心框架;預(yù)報模式系統(tǒng)包含一個動力框架和若干物理過程。通常地,預(yù)報模式動力框架是指預(yù)報模式方程的絕熱、無摩擦的網(wǎng)格尺度動力過程,各個模式物理過程是指不同的源、匯項或不同的次網(wǎng)格尺度物理過程(如深/淺對流過程、輻射過程、大氣邊界層過程和陸面過程)。

    診斷是面向數(shù)值預(yù)報系統(tǒng)的各個環(huán)節(jié)及其組成部分。因此,數(shù)值預(yù)報系統(tǒng)診斷工作的架構(gòu)包括資料同化診斷和預(yù)報模式診斷兩個部分。預(yù)報模式的診斷范圍包括動力框架診斷和各個物理過程診斷。

    因?yàn)橘Y料同化是上接各類觀測、下連預(yù)報模式的一個綜采工作面(coal-face),因此它的診斷范圍是巨大的,而且隨著觀測日益增加在不斷增長擴(kuò)大。資料同化的診斷范圍包括面向觀測(O)的觀測信息監(jiān)控、面向(O-B)和(O-A)統(tǒng)計特征描述的診斷及其觀測質(zhì)量的監(jiān)測和評估、面向得到分析場之求解過程的監(jiān)控和分析、面向同化系統(tǒng)預(yù)先設(shè)定的誤差參數(shù)的診斷、面向分析場(A)的觀測對分析的敏感性評估、面向模式預(yù)報效果的觀測對預(yù)報質(zhì)量的影響評估等。

    1.4 診斷工作組的角色定位(把各個部門聯(lián)合起來,加強(qiáng)跨部門的“交流”“協(xié)調(diào)”和“協(xié)作”)

    隨著數(shù)值預(yù)報的發(fā)展,生產(chǎn)出預(yù)報的所有環(huán)節(jié)任務(wù)更加細(xì)分化和專門化;同時,把各個部門聯(lián)合起來的作用(the unifying role)又隨之日益變得越來越重要。正是在這個背景下,ECMWF成立了診斷工作組(WGD),來加強(qiáng)跨部門的交流、協(xié)調(diào)和合作[3]。WGD的第一個任務(wù)就是明確界定這些角色作用:

    1)合作課題的監(jiān)督(the over-sight of collaborative projects)。一些很好的例子如:現(xiàn)有同化或預(yù)報問題需要協(xié)作開展的診斷;預(yù)報系統(tǒng)升級時所需的各環(huán)節(jié)診斷評估;改善災(zāi)害性天氣預(yù)報的交叉性診斷課題(涉及該類天氣的識別方法,預(yù)報系統(tǒng)對其預(yù)報能力相關(guān)的所用資料的同化和關(guān)鍵物理過程的描述等評估,集合預(yù)報系統(tǒng)是否能夠足夠地表示其發(fā)生的可能性等調(diào)查研究)。除了這些明確主題的合作課題,還可能從中產(chǎn)生未來的診斷領(lǐng)域或開發(fā)新的診斷工具。

    2)診斷工具研發(fā)的戰(zhàn)略協(xié)調(diào)(the strategic coordination of diagnostic developments)。包括:共同感興趣的新診斷工具的研發(fā),如對于物理過程診斷工具的研發(fā),因?yàn)槲磥黼S著模式變得更好,物理過程的改換對模式氣候的影響應(yīng)該會更小并需要更精確的診斷工具,以及物理過程數(shù)量的增加大大增加了相互作用的范圍,這些因素都要求新診斷工具的研發(fā)及相關(guān)環(huán)節(jié)的研發(fā)協(xié)作;使現(xiàn)有診斷工具得到更廣泛應(yīng)用的戰(zhàn)略協(xié)同,如診斷工具的通用腳本語言的使用和更廣泛可操作的界面接口,以及診斷工具的持續(xù)維護(hù)和文檔化;需要確保診斷在計算和存儲資源方面得到滿足。

    正是通過這些角色作用,ECMWF診斷工作組利用診斷研究,理解業(yè)務(wù)數(shù)值預(yù)報的優(yōu)勢與不足,發(fā)現(xiàn)問題并研討確定解決問題的方案,確立了它在模式進(jìn)展中的中心地位,推動數(shù)值預(yù)報系統(tǒng)多年來取得了許多重要的進(jìn)展[6]。

    這些對我們開展診斷工作的啟示是:把各個部門聯(lián)合起來的“融合發(fā)展”是診斷工作組的工作屬性;“能協(xié)作得起來”的客觀要求是數(shù)據(jù)存取格式和圖表繪制顯示的通用性、共用性,而且特別是高度的靈活性,因?yàn)樗鼈兪窃\斷工作的基礎(chǔ)性支撐:診斷中為了發(fā)現(xiàn)、分析問題,需要大量的數(shù)據(jù)和圖表,來以此為基礎(chǔ)進(jìn)行各種各樣的對照比較和數(shù)據(jù)統(tǒng)計。

    上述合作、協(xié)調(diào)和交流的角色作用對我們建立診斷工作組以及如何發(fā)揮作用和從無到有地開展診斷工作提供了直接參考。事實(shí)上,在ECMWF專家的一直關(guān)注下[6],中國氣象局地球系統(tǒng)數(shù)值預(yù)報中心的一個實(shí)體診斷小組已經(jīng)成立;ECMWF專家還建議:應(yīng)該考慮制定一項戰(zhàn)略來確定診斷小組的主要優(yōu)先任務(wù),包括建立一整套日常例行的診斷,旨在監(jiān)測和指導(dǎo)數(shù)值預(yù)報系統(tǒng)的改進(jìn)。

    2 CMA常態(tài)化同化診斷系統(tǒng)的頂層設(shè)計

    CMA同化系統(tǒng)采用四維變分同化方案。變分同化方法是大氣狀態(tài)的后驗(yàn)概率最大時取眾數(shù)的統(tǒng)計最優(yōu)估計方法;導(dǎo)出變分同化目標(biāo)泛函的一個基本假定是“觀測和背景場的誤差無偏且為正態(tài)分布”。CMA四維變分同化系統(tǒng)在業(yè)務(wù)實(shí)施上利用一種增量方法[7],考慮內(nèi)外循環(huán),進(jìn)行目標(biāo)泛函的極小化迭代數(shù)值求解,來減少計算開銷。

    CMA同化診斷系統(tǒng)面向CMA同化系統(tǒng)的各個組成部分和業(yè)務(wù)同化系統(tǒng)研發(fā)的環(huán)節(jié)步驟,理清“現(xiàn)在能做的”和“客觀上必須做的”基本所在,以程序流的全過程和業(yè)務(wù)研發(fā)的內(nèi)在邏輯進(jìn)行頂層設(shè)計。

    2.1 同化診斷的總體架構(gòu)

    首先理清楚任一同化系統(tǒng)運(yùn)行的組成部分和業(yè)務(wù)同化系統(tǒng)研發(fā)的基本環(huán)節(jié)步驟:

    1)同化系統(tǒng)中的程序流表現(xiàn)為觀測輸入、求解過程和結(jié)果輸出。

    2)現(xiàn)行的先進(jìn)同化方法都屬于統(tǒng)計最優(yōu)估計方法,因?yàn)樗鼈兓跀?shù)據(jù)誤差是高斯分布等假定,所以都存在預(yù)先設(shè)定的參數(shù);對于一個業(yè)務(wù)的同化系統(tǒng),在其建立之后,不論是必需的不斷優(yōu)化改進(jìn),還是要考慮系統(tǒng)中參數(shù)的季節(jié)變化,都首先是需要診斷預(yù)先設(shè)定的參數(shù)是否合適。

    3)對于不斷涌現(xiàn)和更替的新資料,是否能進(jìn)入業(yè)務(wù)數(shù)值預(yù)報系統(tǒng),必需對新資料的準(zhǔn)入進(jìn)行診斷評估;面向觀測—預(yù)報交互式數(shù)值預(yù)報,需要診斷評估觀測對分析或預(yù)報的敏感性影響。這些都是基于同化系統(tǒng)的診斷。

    1.1.3 CKD-MBD診斷標(biāo)準(zhǔn) 參照2013年中華醫(yī)學(xué)會腎臟病學(xué)分會頒布《慢性腎臟病礦物質(zhì)和骨異常診治指導(dǎo)》[8]擬定:(1)血清鈣(Ca)、磷(P)、甲狀旁腺激素(iPTH)或維生素D代謝異常;(2)骨轉(zhuǎn)化、礦化、骨量、骨線性生長或骨強(qiáng)度異常;(3)血管或其他軟組織鈣化。符合以上三點(diǎn)之一即可診斷為CKD-MBD。

    由此按照同化系統(tǒng)的程序流和業(yè)務(wù)同化系統(tǒng)研發(fā)的內(nèi)在邏輯,形成同化診斷的總體架構(gòu):

    1)利用同化系統(tǒng)的“輸入信息”進(jìn)行診斷。

    2)利用變分同化的“執(zhí)行過程”進(jìn)行診斷。

    3)利用同化系統(tǒng)的“輸出結(jié)果”進(jìn)行診斷。

    4)面向業(yè)務(wù)同化系統(tǒng)研發(fā)精雕細(xì)刻的“參數(shù)優(yōu)化”進(jìn)行診斷。

    5)面向業(yè)務(wù)同化系統(tǒng)的新觀測資料準(zhǔn)入進(jìn)行診斷評估。

    這些在一起能夠初步形成一套日常例行的同化診斷能力。

    2.2 同化診斷的內(nèi)容及其途徑和作用

    2.2.1 輸入信息的診斷

    這在內(nèi)容上屬于觀測信息監(jiān)控。觀測信息監(jiān)控處在整個數(shù)值預(yù)報系統(tǒng)診斷的前端,也是ECMWF診斷系統(tǒng)的一個關(guān)鍵組成部分。

    觀測信息監(jiān)控包含了觀測資料計數(shù)及其空間分布的監(jiān)控和輸入信息的診斷。資料計數(shù)包括到報的資料數(shù)、被篩除的資料數(shù)、被同化的資料數(shù),空間分布包括各類觀測的總體和不同層(或通道)的空間分布。輸入信息的許多診斷是基于觀測空間下(O-B),如它的統(tǒng)計特征描述,包括均值、均方差、偏差訂正、頻數(shù)分布(對應(yīng)概率密度分布PDF)等。利用(O-B)的信息,在了解數(shù)據(jù)、認(rèn)識數(shù)據(jù)上有諸如以下的許多作用:

    1)用來監(jiān)控資料同化中觀測亦或背景場的誤差是否無偏和正態(tài)分布,以及監(jiān)控觀測資料處理的偏差訂正。

    2)它的時序資料可用來監(jiān)控觀測資料的時間穩(wěn)定性、質(zhì)量控制有效性和觀測質(zhì)量,確保只有始終如一好質(zhì)量的資料在業(yè)務(wù)中被同化使用。

    3)利用它和其他觀測進(jìn)行比較,來監(jiān)控觀測的質(zhì)量。

    事實(shí)上,觀測空間下(O-B)和(O-A)數(shù)據(jù)統(tǒng)計為觀測、同化或預(yù)報等方面的特征描述提供了強(qiáng)有力的基礎(chǔ),許多診斷是基于此,因?yàn)椋∣-B)和(O-A)的監(jiān)控對資料同化過程的上、下游都發(fā)揮作用,除了上面的作用,如對上游:用于同化系統(tǒng)中誤差參數(shù)的預(yù)先設(shè)定即下面2.2.4節(jié)“參數(shù)優(yōu)化的診斷”;如對下游:當(dāng)不同觀測都一致發(fā)現(xiàn)(O-B)的系統(tǒng)性的非零值時可突顯出預(yù)報模式問題。

    2.2.2 執(zhí)行過程的診斷

    這是面向得到分析場的求解過程進(jìn)行的診斷。極小化求解迭代是否收斂和收斂程度,分析值受其直接影響。所以一個基本的方面是需要通過極小化的迭代步數(shù)、目標(biāo)泛函值下降率、目標(biāo)泛函梯度下降程度等的監(jiān)控和分析,來診斷目標(biāo)泛函求解的收斂情況。

    2.2.3 輸出結(jié)果的診斷

    這是面向分析結(jié)果進(jìn)行觀測對分析的敏感性評估。通常地,根據(jù)分析場(xa)與觀測(yo)和背景場(xb)的關(guān)系:xa=xb+W[yo-H(xb)](其中,H為觀測算子,W為最優(yōu)權(quán)重矩陣),利用變分同化系統(tǒng)核心框架中的變量變換和同化系統(tǒng)觀測算子的切線性模式及其伴隨模式,通過伴隨技術(shù)可以相當(dāng)嚴(yán)謹(jǐn)?shù)氐玫礁鱾€觀測對分析場的敏感性。還可以比較簡單地利用變分同化系統(tǒng)目標(biāo)泛函的觀測項(Jo),通過由(O-B)和(O-A)計算的目標(biāo)泛函值來大致地評估各個觀測對分析場的相對貢獻(xiàn)。

    面向分析場的診斷可用來評估不同觀測的貢獻(xiàn)作用多少、是否信息冗余和可能的數(shù)據(jù)質(zhì)量,以便對最有希望的觀測進(jìn)行優(yōu)先排序處理和應(yīng)用于觀測—預(yù)報交互式數(shù)值預(yù)報。

    2.2.4 參數(shù)優(yōu)化的診斷

    同化系統(tǒng)誤差參數(shù)包含于xa=xb+W[yo-H(xb)]中的最優(yōu)權(quán)重矩陣W,所以直接左右了分析結(jié)果xa。因此,診斷預(yù)先設(shè)定的誤差參數(shù)是否合適是同化系統(tǒng)自主研發(fā)在系統(tǒng)建立之后需要考慮的固有步驟,以及應(yīng)用已有或他人同化系統(tǒng)時進(jìn)行參數(shù)優(yōu)化改進(jìn)的必備環(huán)節(jié)。

    可以利用(O-B)和(O-A),通過觀測空間的后驗(yàn)診斷方法[8-9]進(jìn)行新息、背景場誤差、觀測誤差的一致性診斷和分析誤差的診斷:,R表示觀測誤差協(xié)方差矩陣,B表示背景場誤差協(xié)方差矩陣,H表示H的切線性模式,上角標(biāo)T表示矩陣(或向量)的轉(zhuǎn)置,HBHT為觀測空間下的背景場誤差協(xié)方差,HAHT為觀測空間下的分析場誤差協(xié)方差)。

    通過計算出的后驗(yàn)診斷值與同化系統(tǒng)中的預(yù)先設(shè)定值的比較,用來診斷同化系統(tǒng)中預(yù)先設(shè)定的誤差參數(shù)是否合適,以及如果不合適時給出調(diào)優(yōu)改進(jìn)的方向。

    2.2.5 新種類觀測資料準(zhǔn)入的診斷評估

    這是面向模式預(yù)報進(jìn)行觀測對預(yù)報的影響評估。通常地,對現(xiàn)有的所有觀測,通過觀測系統(tǒng)試驗(yàn)(OSE)(考慮逐個觀測排除的一個時段預(yù)報試驗(yàn)(data denial experiments),即逐個觀測排除敏感試驗(yàn)),來診斷評估不同設(shè)置的各種類觀測對預(yù)報質(zhì)量的影響。對于新種類觀測,通過一體化設(shè)計的觀測系統(tǒng)模擬試驗(yàn)(OSSE)和OSE,來評估新種類觀測資料是否準(zhǔn)入。正是這樣的診斷評估保障了使用所有資料的控制預(yù)報試驗(yàn)的預(yù)報質(zhì)量最好。

    3 開展這項工作的主要優(yōu)先任務(wù)和開端性思考

    3.1 主要優(yōu)先任務(wù)

    巧婦難為無米之炊。同化診斷工作的前提和基礎(chǔ)是必需同化診斷必備數(shù)據(jù)的數(shù)據(jù)庫建設(shè)。首要是觀測空間的數(shù)據(jù),其最基本內(nèi)容包含:觀測值yo、觀測誤差標(biāo)準(zhǔn)差σo、質(zhì)量控制標(biāo)識IQC、觀測偏差bo、背景場的觀測相當(dāng)量yb、分析場的觀測相當(dāng)量ya。也由此得到(O-B)和(O-A)數(shù)據(jù)。

    必須用利于通用性、共用性、特別是高度靈活性的觀測空間數(shù)據(jù)庫和模式空間數(shù)據(jù)庫的數(shù)據(jù)存取格式;它們是診斷工作的最基礎(chǔ)性工作和保障性支撐,也是保障“能協(xié)作得起來”的客觀要求。

    在同化診斷工作的內(nèi)容上,首要是針對觀測資料數(shù)統(tǒng)計和利用觀測空間的(O-B)和(O-A)數(shù)據(jù),建立一套日常例行的同化診斷系統(tǒng),來實(shí)現(xiàn)觀測信息監(jiān)控、觀測項(Jo)中觀測對分析作用大小的診斷評估和同化系統(tǒng)參數(shù)的觀測空間后驗(yàn)診斷。

    3.2 開端性思考

    診斷是“基于問題、歸于理解和終于協(xié)作”來解決問題而推動數(shù)值預(yù)報系統(tǒng)的進(jìn)展。診斷的出發(fā)點(diǎn)在“問題導(dǎo)向”,且是以過程為導(dǎo)向的面向各個組成環(huán)節(jié)的問題,診斷的要求是“追求理解”,完成診斷工作依靠“團(tuán)結(jié)協(xié)作”。這些是診斷工作的客觀固有的工作屬性,是診斷工作的應(yīng)有之義和實(shí)際樣子。特別是,把各個部門聯(lián)合起來的“融合發(fā)展”是診斷工作組的工作屬性。開啟診斷工作在觀念上首先需要這些基本共識,來了解這是一份在作風(fēng)上不在說好話、在過程上還是個急不來的工作,而營造和增加工作上的理解氛圍和克服“找事情”“較真”“各自為戰(zhàn)”的主觀性傾向。

    在行動上特別需要把各個部門聯(lián)合起來的“融合發(fā)展”來開展同化診斷工作,例如:

    1)對于不斷更替的衛(wèi)星輻射率新資料,衛(wèi)星觀測信息監(jiān)控的研發(fā)需要同化診斷部門和衛(wèi)星資料同化部門的共同協(xié)作。

    2)對于極小化求解過程的診斷和利用觀測算子的切線性模式及其伴隨模式的觀測對分析或預(yù)報的敏感性影響的診斷評估,需要同化診斷部門和同化系統(tǒng)研發(fā)部門的共同協(xié)作。

    而且,貫徹落實(shí)“充分發(fā)揮組織優(yōu)勢,強(qiáng)化與上下游氣象業(yè)務(wù)單位的合作銜接,統(tǒng)籌用好全國資源”①http://www.cma.gov.cn/2011xwzx/2011xqxxw/2011xqxyw/202203/t20220303_592682.html。,來優(yōu)勢互補(bǔ)、團(tuán)結(jié)協(xié)作地推進(jìn)同化診斷系統(tǒng)的技術(shù)攻關(guān)和業(yè)務(wù)建設(shè),例如:

    1)可以通過CMA數(shù)值預(yù)報中心與國家衛(wèi)星氣象中心、國家氣象信息中心、CMA氣象探測中心的通力協(xié)作,針對數(shù)據(jù)管理,建立通用、共用、特別是高度靈活的觀測空間和模式空間診斷必備數(shù)據(jù)庫存取平臺,筑牢診斷工作的前端基礎(chǔ)。

    2)可以通過CMA數(shù)值預(yù)報中心與國家氣象中心、國家衛(wèi)星氣象中心、華風(fēng)氣象傳媒集團(tuán)有限責(zé)任公司的通力協(xié)作,針對軟件系統(tǒng)的開發(fā)和管理,研討一體化、模塊化的軟件系統(tǒng)封裝技術(shù)、界面接口頂層設(shè)計和避免自說自話的開發(fā)工作方式,建立靈活、豐富的數(shù)值預(yù)報診斷輸出圖表繪制顯示平臺,在技術(shù)上支撐診斷工作的現(xiàn)實(shí)需要和規(guī)范化發(fā)展。

    3)預(yù)報員有著在數(shù)值預(yù)報產(chǎn)品應(yīng)用中通過天氣場景分類的模式產(chǎn)品表現(xiàn)而更具象化地發(fā)現(xiàn)問題的優(yōu)勢,研發(fā)者有著在數(shù)值預(yù)報系統(tǒng)診斷中通過物理和數(shù)值模型的理解而更深入地發(fā)現(xiàn)問題的優(yōu)勢?;趩栴}導(dǎo)向,建立和加強(qiáng)一線研發(fā)人員和全國有經(jīng)驗(yàn)的優(yōu)秀預(yù)報員的聯(lián)系橋梁,通過課題合作機(jī)制和定期交流機(jī)制,共贏地促進(jìn)研發(fā)人員對數(shù)值模式問題的診斷與改進(jìn)和預(yù)報員對數(shù)值模式產(chǎn)品的理解與更好應(yīng)用。

    對于開啟同化系統(tǒng)的診斷工作,觀念上有共識至關(guān)重要,是團(tuán)結(jié)奮斗的思想基礎(chǔ)和前提;行動上通用、共用、特別是高度靈活性的診斷必備數(shù)據(jù)至關(guān)重要,是開展診斷工作的前端基礎(chǔ)任務(wù)和團(tuán)結(jié)協(xié)作的支撐條件。需要遵循診斷工作的實(shí)際樣子,開啟和推動同化診斷系統(tǒng)研發(fā)基于“問題導(dǎo)向、科學(xué)理解和團(tuán)結(jié)協(xié)作”的高質(zhì)量發(fā)展。

    猜你喜歡
    數(shù)值觀測資料
    觀測到恒星死亡瞬間
    軍事文摘(2023年18期)2023-11-03 09:45:42
    用固定數(shù)值計算
    數(shù)值大小比較“招招鮮”
    Party Time
    JUST A THOUGHT
    PAIRS & TWOS
    天測與測地VLBI 測地站周圍地形觀測遮掩的討論
    可觀測宇宙
    太空探索(2016年7期)2016-07-10 12:10:15
    基于Fluent的GTAW數(shù)值模擬
    焊接(2016年2期)2016-02-27 13:01:02
    高分辨率對地觀測系統(tǒng)
    太空探索(2015年8期)2015-07-18 11:04:44
    亚洲欧美激情综合另类| 国产精品永久免费网站| 亚洲男人天堂网一区| 欧美一级a爱片免费观看看 | 免费在线观看完整版高清| 亚洲一卡2卡3卡4卡5卡精品中文| 哪里可以看免费的av片| 亚洲精品在线美女| 午夜日韩欧美国产| 免费观看精品视频网站| 亚洲av片天天在线观看| 真人做人爱边吃奶动态| 精品久久久久久久人妻蜜臀av| 欧美黄色片欧美黄色片| 亚洲精品中文字幕在线视频| 免费在线观看成人毛片| 91成人精品电影| 久久久久久久久免费视频了| 久久中文字幕人妻熟女| 久久婷婷成人综合色麻豆| 国内久久婷婷六月综合欲色啪| 免费看十八禁软件| 日韩精品免费视频一区二区三区| 久久久久免费精品人妻一区二区 | 99国产综合亚洲精品| 中文字幕最新亚洲高清| 亚洲美女黄片视频| av欧美777| 制服诱惑二区| 黄网站色视频无遮挡免费观看| 88av欧美| 欧美+亚洲+日韩+国产| 看片在线看免费视频| 两性夫妻黄色片| 亚洲中文av在线| 亚洲av成人一区二区三| 国产精品九九99| 国内少妇人妻偷人精品xxx网站 | 可以免费在线观看a视频的电影网站| 亚洲av电影在线进入| 国产又色又爽无遮挡免费看| 麻豆成人av在线观看| 国内久久婷婷六月综合欲色啪| 人成视频在线观看免费观看| 91麻豆精品激情在线观看国产| 999久久久精品免费观看国产| 在线观看舔阴道视频| 欧美日韩福利视频一区二区| 免费看十八禁软件| 欧美久久黑人一区二区| av视频在线观看入口| 99精品欧美一区二区三区四区| 麻豆成人av在线观看| www.熟女人妻精品国产| av超薄肉色丝袜交足视频| 亚洲欧美日韩无卡精品| 日本成人三级电影网站| 首页视频小说图片口味搜索| 午夜福利在线观看吧| 国产一区二区在线av高清观看| 久久人人精品亚洲av| 精品国产国语对白av| 亚洲一区高清亚洲精品| 18禁裸乳无遮挡免费网站照片 | 亚洲精品国产一区二区精华液| 久久欧美精品欧美久久欧美| 色综合欧美亚洲国产小说| 亚洲五月婷婷丁香| 久久精品国产亚洲av香蕉五月| 久久久久久久久免费视频了| 国产麻豆成人av免费视频| 亚洲九九香蕉| 国产一级毛片七仙女欲春2 | 成人18禁在线播放| 老鸭窝网址在线观看| 91成年电影在线观看| 国产免费av片在线观看野外av| 777久久人妻少妇嫩草av网站| 老司机福利观看| 国产伦在线观看视频一区| 99久久综合精品五月天人人| 岛国视频午夜一区免费看| 国产爱豆传媒在线观看 | 国产成人精品无人区| 露出奶头的视频| 久久香蕉国产精品| 亚洲无线在线观看| 精品久久久久久久毛片微露脸| 99精品欧美一区二区三区四区| 国产免费av片在线观看野外av| 深夜精品福利| 18禁黄网站禁片免费观看直播| 巨乳人妻的诱惑在线观看| 国产欧美日韩一区二区三| 国产私拍福利视频在线观看| 免费在线观看完整版高清| www.999成人在线观看| 久久婷婷人人爽人人干人人爱| 亚洲精品色激情综合| 久久天堂一区二区三区四区| 国产精品久久久久久精品电影 | 美女 人体艺术 gogo| 欧美日韩福利视频一区二区| 国产熟女午夜一区二区三区| 中国美女看黄片| 色精品久久人妻99蜜桃| 满18在线观看网站| 欧美成人性av电影在线观看| 亚洲男人天堂网一区| 成在线人永久免费视频| 人成视频在线观看免费观看| a级毛片在线看网站| 午夜精品在线福利| 黄色视频不卡| 国内揄拍国产精品人妻在线 | 久热这里只有精品99| 美女免费视频网站| 黄色 视频免费看| 欧美激情 高清一区二区三区| 俺也久久电影网| 午夜福利免费观看在线| 18禁美女被吸乳视频| av在线播放免费不卡| 亚洲自拍偷在线| 又紧又爽又黄一区二区| 三级毛片av免费| 麻豆成人午夜福利视频| 女警被强在线播放| 中文字幕另类日韩欧美亚洲嫩草| 国产爱豆传媒在线观看 | 久久久精品国产亚洲av高清涩受| 两性夫妻黄色片| 亚洲aⅴ乱码一区二区在线播放 | 久久精品人妻少妇| 亚洲专区中文字幕在线| 精品福利观看| 国产免费男女视频| 亚洲人成网站高清观看| 午夜久久久久精精品| 久久久久国产一级毛片高清牌| 欧美人与性动交α欧美精品济南到| 国产麻豆成人av免费视频| 亚洲成av人片免费观看| 后天国语完整版免费观看| 国产成人影院久久av| 18禁裸乳无遮挡免费网站照片 | 高清在线国产一区| 怎么达到女性高潮| 91在线观看av| 美女国产高潮福利片在线看| 侵犯人妻中文字幕一二三四区| 男人舔女人的私密视频| 村上凉子中文字幕在线| 日韩欧美 国产精品| 久久久久亚洲av毛片大全| 日韩免费av在线播放| www.自偷自拍.com| 精品国产一区二区三区四区第35| 午夜福利一区二区在线看| 深夜精品福利| 久久久久久大精品| 国产亚洲精品久久久久5区| av超薄肉色丝袜交足视频| 搡老妇女老女人老熟妇| 大香蕉久久成人网| 别揉我奶头~嗯~啊~动态视频| 热99re8久久精品国产| 日韩有码中文字幕| 欧美在线一区亚洲| 久久狼人影院| 国内揄拍国产精品人妻在线 | 此物有八面人人有两片| 久久久久久人人人人人| 中文亚洲av片在线观看爽| 久热爱精品视频在线9| 国产亚洲精品av在线| 久热爱精品视频在线9| 精品人妻1区二区| 国产精品免费一区二区三区在线| 成人免费观看视频高清| 又大又爽又粗| bbb黄色大片| 精品久久久久久成人av| 一级毛片高清免费大全| 国产亚洲精品一区二区www| 国产亚洲av高清不卡| 国产又爽黄色视频| 免费高清在线观看日韩| 日本a在线网址| 午夜精品久久久久久毛片777| 欧美成人免费av一区二区三区| 国产激情欧美一区二区| 免费在线观看视频国产中文字幕亚洲| 又紧又爽又黄一区二区| 午夜影院日韩av| 岛国视频午夜一区免费看| 琪琪午夜伦伦电影理论片6080| 高清毛片免费观看视频网站| 桃色一区二区三区在线观看| 99精品在免费线老司机午夜| 白带黄色成豆腐渣| av中文乱码字幕在线| 午夜福利一区二区在线看| 久久中文字幕一级| 久久久久九九精品影院| 免费看日本二区| 国产在线精品亚洲第一网站| 精品福利观看| 一区二区日韩欧美中文字幕| 黄色女人牲交| 搡老妇女老女人老熟妇| 亚洲成人久久爱视频| 精品久久久久久久久久免费视频| 久久久久国产精品人妻aⅴ院| 99国产极品粉嫩在线观看| 亚洲欧美一区二区三区黑人| 麻豆国产av国片精品| 色av中文字幕| 午夜免费鲁丝| 婷婷亚洲欧美| 少妇 在线观看| 99re在线观看精品视频| 国产亚洲精品久久久久5区| 中文字幕精品亚洲无线码一区 | 我的亚洲天堂| 国产99白浆流出| avwww免费| 变态另类丝袜制服| 少妇裸体淫交视频免费看高清 | 露出奶头的视频| av福利片在线| 日韩大尺度精品在线看网址| 亚洲国产高清在线一区二区三 | 久久久精品欧美日韩精品| 久久香蕉激情| 夜夜躁狠狠躁天天躁| 午夜视频精品福利| 人人妻人人澡欧美一区二区| 老汉色av国产亚洲站长工具| 欧美av亚洲av综合av国产av| 超碰成人久久| 日日干狠狠操夜夜爽| 搡老妇女老女人老熟妇| 人人妻人人澡人人看| 国产aⅴ精品一区二区三区波| 99久久99久久久精品蜜桃| 黄色成人免费大全| 亚洲成人精品中文字幕电影| 亚洲精品粉嫩美女一区| 18禁裸乳无遮挡免费网站照片 | 欧美久久黑人一区二区| 韩国av一区二区三区四区| 在线观看www视频免费| 国产不卡一卡二| 亚洲av日韩精品久久久久久密| 日韩大码丰满熟妇| 国产久久久一区二区三区| 日本熟妇午夜| 12—13女人毛片做爰片一| 波多野结衣av一区二区av| 黄色女人牲交| 亚洲国产欧美一区二区综合| 国产真实乱freesex| 日本三级黄在线观看| 国产精品亚洲美女久久久| 午夜视频精品福利| 国产欧美日韩一区二区三| 亚洲aⅴ乱码一区二区在线播放 | 欧美色视频一区免费| 国产精品久久电影中文字幕| 亚洲国产欧美网| 不卡一级毛片| 色综合亚洲欧美另类图片| 19禁男女啪啪无遮挡网站| 丝袜人妻中文字幕| 久久久久久亚洲精品国产蜜桃av| 国产精品野战在线观看| 大型黄色视频在线免费观看| 淫秽高清视频在线观看| 免费在线观看亚洲国产| 国产激情欧美一区二区| 夜夜看夜夜爽夜夜摸| 免费av毛片视频| 中文字幕另类日韩欧美亚洲嫩草| 国产国语露脸激情在线看| 两个人看的免费小视频| 欧美精品亚洲一区二区| 日韩 欧美 亚洲 中文字幕| av有码第一页| 1024手机看黄色片| 日韩欧美在线二视频| 女性被躁到高潮视频| 亚洲精品在线观看二区| 欧美激情极品国产一区二区三区| 一区二区三区精品91| 在线播放国产精品三级| 亚洲色图av天堂| 一二三四社区在线视频社区8| 久久香蕉激情| 超碰成人久久| 国产精品二区激情视频| 无人区码免费观看不卡| 色老头精品视频在线观看| 国产精品免费一区二区三区在线| 日韩有码中文字幕| 亚洲在线自拍视频| 神马国产精品三级电影在线观看 | 一个人观看的视频www高清免费观看 | av视频在线观看入口| 黑人欧美特级aaaaaa片| 亚洲全国av大片| 午夜福利在线观看吧| 精品高清国产在线一区| 国产主播在线观看一区二区| 黑丝袜美女国产一区| 亚洲欧美精品综合久久99| 又黄又爽又免费观看的视频| 免费高清视频大片| 欧美人与性动交α欧美精品济南到| 国内精品久久久久精免费| 伊人久久大香线蕉亚洲五| 在线天堂中文资源库| 十分钟在线观看高清视频www| 美国免费a级毛片| 丁香欧美五月| cao死你这个sao货| 国内精品久久久久精免费| 午夜a级毛片| 亚洲五月天丁香| 欧美日韩中文字幕国产精品一区二区三区| 性欧美人与动物交配| 女同久久另类99精品国产91| 亚洲人成77777在线视频| 免费在线观看黄色视频的| 侵犯人妻中文字幕一二三四区| 黄色视频不卡| 97超级碰碰碰精品色视频在线观看| 岛国在线观看网站| 日韩三级视频一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 亚洲三区欧美一区| 久久中文字幕一级| 黑人欧美特级aaaaaa片| 丁香欧美五月| 一级作爱视频免费观看| 日本熟妇午夜| 国产精品免费一区二区三区在线| 嫁个100分男人电影在线观看| 国产精品 欧美亚洲| 一进一出抽搐gif免费好疼| av有码第一页| 十分钟在线观看高清视频www| 女人爽到高潮嗷嗷叫在线视频| 精华霜和精华液先用哪个| 日日干狠狠操夜夜爽| 少妇的丰满在线观看| 韩国av一区二区三区四区| 免费看a级黄色片| a级毛片在线看网站| 成人18禁高潮啪啪吃奶动态图| 精品国产超薄肉色丝袜足j| 天堂动漫精品| 午夜福利在线观看吧| 成人三级做爰电影| 欧美乱色亚洲激情| 久久久久久免费高清国产稀缺| 亚洲 国产 在线| 一个人观看的视频www高清免费观看 | or卡值多少钱| 久久天躁狠狠躁夜夜2o2o| www.精华液| 999精品在线视频| 一个人免费在线观看的高清视频| x7x7x7水蜜桃| 曰老女人黄片| 岛国视频午夜一区免费看| 免费在线观看亚洲国产| 日韩国内少妇激情av| 久久中文看片网| 久久久精品国产亚洲av高清涩受| 777久久人妻少妇嫩草av网站| 岛国视频午夜一区免费看| 成人午夜高清在线视频 | 亚洲国产精品999在线| 日韩av在线大香蕉| 亚洲成国产人片在线观看| ponron亚洲| 久久精品国产亚洲av高清一级| 91老司机精品| 男人舔女人下体高潮全视频| 非洲黑人性xxxx精品又粗又长| 亚洲电影在线观看av| 叶爱在线成人免费视频播放| 久久精品国产99精品国产亚洲性色| 亚洲午夜精品一区,二区,三区| 1024香蕉在线观看| 日韩一卡2卡3卡4卡2021年| 熟妇人妻久久中文字幕3abv| 欧美在线黄色| 色精品久久人妻99蜜桃| 国产成人啪精品午夜网站| 曰老女人黄片| 国产私拍福利视频在线观看| 欧美成狂野欧美在线观看| 一级毛片女人18水好多| 熟女电影av网| 91麻豆精品激情在线观看国产| 99久久国产精品久久久| 日本成人三级电影网站| 国产黄a三级三级三级人| 别揉我奶头~嗯~啊~动态视频| 老司机深夜福利视频在线观看| 亚洲,欧美精品.| 一级片免费观看大全| 国产精品一区二区免费欧美| 国产精品香港三级国产av潘金莲| 久久久精品国产亚洲av高清涩受| 可以在线观看毛片的网站| 国产又爽黄色视频| 久久性视频一级片| 亚洲中文日韩欧美视频| 日韩 欧美 亚洲 中文字幕| a在线观看视频网站| 国产熟女xx| 免费av毛片视频| 欧美日韩福利视频一区二区| 看黄色毛片网站| 99re在线观看精品视频| 视频在线观看一区二区三区| 两个人看的免费小视频| 国产爱豆传媒在线观看 | 亚洲片人在线观看| tocl精华| 香蕉久久夜色| 亚洲国产看品久久| 黄片大片在线免费观看| 国内精品久久久久精免费| 亚洲av成人不卡在线观看播放网| 99国产综合亚洲精品| 99热这里只有精品一区 | 日韩欧美一区二区三区在线观看| 俺也久久电影网| 性色av乱码一区二区三区2| 精品少妇一区二区三区视频日本电影| 国产精品日韩av在线免费观看| 两人在一起打扑克的视频| 色播在线永久视频| 欧美在线黄色| 久久九九热精品免费| 精品国产国语对白av| 国产精品免费一区二区三区在线| 两性夫妻黄色片| 欧美 亚洲 国产 日韩一| 国产一卡二卡三卡精品| 国产精品美女特级片免费视频播放器 | 亚洲真实伦在线观看| 又大又爽又粗| 亚洲一区二区三区不卡视频| 美女高潮喷水抽搐中文字幕| 精品欧美一区二区三区在线| 国产精品综合久久久久久久免费| 日韩欧美免费精品| www.熟女人妻精品国产| 欧美激情久久久久久爽电影| 岛国视频午夜一区免费看| 国产亚洲欧美98| 身体一侧抽搐| 午夜福利在线观看吧| 9191精品国产免费久久| 宅男免费午夜| 久久精品国产99精品国产亚洲性色| 悠悠久久av| 亚洲精华国产精华精| 不卡av一区二区三区| 日韩三级视频一区二区三区| 亚洲欧美精品综合久久99| 麻豆久久精品国产亚洲av| 久久久久久久精品吃奶| 中文字幕最新亚洲高清| 村上凉子中文字幕在线| 久久午夜亚洲精品久久| 热re99久久国产66热| 长腿黑丝高跟| 熟女少妇亚洲综合色aaa.| 亚洲 国产 在线| 这个男人来自地球电影免费观看| 日日夜夜操网爽| 国内揄拍国产精品人妻在线 | 国产片内射在线| 一级a爱片免费观看的视频| 国产亚洲av嫩草精品影院| 久热爱精品视频在线9| 亚洲人成77777在线视频| xxx96com| 成人三级黄色视频| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩国内少妇激情av| 久久久国产精品麻豆| 久久国产精品人妻蜜桃| 美女午夜性视频免费| 18美女黄网站色大片免费观看| 亚洲中文日韩欧美视频| 天堂动漫精品| 国产aⅴ精品一区二区三区波| 午夜免费成人在线视频| 又大又爽又粗| 欧美精品啪啪一区二区三区| 国产精品九九99| www.熟女人妻精品国产| 天天躁狠狠躁夜夜躁狠狠躁| 麻豆国产av国片精品| 亚洲午夜理论影院| 免费高清视频大片| 观看免费一级毛片| 国产真人三级小视频在线观看| 女警被强在线播放| 国产精品亚洲一级av第二区| 露出奶头的视频| 精品一区二区三区av网在线观看| 亚洲精品一区av在线观看| 日本免费a在线| 精品久久久久久久人妻蜜臀av| 黄色女人牲交| 午夜福利成人在线免费观看| 亚洲国产看品久久| 首页视频小说图片口味搜索| 久久国产精品影院| 1024视频免费在线观看| 免费无遮挡裸体视频| 91成年电影在线观看| 精品国内亚洲2022精品成人| 欧美激情 高清一区二区三区| 夜夜夜夜夜久久久久| 12—13女人毛片做爰片一| www.www免费av| 好看av亚洲va欧美ⅴa在| 一二三四社区在线视频社区8| 亚洲第一电影网av| 黄色女人牲交| 色综合欧美亚洲国产小说| 午夜福利在线在线| 别揉我奶头~嗯~啊~动态视频| 久久香蕉精品热| 久久亚洲真实| 国产精品亚洲美女久久久| 欧美性猛交黑人性爽| 脱女人内裤的视频| 免费在线观看完整版高清| 亚洲国产精品久久男人天堂| 丰满的人妻完整版| 欧美中文综合在线视频| 村上凉子中文字幕在线| 精品一区二区三区四区五区乱码| ponron亚洲| 国产精品美女特级片免费视频播放器 | 在线视频色国产色| 麻豆国产av国片精品| 妹子高潮喷水视频| 香蕉国产在线看| 亚洲精品在线美女| 亚洲精华国产精华精| 在线十欧美十亚洲十日本专区| 欧美性猛交╳xxx乱大交人| 热re99久久国产66热| 国产三级黄色录像| 欧美最黄视频在线播放免费| 88av欧美| 欧美一级a爱片免费观看看 | 色综合亚洲欧美另类图片| 午夜成年电影在线免费观看| avwww免费| 熟女少妇亚洲综合色aaa.| 国产在线观看jvid| 亚洲成a人片在线一区二区| 国产精品久久久久久亚洲av鲁大| tocl精华| 老汉色∧v一级毛片| 真人一进一出gif抽搐免费| 亚洲一区二区三区色噜噜| 国产成人影院久久av| 久久香蕉精品热| 欧美乱妇无乱码| 亚洲中文字幕一区二区三区有码在线看 | 99久久精品国产亚洲精品| 午夜免费成人在线视频| 欧美色欧美亚洲另类二区| 成人国产一区最新在线观看| 美女 人体艺术 gogo| 男人舔奶头视频| 精品一区二区三区四区五区乱码| 色在线成人网| www国产在线视频色| 国产精品久久久av美女十八| 久久精品国产99精品国产亚洲性色| 无遮挡黄片免费观看| xxx96com| 波多野结衣高清作品| 一本精品99久久精品77| 国产高清激情床上av| 久久国产精品人妻蜜桃| 国产午夜精品久久久久久| 国产精品一区二区免费欧美| 一区福利在线观看| 成人国语在线视频| 午夜福利欧美成人| 亚洲av电影不卡..在线观看| 搡老岳熟女国产| www.精华液| 久久精品国产清高在天天线| 一卡2卡三卡四卡精品乱码亚洲| 国产高清激情床上av| 久久午夜综合久久蜜桃| 男女床上黄色一级片免费看| 日韩大尺度精品在线看网址| 日本 av在线|