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

    二戰(zhàn)后美國三個大氣模式中心的美日組合

    2014-03-02 09:31:30Th聞新宇朱清照
    關(guān)鍵詞:環(huán)流大氣日本

    Th 聞新宇 朱清照

    二戰(zhàn)后美國三個大氣模式中心的美日組合

    Th 聞新宇 朱清照

    著名計算機科學(xué)家高德納(Donald E. Knuth)曾經(jīng)這樣定義科學(xué)與藝術(shù)的邊界:“科學(xué)是那些我們理解得足夠清楚,以至于可以描述成計算機語言的知識;而藝術(shù)則是其他還未達(dá)到這種程度的認(rèn)知。(Science is what we understand well enough to explain to a computer; Art is everything else we do.)”從這個意義上講,現(xiàn)代氣象學(xué)的發(fā)端似乎不是20世紀(jì)初以V. Bjerknes為首的挪威學(xué)派,而應(yīng)是二戰(zhàn)后的1950年。

    從1950年Charney首次在ENIAC上成功試驗了正壓渦度方程開始,發(fā)展數(shù)值模式就成為了那個時代大氣科學(xué)領(lǐng)域最激動人心的主旋律。Rossby于1947年從美國回到瑞典斯德哥爾摩大學(xué),領(lǐng)導(dǎo)發(fā)展了全世界第一個業(yè)務(wù)化的正壓預(yù)報模式,在瑞典BESK計算機上一周三次地預(yù)報北大西洋地區(qū)的天氣狀況。然而由于諸多原因,歐洲在模式研發(fā)的道路上雖然贏在了起跑線上,卻輸在了途中。美國最終以蓬勃向上的絕對優(yōu)勢,領(lǐng)導(dǎo)了二戰(zhàn)后數(shù)值模式的研發(fā)方向。特別地,美國在1950—1970年代,先后誕生了三個具有里程碑意義的模式研發(fā)中心:位于科羅拉多州博爾德市的美國國家大氣科學(xué)研究中心(NCAR)、位于加州洛城的加州大學(xué)洛杉磯分校(UCLA)和位于新澤西州普林斯頓市的地球流體力學(xué)實驗室(GFDL)。其中最耐人尋味的是,這三個模式中心中最早領(lǐng)導(dǎo)模式研發(fā)工作的,都是一個美國學(xué)者和一個日本學(xué)者的組合:NCAR的Warren Washington和栗原秋良(Akira Kasahara)(圖1)、UCLA的Yale Mintz和荒川昭夫(Akio Arakawa),以及GFDL的Joseph Smagorinsky和真鍋淑郎(Syukuro Manabe)。這種固定的“美日組合”究竟有哪些優(yōu)勢?這種優(yōu)勢如何推動了大氣環(huán)流模式的早期發(fā)展?這樣的組合在特殊的歷史年代又有怎樣的必然性?當(dāng)我們回望這段歷史,仔細(xì)考察這三個模式中心早期的發(fā)展經(jīng)歷時,就會更深入地理解這些看似巧合的故事其實都孕育在某種必然中。(圖2)。

    一、國家大氣科學(xué)研究中心(NCAR)

    NCAR作為美國高校大氣科學(xué)聯(lián)盟(UCAR)下屬的研究中心始建于1960年,其紅色建筑由華裔著名建筑師貝聿明設(shè)計,坐落在美麗的博爾德市西南郊Mesa路盡頭的小山上。在NCAR發(fā)展早期,領(lǐng)導(dǎo)其模式研發(fā)工作的是美國黑人氣象學(xué)家Warren Washington和日本學(xué)者笠原秋良(Akira Kasahara)。

    圖1 1975年Warren Washington和笠原秋良的合影(見于NCAR檔案)

    笠原秋良于1926年生于日本東京,1944年考入東京大學(xué)修習(xí)氣象學(xué)專業(yè),1948年畢業(yè)后師從日本著名氣象學(xué)家正野重方(Shigekata Shyono)攻讀研究生。當(dāng)時的日本處于戰(zhàn)后恢復(fù)期,美國盟軍在東京建立圖書館大肆傳播西方文化和價值觀,也向東大提供了許多美國最新的氣象學(xué)雜志,這些專業(yè)雜志使得包括笠原在內(nèi)的很多日本學(xué)者,了解到美國正在積極發(fā)展一種稱為“數(shù)值模式”的計算機與程序相結(jié)合的設(shè)備。正野重方敏銳地意識到這一革命性的新方向?qū)τ趶浐蟿恿W(xué)與經(jīng)驗天氣預(yù)報之間的鴻溝具有重要意義,并開始思考如何發(fā)展日本自己的數(shù)值模式。正野教授當(dāng)時與自己的研究助理、學(xué)生和日本氣象廳(JMA)的相關(guān)研究人員一起成立了一個數(shù)值天氣預(yù)報小組,定

    期討論相關(guān)問題。笠原曾生動地回憶道:那是1949年的一天,正野教授興奮地舉著Charney的論文沖進(jìn)教室,激動地說:“看!氣象學(xué)將因這篇文章而真正實現(xiàn)現(xiàn)代化!”1952年,正野重方的研究助理岸保勘三郎(Kanzaburo Gambo)因在長波頻散和數(shù)值預(yù)報理論方面的工作受邀去美國普林斯頓大學(xué)高級研究所與Charney合作兩年,期間他經(jīng)常寫信給笠原,告訴他Charney那邊在模式發(fā)展方面的工作進(jìn)展。1954年,岸??安┦炕氐綎|京和正野教授共同組建了更大規(guī)模的數(shù)值預(yù)報小組,共包括25名研究生和日本氣象廳雇員,他們每月都要集中討論和匯報進(jìn)展。這個活躍而充滿激情的小組,在整個1950年代一直活躍在日本數(shù)值預(yù)報領(lǐng)域的最前沿。后人也稱之為“東京NWP小組”。

    正是由于出身于正野教授的NWP小組,笠原接受了廣泛的學(xué)術(shù)熏陶和良好的專業(yè)訓(xùn)練。1956年,笠原受Platzman的邀請到芝加哥大學(xué)參與熱帶氣旋的研究,并在那里工作了6年。當(dāng)時芝加哥大學(xué)的計算設(shè)備并不好,笠原經(jīng)常要跑到位于首都華盛頓的美國氣象局去完成計算任務(wù)。1961年,美國國家自然科學(xué)基金會(NSF)和日本學(xué)術(shù)振興會(JSPS)達(dá)成協(xié)議,共同促進(jìn)美日在尖端科學(xué)領(lǐng)域的交流與合作。兩年后,笠原秋良受此計劃資助前往NCAR工作,NCAR則承諾將提供最好的研究環(huán)境和計算設(shè)備。在笠原到達(dá)NCAR工作后的第三個月,他遇到了同樣新加入NCAR的黑人青年:沃倫?華盛頓(Warren Washington)。

    沃倫?華盛頓于1936年生于美國俄勒岡州波特蘭市,1960年畢業(yè)于俄勒岡州立大學(xué)獲氣象學(xué)碩士學(xué)位,1964年畢業(yè)于賓夕法尼亞州立大學(xué)獲氣象學(xué)博士學(xué)位,成為美國歷史上第二位獲得大氣科學(xué)博士學(xué)位的非洲裔美國人。據(jù)華盛頓自己回憶,1963年8月28日,他正在驅(qū)車前往NCAR報到的路上,在汽車旅館的電視里看到了馬丁?路德?金(Martin Luther King)的著名演講“I Have a Dream”,這段演說極大地激發(fā)了他作為一位非裔美國人在科學(xué)研究領(lǐng)域中爭取平等的熱情,這對他后半生的事業(yè)都產(chǎn)生了很大的影響。

    在華盛頓加入NCAR不久后的一天,他主動向笠原提出共同建立一個大氣環(huán)流模式(GCM)的想法,這讓笠原感到非常震驚,因為那時華盛頓還只是一個剛剛畢業(yè)的青年學(xué)生。盡管笠原一直關(guān)注著GCM的發(fā)展,但此時NCAR在模式方面已經(jīng)落在了后面。GFDL、UCLA和LLNL 氣象學(xué)家Chuck Leith都已經(jīng)發(fā)展了自己的數(shù)值模式。笠原認(rèn)為1917年理查森的手算試驗不能算是完全的失敗,他所使用的控制方程在本質(zhì)上應(yīng)該可以用計算機計算并得到更精確的解。于是笠原提出,NCAR的模式研發(fā)要采用理查森的原始方程組,全球覆蓋,這可以使NCAR的大氣環(huán)流模式研發(fā)趨向科學(xué)意義上的真正完善和成熟。

    圖2 1950—1990年GFDL、NCAR和UCLA有關(guān)模式發(fā)展和模式傳播的示意圖(修改自Randall D(2000)第二章Fig 1)

    1964年5月,華盛頓和笠原秋良共同向NCAR主

    任湯普森(Thompson)提出建立一個氣候模式的想法,并得到了支持。他們最終決定購買一臺先進(jìn)的CDC3600計算機用于數(shù)值計算。在模式發(fā)展初期,他們兩人主要專注于科學(xué)方面的研究,由大量程序員負(fù)責(zé)編寫和完善代碼。整個團(tuán)隊幾乎每天都要開會討論模式的輸出結(jié)果并確定第二天的改進(jìn)方向,模式在這種不斷地修改中也變得越來越穩(wěn)定。最終在科學(xué)家和程序員的共同努力下,笠原和華盛頓共同署名于1967年發(fā)表了基于原始方程、覆蓋全球的GCM論文,在這篇文章里他們系統(tǒng)介紹了NCAR比較復(fù)雜的氣候模式的動力框架、物理過程和初步的模擬結(jié)果。作為NCAR第一代大氣環(huán)流模式,這個5?分辨率的兩層模式也被稱為Kasahara-Washington Model或NCAR 1。之后的1970年,他們進(jìn)一步大幅改進(jìn)了模式的靈活性,發(fā)布了第二代的NCAR 2,這一代模式默認(rèn)分辨率仍為水平5?,但垂直增加到6層。1973年他們又發(fā)布了NCAR 3,重點改進(jìn)了差分算法。縱觀Kasahara-Washington大氣環(huán)流模式的發(fā)展歷程,他們始終在計算方法穩(wěn)定性、極點處理,以及如何控制計算機舍入誤差方面給予了更多的關(guān)注。

    1970年代后期,NCAR逐漸放棄了Kasahara-Washington模式系列,轉(zhuǎn)而開發(fā)給美國所有大學(xué)使用的新一代氣候模式CCM(Community Climate Model),這一方面是出于更好地服務(wù)UCAR的需求,另一方面也是因為澳大利亞數(shù)值預(yù)報研究中心(ANMRC)發(fā)展了具有里程碑意義的譜模式。CCM系列一直發(fā)展至今,經(jīng)歷了CCM0B,CCM0A,CCM1,CCM2,CCM3,CSM1,CCSM2,CCSM3,CCSM4,是今天CESM1一脈的源頭。

    從1964年至今,NCAR一直努力延續(xù)著模式研發(fā)的傳統(tǒng),使得今日NCAR的全系列數(shù)值模式,包括氣候模式CESM、化學(xué)模式MOZART、天氣模式WRF等,被世界范圍內(nèi)的眾多科學(xué)家所廣泛使用。今天NCAR的諸多數(shù)值模式都已成為美國模式領(lǐng)域最重要的風(fēng)向標(biāo)。

    二、加州大學(xué)洛杉磯分校(UCLA)

    加州大學(xué)洛杉磯分校的大氣科學(xué)系是由挪威學(xué)派的J. Bjerknes(Jacob Bjerknes)于二戰(zhàn)期間建立的。1940年4月,挪威剛剛被德軍占領(lǐng),J. Bjerknes便接受美國的邀請,著手在UCLA建立氣象系為美國空軍培養(yǎng)氣象預(yù)報員。這里成為當(dāng)時美國5個基于大學(xué)的緊急軍事氣象訓(xùn)練基地之一。從那時起,UCLA氣象系逐漸成為美國最重要的大氣科學(xué)研究中心之一,在動力氣象、大氣環(huán)流模式、海氣相互作用、邊界層氣象、平流層物理化學(xué)等很多領(lǐng)域都保持了領(lǐng)先的研究水平。

    荒川昭夫(Akio Arakawa)生于1927年,中學(xué)時他就對自然科學(xué),特別是物理很感興趣。二戰(zhàn)結(jié)束后的1947年,近20歲的荒川考取了東京大學(xué)物理系,一直接受純粹的物理學(xué)訓(xùn)練。1950年荒川本科畢業(yè),也曾考慮過繼續(xù)讀研究生,但由于當(dāng)時剛剛訂婚,迫于生計需要他還是選擇了工作。與他同班畢業(yè)的有30名同學(xué),但當(dāng)時日本經(jīng)濟(jì)慘淡,根本沒有招收物理系畢業(yè)生的工作崗位,全班只有3個工作機會,其中之一便是去日本氣象廳?;拇◥貉a了氣象學(xué)知識,沒想到面試時主要被問了很多物理和數(shù)學(xué)問題,他被順利錄取了。對于自己的教育背景,荒川曾這樣回憶道:我從來沒讀過研究生,也沒有獲得過碩士、博士學(xué)位(實際上1961年他被授予東京大學(xué)名譽博士學(xué)位),我本科也沒接受過氣象學(xué)訓(xùn)練。在步入JMA工作之前,我對氣象學(xué)一無所知,也毫無興趣。

    荒川自己也沒想到,他畢業(yè)后進(jìn)入JMA后的第一份工作,卻是隨船出海。在經(jīng)過半年的培訓(xùn)后,荒川和另外幾名入職新人被送上了JMA氣象監(jiān)測船。在之后一年的時間里,他隨船往返于日本東北部X-Ray測站(152?E,39?N,測美日之間的天氣狀況)和西南部的Tango測站(135?E,29?N,測臺風(fēng))。船上的生活并不緊張,由于他們幾個新人是被作為未來JMA中層領(lǐng)導(dǎo)人員進(jìn)行培養(yǎng)的,因此荒川有很多時間觀察中緯度海面風(fēng)暴路徑上的鋒面氣旋云系和低緯度的強對流云團(tuán),而且還能經(jīng)常見習(xí)其他15名真正的氣象員每天開展日常觀測、手繪天氣圖、并進(jìn)行相關(guān)天氣討論等活動。正是這一年內(nèi)6~7次的隨船出海,使荒川真正愛上了氣象學(xué)。

    一年期的出海結(jié)束后,荒川面臨兩個選擇,或者去JMA做統(tǒng)計員,或者繼續(xù)隨船做海洋觀測。此時JMA下屬的氣象研究所(MRI)主任荒川秀?。℉idetoshi Arakawa)教授突然主動邀請荒川昭夫加入MRI從事天氣預(yù)報、分析工作。當(dāng)時MRI沒有計算機,所有預(yù)報工作都依賴在500hPa天氣圖上應(yīng)用挪威人Fjortoft(1952)的勻速平流公式。而1950年代中期,國際氣象學(xué)界最前沿的領(lǐng)域是對大氣環(huán)流的研究,比如Fultz的轉(zhuǎn)盤試驗,這深深吸引了荒川的關(guān)注。1958年荒川將國際上對大氣環(huán)流的研究進(jìn)行了總結(jié),用日語出版了單行本《大氣環(huán)流的現(xiàn)代理論》,向日本氣象學(xué)界介紹國際上有關(guān)大氣環(huán)流的最新進(jìn)展,這本書直到今天仍被很多看過的人大加稱贊。

    到了1959年,MRI幾經(jīng)權(quán)衡從美國購買了IBM

    704計算機,這比MRI原來的富士通電腦能力強很多?;拇ㄒ矎臉I(yè)務(wù)預(yù)報部門轉(zhuǎn)到了新成立的數(shù)值預(yù)報(NWP)部門。新部門沒有程序員,荒川就從零學(xué)起,對這臺計算機的硬件和剛剛誕生的Fortran語言(1957年,F(xiàn)ortran語言誕生在IBM 704機器上)逐漸有了全面的了解,并開始用計算機解決工作中的具體問題。要知道,當(dāng)時即使IBM駐日本分公司里也沒幾個人懂Fortran語言。同時,開始接觸數(shù)值計算的荒川也開始規(guī)律地參加“東京NWP小組”的每月討論,并受到了正野重方教授的賞識。1960年,國際數(shù)值預(yù)報大會(NWP Symposium)在日本東京召開(亦稱東京NWP大會)。這是一次具有里程碑意義的大會,當(dāng)時國際上著名的有關(guān)NWP和大氣環(huán)流方面的權(quán)威專家?guī)缀跸?shù)到場,包括:Bolin,Charney,Doos,Eliassen,F(xiàn)jortoft,Gates,郭曉嵐,Lorenz,Mintz,Namias,Pfeffer,Phillips,Platzman,Shuman,Smagorinsky,Wurtele等。在這次大會上,東大的正野教授與Mintz有過深入的交流,Mintz表明他想發(fā)展能模擬季節(jié)變化、可長期穩(wěn)定積分的大氣環(huán)流模式的宏大愿望,正野便向他推薦了荒川昭夫。Mintz在拜讀了一些荒川的日語論文的英譯本后,決定邀請荒川去UCLA進(jìn)行為期2年的訪問研究,其核心目的就是幫助他發(fā)展大氣環(huán)流模式。

    Yale Mintz(1916—1991年)于1937年畢業(yè)于達(dá)特茅斯學(xué)院的人類學(xué)系,1942年于哥倫比亞大學(xué)獲得地質(zhì)學(xué)的碩士學(xué)位,然后到紐約大學(xué)氣象系工作了2年,又考到美國西岸的UCLA師從創(chuàng)建氣象系的J. Bjerknes。Mintz是個熱情奔放、興趣廣泛、智力過剩的年輕人,J. Bjerknes曾生動地回憶道:他來報到的時候,除了人之外,就是滿滿一大車的書。1949年,Mintz獲得氣象學(xué)博士學(xué)位,這是該校歷史上第二個獲此學(xué)位的人(第一個是Charney)。之后,Mintz留校工作,并作為J. Bjerknes的得力助手參與很多UCLA氣象系的建設(shè)和管理工作。

    1950年代,天氣預(yù)報多是局地、大陸尺度的預(yù)報,剛剛興起的NWP預(yù)報通常也只能報2天。而受到逐漸興起的有關(guān)大氣環(huán)流的研究的啟發(fā),Mintz想發(fā)展能長期穩(wěn)定積分的GCM,并把預(yù)報推廣到全球。1961年荒川受邀來到UCLA,他們兩人便開始了一段長達(dá)20年的合作。他們兩人的工作風(fēng)格其實天差地別,Mintz由于教育背景的限制并不是一個徹底的氣象學(xué)家,但他眼界寬廣,因此他只關(guān)心宏大圖景,專注于不斷拓展GCM的功能。他總覺得荒川工作進(jìn)度很慢,而荒川總是不得不說服Mintz慢下來,并告訴他目前工作的難點在哪里?;拇ㄔ貞浀溃骸拔覀兠刻於家懻?個小時,或者是在電腦前,或者是在電話旁,他每天都有很多新想法,我通常要否定一多半,而剩下的都是非常好的點子?!庇捎诋?dāng)時UCLA比那些國家實驗室(如NCAR,GFDL,LLNL等)小很多,甚至連程序員都沒有,而Mintz又有野心要與他們競爭,所以荒川幾乎承擔(dān)了全部的代碼開發(fā)工作,而且每個周末都要加班趕工。當(dāng)時UCLA剛剛購買了一臺很先進(jìn)的IBM 709,荒川憑借自己在日本積累的編程經(jīng)驗,很小心地編寫Fortran程序卡片,實際的工作進(jìn)展并不算慢。

    1961—1963年,荒川和Mintz主要解決了平流方程中常見的雅可比算子(Arakawa’s Jacobian operators)難題,并順利發(fā)展出了第一代UCLA兩層模式,也稱Mintz-Arakawa Model。到了1963年,荒川由于學(xué)者(J)簽證到期不得不返回日本,并且2年之內(nèi)不得再入境,荒川便回到了JMA工作。離開荒川的Mintz失落不已,專門請了程序員Dennis Subsay才能繼續(xù)開展工作,同時急切地每個月給荒川寫一封信,希望他能盡快返美。1965年,荒川再次赴美,并獲得了UCLA的永久教席。他們兩人的再次合作爆發(fā)了前所未有的能量,如向模式里加入水汽、完善格點設(shè)計、改進(jìn)守恒性更好的差分算法等。

    UCLA的大氣環(huán)流模式在1965年之后發(fā)展很快,但兩人都沒有時間停下來撰寫論文,以至于有關(guān)UCLA模式的論文是三個模式中心發(fā)表最晚的;不僅如此,荒川和Mintz湊巧一生都沒有共同署名的文章發(fā)表。但是,這沒有妨礙UCLA的模式成為世界上傳播最廣泛的模式。這一方面是由于UCLA作為大學(xué)相比于那些國家實驗室更自由、更友好的開放態(tài)度所致。事實上,1965年之后,全世界各地的學(xué)生和訪問學(xué)者紛紛去UCLA交流訪問,他們畢業(yè)或回國后又把模式代碼帶到了世界各地。另一方面是由于IBM公司對UCLA提供了巨大的技術(shù)支持,IBM大規(guī)??茖W(xué)計算部的工作人員每周從加州圣何塞(San Jose)到UCLA來聽課、調(diào)程序,并最終幫助他們撰寫了高質(zhì)量的文檔,這些文檔里甚至包含了完整的源代碼,這使得世界各地的科學(xué)家可以輕易地復(fù)制UCLA的模式。

    第一代UCLA的大氣環(huán)流模式是原型模式,大體于1965年成型,沒有水汽過程,水平分辨率7?×9?,垂直兩層,物理過程方面也很不完善,荒川稱之為UCLA I。第二代模式UCLA II于1960年代后期完成,這是第一個可以稱得上是“產(chǎn)品級”的代碼,水平分辨率4 ?×5 ?,使用B格點,同樣只有兩層。UCLA II流

    傳很廣,是因為1970年Mintz的研究生Lawrence Gates畢業(yè)后將模式帶到了蘭德公司(RAND Corp.,美國智庫),并將其應(yīng)用在很多美國國防部的項目中,這個蘭德公司的版本后來又被進(jìn)一步帶到了俄勒岡州立大學(xué)。第二代模式的垂直擴展版UCLA II-3L在第二代模式的基礎(chǔ)上,主要把垂直分層增加到3層,這一版本后來又增加垂直分層后擴散到了三個NASA中心:1972年,一個9層版本傳播到NASA/GISS;幾年后又傳播到NASA/GLAS和NASA/GLA。第三代模式UCLA III重新設(shè)計了混網(wǎng)格,采用C格點,這也是后來UCLA模式一直沿用的格點設(shè)計,垂直方向增加到6或12層。1970年代中期,UCLA III被傳播到美國海軍和日本氣象廳MRI。第四代模式UCLA IV從1980年代早期開始研發(fā),重新設(shè)計了σ垂直坐標(biāo),這一版標(biāo)志著UCLA模式基本走向成熟。這一版代碼被傳播到美國海軍、NASA/GLA、DOE/LLNL,甚至中國氣象局,其中荒川的學(xué)生David Randall將這一版的模式代碼帶到科羅拉多州立大學(xué)(CSU),CSU成為后來最重要的模式發(fā)展分支之一。

    Mintz于1977年退休。他在1980年代一直在NASA/GSFC擔(dān)任科學(xué)顧問,但轉(zhuǎn)而關(guān)注陸面過程和水循環(huán)。他于1990年獲得Rossby獎,于1991年4月27日在以色列死于癌癥?;拇ㄔ缭?977年就獲得了Rossby獎,于1998年從UCLA退休,但他一直留在UCLA繼續(xù)工作至今。他的學(xué)生David Randall將荒川退休紀(jì)念會上的很多科學(xué)家的發(fā)言結(jié)集出版(General Circulation ModelDevelopment: Past, Present, and Future),這是一份全面回顧UCLA模式發(fā)展的文獻(xiàn)。

    三、地球流體力學(xué)實驗室(GFDL)

    GFDL的前身是美國氣象局于1955年在von Neumann倡議下建立的大氣環(huán)流研究部(General Circulation Research Section,GCRS),位于馬里蘭州的Suitland市,距離氣象局的聯(lián)合天氣預(yù)報小組(Joint NWP Unit)非常近,Joseph Smagorinsky受邀擔(dān)任主任。當(dāng)時這個部門的成立主要是為了與von Neumann繼續(xù)合作開發(fā)濾波模式。1959年該部門改名為大氣環(huán)流研究實驗室(General Circulation Research Laboratory,GCRL)并遷至首都華盛頓,從這時起,Smagorinsky開始獨立研發(fā)9層原始方程模式。1963年,GCRL更名為現(xiàn)在的名字GFDL,并于1968年再次遷至現(xiàn)在的普林斯頓。今天的GFDL已經(jīng)成為美國最重要的大氣科學(xué)研究中心之一,其在氣候動力和模式模擬方面都具有領(lǐng)先地位。

    真鍋淑郎(Syukuro Manabe)于1931年生于日本愛媛縣,并分別于1953,1955和1958年獲得東京大學(xué)的本科、碩士和博士學(xué)位。同他的學(xué)長笠原秋良一樣,他研究生階段的導(dǎo)師也是著名的正野重方教授。在研究生階段,真鍋就在正野教授和岸??敖淌诘闹笇?dǎo)下,在日本氣象學(xué)報上發(fā)表了幾篇論文,其聰穎、活躍的性格深得正野教授喜愛。特別地,他在“東京NWP小組”里也異常活躍,這都為他今后從事數(shù)值模式研發(fā)方面的工作打下了堅實基礎(chǔ)。受導(dǎo)師正野教授的推薦,真鍋在博士畢業(yè)的1958年就收到了來自美國Smagorinsky的邀請,畢業(yè)后直接飛赴美國加入Smagorinsky領(lǐng)導(dǎo)的GCRS。

    Joseph Smagorinsky于1924年生于美國一個油畫商之家。二戰(zhàn)時他正在紐約大學(xué)讀大二,他毅然選擇參加空軍,并先后被派到布朗大學(xué)學(xué)習(xí)數(shù)學(xué)和物理,以及MIT跟隨Lorenz學(xué)習(xí)動力氣象。在這期間,他能親身感受到戰(zhàn)爭對氣象信息的巨大渴求,也積累了很多實際的氣象觀測經(jīng)驗。于是二戰(zhàn)結(jié)束后,他決定繼續(xù)深造氣象學(xué)。在一次由Charney主講的天氣預(yù)報課程上,Smagorinsky在提問環(huán)節(jié)問了一系列令人印象深刻的問題,Charney認(rèn)可他的才華,便邀請他到普林斯頓的高級研究所去用計算機來驗證他的一系列想法。1950年4月,年輕的Smagorinsky參加了Charney領(lǐng)導(dǎo)的著名的ENIAC試驗,這使得Smagorinsky成為最早一批親身參與1950年ENIAC數(shù)值試驗的學(xué)者之一。之后他便正式加入普林斯頓高級研究所,與Charney和von Neumann合作,共同推進(jìn)數(shù)值試驗方面的工作。1953年,時年29歲的Smagorinsky得到了一個美國氣象局下屬的聯(lián)合數(shù)值天氣預(yù)報中心的工作機會。兩年后,美國氣象局又在von Neumann的建議下新成立了一個大氣環(huán)流研究部(GCRS),由Smagorinsky直接領(lǐng)導(dǎo)。如上文所述,之后,他領(lǐng)導(dǎo)的這個機構(gòu)幾度搬家,并最終以GFDL的名字遷至新澤西州的普林斯頓。

    Smagorinsky在1955—1956年,通過與von Neumann、Charney和Phillips合作,已經(jīng)發(fā)展出兩層半球模式;從1959年開始,他又與真鍋合作著手發(fā)展9層半球原始方程模式。盡管當(dāng)時他們的模式使用了過多的黏性耗散,使得模擬的結(jié)果看起來更平滑漂亮,但計算不穩(wěn)定一直是GFDL模式頭疼的問題。在GFDL,由于有大量的程序員的支持,真鍋得以全身心投入科學(xué)方面的思考中,特別是對很多物理過程如何在模式中表達(dá),他做了海量的調(diào)查研究工作。比如真鍋最早發(fā)展的對流調(diào)整方案,就是對輻射傳輸和積云對流過程簡潔的高度抽象。除了模式研發(fā)本身,真鍋還最早提出了使用模式進(jìn)行CO2加倍試驗,這成為

    今天利用模式進(jìn)行全球變暖研究的起點。Smagorinsky并不具體負(fù)責(zé)模式代碼的發(fā)展,他更多地關(guān)注方向性的問題。比如他最早意識到氣候模式不能只有大氣分量,海洋在長期氣候模擬中起著舉足輕重的作用。為此他于1961年邀請海洋模式方面的專家Kirk Bryan加入GFDL(當(dāng)時叫GCRL),這使得GFDL的模式最早于1970年代就成為了耦合的海氣模式。

    1960年代中期以后,Smagorinsky的精力更多地轉(zhuǎn)移到WMO下屬的全球大氣研究計劃(GARP)的規(guī)劃中,使得真鍋成為GFDL模式研發(fā)方面實際意義上的領(lǐng)軍人物。從這時起直到1998年真鍋從GFDL退休,他成為世界上領(lǐng)導(dǎo)模式發(fā)展時間最長、最富有進(jìn)取精神的領(lǐng)袖。在漫長的30多年的歲月中,一大批圍繞在真鍋身邊的合作者也逐漸成為模式研發(fā)領(lǐng)域的翹楚,包括:Kurihara,Strickler、Wetherald、Holloway、Stouffer和Bryan等。

    GFDL模式的命名不像其他單位那么自成一體,主要包括如下系列:MARKFORT系列,由Smagorinsky最早領(lǐng)導(dǎo)開發(fā)的9層半球原始方程模式,1960年代很多有影響力的論文都是基于這個系列;Zodiac系列,這是GFDL第二代模式,主要引入了Kurihara設(shè)計的新一代全球減網(wǎng)格,這一代模式一直使用到1970年代;Sector系列,這是一個東西方向只有60個經(jīng)度的非全球GCM,且海陸分布也高度理想化;SKYHI系列,是從1975年之后開始發(fā)展的垂直高分辨率GCM,包含對流層、平流層和中層;GFDL譜模式系列,這是GFDL在1970年代中期以后,從澳大利亞數(shù)值氣象研究中心(ANMRC)引入的,當(dāng)時ANMRC的Bourke等人在1970年代早期發(fā)展的譜模式大有席卷全球之勢;Supersource系列,1970年代后期,Holloway等人重寫GFDL的譜模式代碼,大大增加模塊化和用戶可選項,使其成為GFDL模式延續(xù)至今的一個譜模式系列,直到2000年之后,才被林山建等發(fā)展的有限體積動力框架所取代。

    Smagorinsky于2005年9月去世,享年81歲。Smagorinsky一生著作等身、獲獎無數(shù),并于1972年獲得Rossby獎。真鍋于1998年退休,之后曾回到日本擔(dān)任全球變化尖端研究院的科學(xué)顧問,但由于深感不適應(yīng)日本僵化的學(xué)術(shù)氛圍,他于2002年又返回GFDL繼續(xù)從事研究工作至今。

    四、對美日組合的反思

    1950—1970年這20余年間,美國在模式研發(fā)領(lǐng)域,除了以上所述的三大中心外,還曾誕生過三個富有個人英雄主義色彩的開拓者,比如1950年領(lǐng)導(dǎo)ENIAC試驗的Charney;1956年全世界第一個真正意義使用原始方程開展GCM試驗的Phillips;在LLNL發(fā)展出LAM模式、并用電影膠片展示其漂亮結(jié)果的Cecil Leith(外號“Chuck”Leith)。他們都曾經(jīng)對模式發(fā)展起到過重要的推動作用,但最終遺憾地沒能留下代碼遺產(chǎn)。

    本文所記述的三大模式中心NCAR、UCLA和GFDL,為什么能站在那個年代模式發(fā)展的最前沿,為什么又恰巧都是一個美國人和一個日本人的“美日組合”,卻是很值得深入探討的問題。除卻一些泛泛的原因外,我們想提出至少四個最重要的理解角度。

    第一,二戰(zhàn)后特殊的日本經(jīng)濟(jì)狀況與美日戰(zhàn)略同盟。二戰(zhàn)后日本作為戰(zhàn)敗國經(jīng)濟(jì)十分困窘,據(jù)“東京NWP小組”的回憶,當(dāng)時東京大學(xué)的老師們薪水很低(除了岸??叭晒べY較高之外),他們不得不利用課余時間去中學(xué)代課。當(dāng)時為了維持NWP小組每月一次的活動,他們還到外面去拉贊助,比如當(dāng)時的《朝日新聞》出版社就曾經(jīng)資助他們2800美元(他們自己都不知道為什么會從《朝日新聞》得到這筆慷慨的資助)。困窘的經(jīng)濟(jì)使得優(yōu)秀的日本青年不得不直面生活的挑戰(zhàn),荒川就是因此接受JMA的工作走上氣象監(jiān)測船的。當(dāng)時更多的日本優(yōu)秀青年學(xué)者都愿意赴美工作,他們大都承認(rèn)有相當(dāng)一部分原因是迫于生計。而恰巧當(dāng)時受冷戰(zhàn)影響,美日結(jié)成了特殊的同盟關(guān)系,美國不僅可以賣給JMA當(dāng)時最先進(jìn)的IBM計算機,還把國際性的NWP會議放到日本東京去開,兩國甚至達(dá)成協(xié)議促進(jìn)尖端領(lǐng)域的科技合作……這都反映出當(dāng)時的美日關(guān)系是親密無間的。正是在這樣的時代背景下,很多優(yōu)秀的日本年輕學(xué)生紛紛踏上去往美利堅的航船,在遙遠(yuǎn)的大洋彼岸繼續(xù)尋找他們的科學(xué)之路。

    第二,計算機成為戰(zhàn)后尖端科學(xué)發(fā)展的有力武器。美國在戰(zhàn)后最先發(fā)展出世界領(lǐng)先的計算機工業(yè),誕生了像IBM、Cray等先進(jìn)的計算機制造商,這為大氣環(huán)流模式的研發(fā)提供了必備的硬件平臺。當(dāng)時的瑞典、日本都有自己的計算機品牌,但由于速度不夠快,都只能做NWP的短期預(yù)報。比如當(dāng)時日本的NWP就可以完成2天的計算,卻無法更長。只有美國生產(chǎn)的高速計算機,才能真正支撐起大規(guī)模的GCM計算。這種硬件上的領(lǐng)先,為美國在GCM研發(fā)領(lǐng)域領(lǐng)先世界打出了足夠的提前量。使得只有美國才最終成為了GCM研發(fā)最好的土壤。

    第三,美日在人才培養(yǎng)方面的天然差別與互補。日本在戰(zhàn)后把教育放在了最重要的基礎(chǔ)地位,并于1947年重新頒布《教育基本法》,強調(diào)教育對完善人

    格、提升個人價值和培養(yǎng)自立精神的巨大作用,在學(xué)制、教材、教學(xué)方法上也努力與世界全面接軌。這使得日本即使在最貧窮的戰(zhàn)后恢復(fù)時期,也保證了學(xué)生能學(xué)到全世界最先進(jìn)的知識,并以日本人特有的嚴(yán)謹(jǐn)、扎實的學(xué)習(xí)習(xí)慣內(nèi)化到年輕人的頭腦中。這一切,使得當(dāng)時的日本年輕人,在個人知識層面似乎并沒有因戰(zhàn)爭而經(jīng)歷任何“斷層”。諸如荒川昭夫等人在來到美國后,可以直接上手最先進(jìn)的計算機,并用最新的Fortran語言進(jìn)行編程,再加上日本人特有的嚴(yán)謹(jǐn)、求實、腳踏實地,使得他們可以很好地完成與編程相關(guān)的工作。而美國的教育體制則更注重對人宏觀判斷力和領(lǐng)導(dǎo)力的培養(yǎng),這導(dǎo)致三個模式中心無一例外都是美國人擔(dān)任宏觀領(lǐng)導(dǎo)角色,負(fù)責(zé)設(shè)計模式發(fā)展的大方向,因為他們對總體的把控力和對未來的想象力都大大領(lǐng)先于日本人。如此一來,美日組合脫胎于各自的教育體系,盡管有著與生俱來的差別和短長,但卻具有天然的互補性。

    第四,一個人,一個小組,一次國際會議,成為歷史的關(guān)鍵節(jié)點和契機。這個人就是東京大學(xué)的正野重方教授;這個小組就是由正野重方與岸??叭晒餐I(lǐng)導(dǎo)的“東京NWP小組”;這次會議就是1960年在日本東京召開的國際NWP大會。正是由于正野教授經(jīng)常鼓勵自己的學(xué)生關(guān)注模式發(fā)展前沿、鼓勵他們大膽實踐、推薦他們到美國最頂尖的實驗室去合作,正是由于有東京NWP小組這樣的研討氛圍,才孕育出訓(xùn)練良好的笠原秋良、荒川昭夫和真鍋淑郎這一批模式專家。不僅這三人,還有很多同時代的日本氣象學(xué)家皆出自這一血脈,比如Kurihara和Matsuno等。正是由于有1960年“巨星”云集的東京NWP大會這樣的平臺,使得東京NWP小組得以向世界展示其研究水平和人才儲備,并因此走出了像荒川昭夫等許多優(yōu)秀的年輕人。因此從某種意義上說,美國三大模式中心的美日組合現(xiàn)象,應(yīng)當(dāng)歸功于東京大學(xué)正野重方教授的遠(yuǎn)見和推動,應(yīng)當(dāng)歸功于一個人、一個小組、一次國際會議這樣的關(guān)鍵性樞紐。

    20世紀(jì)中期,模式發(fā)展的浪潮將日漸成熟的動力學(xué)應(yīng)用在傳統(tǒng)的天氣學(xué)分析之上,徹底改變了天氣預(yù)報和氣候預(yù)測的面貌,使得大氣科學(xué)最終站上了現(xiàn)代科學(xué)的舞臺。一個時代漸漸遠(yuǎn)去,但那些五彩斑斕,甚至有些激動人心的故事并沒有因此被人們忘記。細(xì)細(xì)品味他們背后的種種偶然,感慨之余不免發(fā)現(xiàn)其實一切又在必然中。

    致謝:感謝中科院戰(zhàn)略先導(dǎo)專項(XDA05080801)和國家自然科學(xué)基金(41005035, 41130962,41130105)對本文的支持。并特別感謝日本早稻田大學(xué)水上弘子女士對本文日文資料的翻譯和整理。

    (作者單位:北京大學(xué)物理學(xué)院大氣與海洋科學(xué)系)

    深入閱讀

    Anon. “Suki” Manabe, Pioneer of Climate modeling. IPRC Climate, 2005, 5(2): 11-15. http://iprc.soest.hawaii.edu/ newsletters/iprc_climate_vol5_no2.pdf

    Johnson D R. 1996. On the scientif i c contributions and insight of Professor Yale Mintz. J Climate, 9: 3211-3224.

    Lewis J M. 1993. Meteorologists from the University of Tokyo: Their exodus to the United States following World War II. Bull Amer Meteor Soc, 74(7): 1351-1360.

    Randall D. 2000. General Circulation Model Development: Past, Present, and Future. Utah: Academic Press.

    Smagorinsky J. 1983. The beginnings of numerical weather prediction and general circulation modeling: Early recollections. Advances in Geophysics, 25: 36.

    猜你喜歡
    環(huán)流大氣日本
    大氣的呵護(hù)
    軍事文摘(2023年10期)2023-06-09 09:15:06
    日本元旦是新年
    華人時刊(2022年3期)2022-04-26 14:29:08
    內(nèi)環(huán)流控溫技術(shù)應(yīng)用實踐與發(fā)展前景
    探尋日本
    中華手工(2021年2期)2021-09-15 02:21:08
    《黃金時代》日本版
    電影(2019年3期)2019-04-04 11:57:16
    熱鹽環(huán)流方程全局弱解的存在性
    謎底大揭秘
    大氣古樸揮灑自如
    大氣、水之后,土十條來了
    去日本怎么玩?它告訴你
    寂寞人妻少妇视频99o| 日韩 亚洲 欧美在线| 日日摸夜夜添夜夜爱| 久久久久人妻精品一区果冻| 久久国产亚洲av麻豆专区| 我的女老师完整版在线观看| 日韩av不卡免费在线播放| 国产午夜精品一二区理论片| 国产欧美日韩一区二区三区在线| 内地一区二区视频在线| 免费日韩欧美在线观看| 天堂8中文在线网| 日韩不卡一区二区三区视频在线| 大香蕉久久网| 狠狠精品人妻久久久久久综合| 亚洲欧美中文字幕日韩二区| 久久久精品94久久精品| 色婷婷av一区二区三区视频| 国产极品天堂在线| 国产av码专区亚洲av| 亚洲精华国产精华液的使用体验| tube8黄色片| 99热全是精品| 国产精品嫩草影院av在线观看| 毛片一级片免费看久久久久| 久久午夜综合久久蜜桃| 亚洲av.av天堂| 久久久久精品久久久久真实原创| 蜜臀久久99精品久久宅男| 国产精品久久久久久久电影| 女性生殖器流出的白浆| 秋霞伦理黄片| 精品亚洲成国产av| 国产一级毛片在线| 亚洲人与动物交配视频| 午夜福利影视在线免费观看| 两个人看的免费小视频| 高清av免费在线| 精品国产一区二区久久| av免费观看日本| 制服丝袜香蕉在线| 免费av不卡在线播放| 中国美白少妇内射xxxbb| 久久精品aⅴ一区二区三区四区 | 国产成人免费观看mmmm| 汤姆久久久久久久影院中文字幕| 女人久久www免费人成看片| 久久精品aⅴ一区二区三区四区 | 高清不卡的av网站| 免费av中文字幕在线| 亚洲av.av天堂| 日本wwww免费看| 国产一区二区三区av在线| 免费观看在线日韩| 99国产综合亚洲精品| 捣出白浆h1v1| 蜜桃在线观看..| 久久久久久久亚洲中文字幕| 美女视频免费永久观看网站| 一级黄片播放器| 在线免费观看不下载黄p国产| 午夜福利影视在线免费观看| 免费观看性生交大片5| 精品一区在线观看国产| 久久这里只有精品19| 免费观看a级毛片全部| av在线观看视频网站免费| 日韩av不卡免费在线播放| 你懂的网址亚洲精品在线观看| 精品人妻熟女毛片av久久网站| 国产午夜精品一二区理论片| 成年人免费黄色播放视频| 欧美激情极品国产一区二区三区 | 日本-黄色视频高清免费观看| 久久久久久久久久久免费av| 亚洲在久久综合| 国产伦理片在线播放av一区| 久久国产精品男人的天堂亚洲 | 看非洲黑人一级黄片| freevideosex欧美| 国产黄色免费在线视频| 亚洲一码二码三码区别大吗| 肉色欧美久久久久久久蜜桃| 香蕉丝袜av| 久久女婷五月综合色啪小说| 欧美最新免费一区二区三区| 久久久久久久久久久免费av| 国产免费一级a男人的天堂| 国产精品无大码| 一级毛片电影观看| av不卡在线播放| 纯流量卡能插随身wifi吗| 午夜影院在线不卡| av有码第一页| av天堂久久9| 国产成人a∨麻豆精品| av卡一久久| 久久97久久精品| 久久午夜福利片| 国产精品熟女久久久久浪| 国产精品一区www在线观看| 亚洲欧美日韩另类电影网站| 国产1区2区3区精品| 午夜老司机福利剧场| 国产亚洲精品久久久com| 婷婷色综合www| 亚洲五月色婷婷综合| 婷婷色av中文字幕| 久久人人97超碰香蕉20202| 久久精品国产a三级三级三级| av在线app专区| videos熟女内射| 精品视频人人做人人爽| 亚洲精品一二三| 90打野战视频偷拍视频| 亚洲成色77777| 欧美xxⅹ黑人| 亚洲av国产av综合av卡| 如日韩欧美国产精品一区二区三区| 黑人高潮一二区| av国产精品久久久久影院| 国产精品一二三区在线看| 宅男免费午夜| 日韩av在线免费看完整版不卡| 一级毛片 在线播放| www.熟女人妻精品国产 | 精品卡一卡二卡四卡免费| 国产在视频线精品| 黄色怎么调成土黄色| 久久ye,这里只有精品| 香蕉国产在线看| 日韩伦理黄色片| 亚洲精品乱久久久久久| 色哟哟·www| 亚洲av电影在线观看一区二区三区| 人人妻人人添人人爽欧美一区卜| 欧美日韩视频高清一区二区三区二| 久久97久久精品| 精品久久久精品久久久| av女优亚洲男人天堂| 超碰97精品在线观看| 国产一区有黄有色的免费视频| 精品视频人人做人人爽| 日韩成人伦理影院| 久久鲁丝午夜福利片| 肉色欧美久久久久久久蜜桃| 亚洲一区二区三区欧美精品| 欧美bdsm另类| 日韩中文字幕视频在线看片| 亚洲内射少妇av| 久久人人爽av亚洲精品天堂| 免费观看性生交大片5| 亚洲精品av麻豆狂野| av国产久精品久网站免费入址| 久久国内精品自在自线图片| 国产成人aa在线观看| 女的被弄到高潮叫床怎么办| av在线app专区| 久久久精品免费免费高清| 高清视频免费观看一区二区| 中国美白少妇内射xxxbb| 亚洲国产日韩一区二区| 香蕉丝袜av| 日韩av不卡免费在线播放| 波多野结衣一区麻豆| 最新中文字幕久久久久| av免费在线看不卡| 一级片'在线观看视频| 亚洲天堂av无毛| 18在线观看网站| 极品人妻少妇av视频| 一级毛片电影观看| av在线播放精品| av不卡在线播放| a级毛色黄片| 亚洲av在线观看美女高潮| 日韩视频在线欧美| 国产成人一区二区在线| 成人免费观看视频高清| 人成视频在线观看免费观看| 日本午夜av视频| 日韩大片免费观看网站| 97在线视频观看| 亚洲精品久久午夜乱码| 美国免费a级毛片| 激情视频va一区二区三区| av黄色大香蕉| 免费看不卡的av| 亚洲中文av在线| 欧美日韩一区二区视频在线观看视频在线| 中文天堂在线官网| 观看av在线不卡| 国产成人精品婷婷| 亚洲情色 制服丝袜| 国产成人免费无遮挡视频| 黄色视频在线播放观看不卡| 伊人久久国产一区二区| 欧美97在线视频| 99国产精品免费福利视频| 人体艺术视频欧美日本| 国产国语露脸激情在线看| 成人国语在线视频| 成人18禁高潮啪啪吃奶动态图| 久久国产亚洲av麻豆专区| 久久国内精品自在自线图片| 麻豆精品久久久久久蜜桃| 伊人久久国产一区二区| 在线观看一区二区三区激情| 国产一区有黄有色的免费视频| 久久av网站| 高清黄色对白视频在线免费看| 超色免费av| 精品久久久久久电影网| 男女国产视频网站| 王馨瑶露胸无遮挡在线观看| 看免费av毛片| 国产精品一二三区在线看| 久久国产精品男人的天堂亚洲 | 一本色道久久久久久精品综合| 亚洲国产精品一区三区| 午夜91福利影院| 亚洲五月色婷婷综合| 在线观看人妻少妇| 十八禁网站网址无遮挡| 精品福利永久在线观看| 菩萨蛮人人尽说江南好唐韦庄| 精品第一国产精品| 男人爽女人下面视频在线观看| 亚洲伊人色综图| 日韩成人伦理影院| 成人毛片a级毛片在线播放| 人人妻人人澡人人看| 国产成人av激情在线播放| 欧美xxxx性猛交bbbb| 老女人水多毛片| 亚洲中文av在线| 免费人妻精品一区二区三区视频| 免费不卡的大黄色大毛片视频在线观看| 9色porny在线观看| 中文乱码字字幕精品一区二区三区| 免费黄网站久久成人精品| 97超碰精品成人国产| 日本免费在线观看一区| 中文字幕人妻熟女乱码| 国国产精品蜜臀av免费| 好男人视频免费观看在线| 黄色 视频免费看| 少妇被粗大的猛进出69影院 | 日本午夜av视频| 亚洲内射少妇av| 咕卡用的链子| 性色av一级| 久热这里只有精品99| 成人亚洲欧美一区二区av| 欧美日韩精品成人综合77777| 一边亲一边摸免费视频| 亚洲婷婷狠狠爱综合网| 激情视频va一区二区三区| 久久狼人影院| 中文字幕人妻丝袜制服| 观看av在线不卡| av一本久久久久| 五月天丁香电影| 一个人免费看片子| av有码第一页| 欧美精品一区二区免费开放| 欧美国产精品一级二级三级| 我要看黄色一级片免费的| 超碰97精品在线观看| 久久av网站| 欧美亚洲 丝袜 人妻 在线| 男女边吃奶边做爰视频| 欧美xxxx性猛交bbbb| 久久婷婷青草| 国产福利在线免费观看视频| 成人国语在线视频| 国产精品久久久久久av不卡| 欧美 日韩 精品 国产| 蜜桃在线观看..| 高清黄色对白视频在线免费看| 爱豆传媒免费全集在线观看| av不卡在线播放| 国产乱人偷精品视频| 在线天堂中文资源库| 成人手机av| 免费av中文字幕在线| av女优亚洲男人天堂| 国产亚洲av片在线观看秒播厂| 国产精品.久久久| 中文字幕人妻熟女乱码| 久久综合国产亚洲精品| 最近2019中文字幕mv第一页| 美女中出高潮动态图| 精品熟女少妇av免费看| 日本wwww免费看| 日日摸夜夜添夜夜爱| 免费看av在线观看网站| 日韩在线高清观看一区二区三区| 高清黄色对白视频在线免费看| 美女国产视频在线观看| 午夜福利视频精品| 欧美日韩亚洲高清精品| 亚洲,欧美,日韩| 欧美亚洲 丝袜 人妻 在线| 国产高清国产精品国产三级| 9热在线视频观看99| 久久精品aⅴ一区二区三区四区 | 国产欧美日韩一区二区三区在线| 伦精品一区二区三区| av在线观看视频网站免费| 男女边吃奶边做爰视频| 熟女人妻精品中文字幕| 久热这里只有精品99| 欧美精品一区二区大全| 有码 亚洲区| 国产老妇伦熟女老妇高清| 国产日韩欧美在线精品| 全区人妻精品视频| 91国产中文字幕| 男女边吃奶边做爰视频| 国产 精品1| 成年女人在线观看亚洲视频| 天天操日日干夜夜撸| 国产男人的电影天堂91| 热99国产精品久久久久久7| 久久国产精品男人的天堂亚洲 | 久久 成人 亚洲| 女人被躁到高潮嗷嗷叫费观| 国产欧美亚洲国产| 日韩中字成人| 高清不卡的av网站| 三上悠亚av全集在线观看| 日韩,欧美,国产一区二区三区| 国产一区二区三区av在线| 日本av手机在线免费观看| 国产欧美日韩一区二区三区在线| 黄色一级大片看看| 婷婷色综合大香蕉| 亚洲高清免费不卡视频| 高清毛片免费看| 看十八女毛片水多多多| 免费在线观看黄色视频的| 久久精品久久精品一区二区三区| 国产精品一区www在线观看| 久久久久久久亚洲中文字幕| 永久网站在线| 97在线视频观看| 国产亚洲av片在线观看秒播厂| 久久99精品国语久久久| 国产1区2区3区精品| 男人添女人高潮全过程视频| 香蕉国产在线看| 好男人视频免费观看在线| av黄色大香蕉| 精品国产露脸久久av麻豆| 18禁裸乳无遮挡动漫免费视频| www日本在线高清视频| 免费人成在线观看视频色| 美女大奶头黄色视频| 男女国产视频网站| 国产精品国产三级国产专区5o| 欧美老熟妇乱子伦牲交| 日韩一本色道免费dvd| 国产亚洲午夜精品一区二区久久| 狂野欧美激情性bbbbbb| 成年人免费黄色播放视频| 久久这里只有精品19| 黄色毛片三级朝国网站| 极品少妇高潮喷水抽搐| 久久99一区二区三区| 国产精品一区二区在线不卡| 狠狠精品人妻久久久久久综合| 国产精品一区二区在线观看99| 少妇人妻久久综合中文| 我要看黄色一级片免费的| 国产av一区二区精品久久| 免费高清在线观看视频在线观看| 精品第一国产精品| 成年av动漫网址| 午夜福利网站1000一区二区三区| 国产高清国产精品国产三级| 欧美日韩亚洲高清精品| 欧美最新免费一区二区三区| 亚洲av.av天堂| 美女国产高潮福利片在线看| 中国美白少妇内射xxxbb| 日韩伦理黄色片| a级毛片黄视频| 伦精品一区二区三区| 亚洲欧美成人综合另类久久久| 性色avwww在线观看| 精品一品国产午夜福利视频| 亚洲精品久久久久久婷婷小说| 啦啦啦中文免费视频观看日本| 国产在线一区二区三区精| 纵有疾风起免费观看全集完整版| 狂野欧美激情性xxxx在线观看| 国产成人欧美| 国产有黄有色有爽视频| 精品一区在线观看国产| 国产精品久久久久久久久免| 久久久久久久久久人人人人人人| 日韩三级伦理在线观看| 狂野欧美激情性bbbbbb| 久久精品国产亚洲av涩爱| 国产成人91sexporn| 少妇人妻精品综合一区二区| 久久精品国产亚洲av涩爱| 搡女人真爽免费视频火全软件| 我要看黄色一级片免费的| 秋霞伦理黄片| 五月伊人婷婷丁香| 久久精品国产综合久久久 | 五月伊人婷婷丁香| 国产在视频线精品| 国产又色又爽无遮挡免| 国产极品粉嫩免费观看在线| 曰老女人黄片| 欧美激情极品国产一区二区三区 | 国产日韩欧美亚洲二区| 日韩av不卡免费在线播放| 日日爽夜夜爽网站| 久久国产精品男人的天堂亚洲 | 在线免费观看不下载黄p国产| 精品酒店卫生间| 九色成人免费人妻av| 欧美成人精品欧美一级黄| 欧美精品一区二区大全| 日日摸夜夜添夜夜爱| 国产淫语在线视频| 中文字幕人妻熟女乱码| 老司机影院成人| 捣出白浆h1v1| 一级毛片电影观看| 男女午夜视频在线观看 | a 毛片基地| 欧美激情极品国产一区二区三区 | 国产色爽女视频免费观看| 两性夫妻黄色片 | 国产亚洲av片在线观看秒播厂| 国产精品国产三级专区第一集| 毛片一级片免费看久久久久| 午夜免费男女啪啪视频观看| 日韩不卡一区二区三区视频在线| 精品一区二区三区视频在线| 九草在线视频观看| 大片电影免费在线观看免费| 女人被躁到高潮嗷嗷叫费观| 亚洲av欧美aⅴ国产| 久久久久久久精品精品| 国产免费又黄又爽又色| 欧美人与性动交α欧美精品济南到 | 老司机影院成人| 久久久久精品人妻al黑| 一级毛片 在线播放| 国产无遮挡羞羞视频在线观看| av电影中文网址| 成人午夜精彩视频在线观看| 日本vs欧美在线观看视频| 一二三四在线观看免费中文在 | 国产精品99久久99久久久不卡 | 国产成人av激情在线播放| 国产精品蜜桃在线观看| 久久久久久久久久久久大奶| 亚洲色图综合在线观看| 久久精品aⅴ一区二区三区四区 | 国产成人精品福利久久| 国产亚洲午夜精品一区二区久久| 亚洲国产成人一精品久久久| 精品福利永久在线观看| av卡一久久| 色5月婷婷丁香| 亚洲成国产人片在线观看| 国产精品蜜桃在线观看| 亚洲精品日韩在线中文字幕| 午夜福利视频精品| 欧美老熟妇乱子伦牲交| 大码成人一级视频| 狠狠精品人妻久久久久久综合| 久久久久久久久久久久大奶| 午夜91福利影院| 久久久久精品人妻al黑| 欧美另类一区| 亚洲欧美色中文字幕在线| 看非洲黑人一级黄片| 一二三四在线观看免费中文在 | 日韩伦理黄色片| 欧美成人午夜精品| 国产免费一级a男人的天堂| 丰满少妇做爰视频| 不卡视频在线观看欧美| 男女国产视频网站| 又黄又爽又刺激的免费视频.| 亚洲,欧美,日韩| 免费av中文字幕在线| 亚洲欧美色中文字幕在线| av片东京热男人的天堂| 午夜精品国产一区二区电影| 超色免费av| 人妻少妇偷人精品九色| 久久国产精品大桥未久av| 色视频在线一区二区三区| 国产乱人偷精品视频| 观看av在线不卡| 一边摸一边做爽爽视频免费| 日本-黄色视频高清免费观看| av卡一久久| 国产一区二区激情短视频 | 日本黄色日本黄色录像| 搡老乐熟女国产| 成人综合一区亚洲| 亚洲综合色惰| 国产深夜福利视频在线观看| 日韩中字成人| 另类精品久久| 亚洲一区二区三区欧美精品| 在线 av 中文字幕| 欧美最新免费一区二区三区| 观看美女的网站| 国产不卡av网站在线观看| 亚洲精品色激情综合| 高清毛片免费看| 免费观看a级毛片全部| 国产精品国产三级国产专区5o| 母亲3免费完整高清在线观看 | 七月丁香在线播放| 纯流量卡能插随身wifi吗| 国产精品无大码| 久久精品夜色国产| 国产女主播在线喷水免费视频网站| 99九九在线精品视频| 国产一区二区激情短视频 | 亚洲国产精品999| 大码成人一级视频| videos熟女内射| 九九爱精品视频在线观看| 国产黄频视频在线观看| 最近2019中文字幕mv第一页| 欧美日韩一区二区视频在线观看视频在线| 成人18禁高潮啪啪吃奶动态图| h视频一区二区三区| 黑人高潮一二区| 十八禁高潮呻吟视频| 26uuu在线亚洲综合色| 妹子高潮喷水视频| 亚洲av欧美aⅴ国产| 女性被躁到高潮视频| 天堂中文最新版在线下载| 性色avwww在线观看| 亚洲经典国产精华液单| 一区二区三区四区激情视频| 久久av网站| 美女xxoo啪啪120秒动态图| 久久毛片免费看一区二区三区| 日本猛色少妇xxxxx猛交久久| 日韩视频在线欧美| 欧美精品国产亚洲| 欧美精品高潮呻吟av久久| av片东京热男人的天堂| 伊人久久国产一区二区| 成人国产麻豆网| 18禁观看日本| 草草在线视频免费看| 日韩制服丝袜自拍偷拍| 成年美女黄网站色视频大全免费| 美国免费a级毛片| 九九在线视频观看精品| 日韩欧美精品免费久久| 人体艺术视频欧美日本| 午夜老司机福利剧场| 精品国产一区二区三区久久久樱花| 美女内射精品一级片tv| 欧美最新免费一区二区三区| 久久精品国产综合久久久 | 欧美国产精品一级二级三级| 亚洲三级黄色毛片| 久久人人爽人人爽人人片va| av在线老鸭窝| 少妇的丰满在线观看| 国产日韩欧美在线精品| av又黄又爽大尺度在线免费看| 国产xxxxx性猛交| 高清av免费在线| 久久精品aⅴ一区二区三区四区 | 狠狠婷婷综合久久久久久88av| 美女福利国产在线| 日韩熟女老妇一区二区性免费视频| 国产1区2区3区精品| 女人久久www免费人成看片| 久久精品久久精品一区二区三区| 伦精品一区二区三区| 午夜免费鲁丝| 精品人妻熟女毛片av久久网站| 亚洲情色 制服丝袜| 欧美老熟妇乱子伦牲交| 亚洲av电影在线观看一区二区三区| 久久久欧美国产精品| 成人18禁高潮啪啪吃奶动态图| 国产毛片在线视频| 人人澡人人妻人| 十八禁网站网址无遮挡| 日日爽夜夜爽网站| 成人亚洲精品一区在线观看| 日韩欧美精品免费久久| 欧美日韩一区二区视频在线观看视频在线| 欧美国产精品va在线观看不卡| 高清黄色对白视频在线免费看| 五月天丁香电影| 日韩人妻精品一区2区三区| 少妇 在线观看| 22中文网久久字幕| 水蜜桃什么品种好| 国产精品成人在线| 18禁动态无遮挡网站| 色婷婷av一区二区三区视频|