• <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)流方程全局弱解的存在性
    謎底大揭秘
    大氣古樸揮灑自如
    大氣、水之后,土十條來了
    去日本怎么玩?它告訴你
    xxxwww97欧美| 美女cb高潮喷水在线观看| 成熟少妇高潮喷水视频| 少妇的逼水好多| 内地一区二区视频在线| 春色校园在线视频观看| 草草在线视频免费看| 亚洲人与动物交配视频| 最好的美女福利视频网| 国产探花极品一区二区| 亚洲人与动物交配视频| 亚洲欧美日韩高清专用| 色av中文字幕| 在线观看舔阴道视频| 99热这里只有是精品50| 国产探花极品一区二区| 国产一区二区在线观看日韩| 午夜亚洲福利在线播放| 午夜精品在线福利| 成人美女网站在线观看视频| 国产午夜福利久久久久久| 此物有八面人人有两片| 在线天堂最新版资源| 国产精品福利在线免费观看| 日本一本二区三区精品| 变态另类丝袜制服| 精品人妻1区二区| 午夜久久久久精精品| 高清日韩中文字幕在线| 日日夜夜操网爽| 国产高清视频在线观看网站| 免费在线观看影片大全网站| 一级a爱片免费观看的视频| 亚洲国产欧美人成| 久久午夜亚洲精品久久| 日韩大尺度精品在线看网址| 变态另类成人亚洲欧美熟女| 精品国产三级普通话版| 国产高清有码在线观看视频| 国产人妻一区二区三区在| 99热这里只有是精品在线观看| 欧美国产日韩亚洲一区| 日本成人三级电影网站| 国产精品人妻久久久影院| 听说在线观看完整版免费高清| av女优亚洲男人天堂| 国产 一区 欧美 日韩| 老熟妇乱子伦视频在线观看| 国产精品av视频在线免费观看| 国产一区二区三区在线臀色熟女| 日本撒尿小便嘘嘘汇集6| 无遮挡黄片免费观看| 国产高清视频在线观看网站| 无人区码免费观看不卡| 搞女人的毛片| 18禁黄网站禁片免费观看直播| 亚洲av成人av| 99久久精品国产国产毛片| 亚洲一区高清亚洲精品| 一个人看视频在线观看www免费| 国产精品综合久久久久久久免费| 久久久久久国产a免费观看| 国内精品宾馆在线| 中文字幕久久专区| 1000部很黄的大片| 变态另类丝袜制服| 夜夜夜夜夜久久久久| 啦啦啦韩国在线观看视频| 人妻夜夜爽99麻豆av| 亚洲av美国av| 精品一区二区三区视频在线| av在线老鸭窝| 国产亚洲精品久久久久久毛片| 深夜a级毛片| 噜噜噜噜噜久久久久久91| 1024手机看黄色片| 国产激情偷乱视频一区二区| .国产精品久久| av中文乱码字幕在线| 女生性感内裤真人,穿戴方法视频| 日韩亚洲欧美综合| 少妇的逼水好多| 99国产极品粉嫩在线观看| 精品久久久久久久久av| 在线看三级毛片| 免费人成在线观看视频色| 国内少妇人妻偷人精品xxx网站| 色哟哟·www| 能在线免费观看的黄片| 国产精品伦人一区二区| 内地一区二区视频在线| 日韩强制内射视频| 欧美色视频一区免费| 91久久精品国产一区二区成人| 亚洲五月天丁香| 国产真实乱freesex| 国产黄色小视频在线观看| 少妇人妻一区二区三区视频| 老熟妇仑乱视频hdxx| 麻豆久久精品国产亚洲av| 韩国av一区二区三区四区| 国产探花在线观看一区二区| 99久久成人亚洲精品观看| 国产在线精品亚洲第一网站| 男女那种视频在线观看| 亚洲一区二区三区色噜噜| 久久久久久伊人网av| 九九热线精品视视频播放| 亚洲精品一区av在线观看| 听说在线观看完整版免费高清| 国产伦精品一区二区三区视频9| 久久国产乱子免费精品| 亚洲国产欧美人成| 国产精品女同一区二区软件 | 欧美最新免费一区二区三区| 黄色欧美视频在线观看| 淫妇啪啪啪对白视频| 国产午夜精品论理片| 亚洲专区国产一区二区| 欧美三级亚洲精品| 99久久九九国产精品国产免费| 午夜日韩欧美国产| 亚洲在线自拍视频| 夜夜夜夜夜久久久久| 国产亚洲av嫩草精品影院| 深夜a级毛片| 国产三级在线视频| 午夜免费成人在线视频| 久久精品人妻少妇| 国产精品久久久久久久电影| 动漫黄色视频在线观看| avwww免费| 国产成人福利小说| 成人三级黄色视频| 精品久久久久久久久亚洲 | 偷拍熟女少妇极品色| 嫩草影视91久久| 日韩欧美精品v在线| 日韩精品有码人妻一区| 五月玫瑰六月丁香| 久久精品国产亚洲网站| 久久人人精品亚洲av| 国产精品精品国产色婷婷| 桃色一区二区三区在线观看| 免费无遮挡裸体视频| 91麻豆精品激情在线观看国产| 欧美色视频一区免费| 欧美高清成人免费视频www| 色在线成人网| 精品久久久久久久末码| 成人精品一区二区免费| 搡女人真爽免费视频火全软件 | 成人无遮挡网站| 观看免费一级毛片| 日日啪夜夜撸| 国模一区二区三区四区视频| 看黄色毛片网站| 精华霜和精华液先用哪个| 国产精品久久久久久av不卡| 国产精品久久视频播放| 国产精品一区二区免费欧美| 国产真实伦视频高清在线观看 | 亚洲成人久久性| 琪琪午夜伦伦电影理论片6080| 热99re8久久精品国产| 国产精品亚洲美女久久久| 99久久精品热视频| 精品无人区乱码1区二区| 日韩国内少妇激情av| 国产高清激情床上av| 国产亚洲av嫩草精品影院| 日韩在线高清观看一区二区三区 | 久99久视频精品免费| 中文资源天堂在线| 免费电影在线观看免费观看| 国产亚洲欧美98| 三级男女做爰猛烈吃奶摸视频| 欧美黑人欧美精品刺激| 老司机福利观看| 999久久久精品免费观看国产| 天堂√8在线中文| 国产亚洲欧美98| 精品不卡国产一区二区三区| 淫秽高清视频在线观看| 国产 一区 欧美 日韩| 神马国产精品三级电影在线观看| 超碰av人人做人人爽久久| 99在线视频只有这里精品首页| 男人狂女人下面高潮的视频| 免费黄网站久久成人精品| 成人精品一区二区免费| 亚洲av中文av极速乱 | 亚洲欧美精品综合久久99| 国产伦在线观看视频一区| 不卡一级毛片| 少妇被粗大猛烈的视频| 琪琪午夜伦伦电影理论片6080| 麻豆国产97在线/欧美| 成年女人永久免费观看视频| 久久天躁狠狠躁夜夜2o2o| 婷婷丁香在线五月| 亚洲精品一卡2卡三卡4卡5卡| 五月玫瑰六月丁香| 日本在线视频免费播放| 一区二区三区高清视频在线| 国产一区二区激情短视频| 能在线免费观看的黄片| 午夜福利成人在线免费观看| 免费人成在线观看视频色| 国语自产精品视频在线第100页| 国产乱人伦免费视频| 不卡一级毛片| 久久久久久久久大av| АⅤ资源中文在线天堂| 欧美成人一区二区免费高清观看| 久久午夜亚洲精品久久| 此物有八面人人有两片| 香蕉av资源在线| 午夜老司机福利剧场| 成人av在线播放网站| 亚洲不卡免费看| 伊人久久精品亚洲午夜| 无遮挡黄片免费观看| 久久精品国产自在天天线| 国产精品国产三级国产av玫瑰| 女生性感内裤真人,穿戴方法视频| 超碰av人人做人人爽久久| 亚洲人与动物交配视频| 欧美潮喷喷水| 久久精品夜夜夜夜夜久久蜜豆| 俺也久久电影网| 在线观看舔阴道视频| 久久精品国产清高在天天线| 久久草成人影院| 国产熟女欧美一区二区| 午夜影院日韩av| 久久久久久久久中文| 中文字幕人妻熟人妻熟丝袜美| 亚洲av中文字字幕乱码综合| 又黄又爽又免费观看的视频| 国产精品久久久久久亚洲av鲁大| 美女高潮喷水抽搐中文字幕| 嫩草影院入口| 免费在线观看成人毛片| 国产大屁股一区二区在线视频| 亚洲av五月六月丁香网| 中文字幕熟女人妻在线| 欧美潮喷喷水| 国产91精品成人一区二区三区| 99热这里只有是精品在线观看| 成年女人看的毛片在线观看| 亚洲av第一区精品v没综合| 精品久久久噜噜| 自拍偷自拍亚洲精品老妇| 亚洲四区av| 不卡视频在线观看欧美| 亚洲最大成人av| 十八禁国产超污无遮挡网站| 一进一出抽搐gif免费好疼| 亚洲一级一片aⅴ在线观看| 亚洲aⅴ乱码一区二区在线播放| 久久久午夜欧美精品| 乱码一卡2卡4卡精品| 黄色配什么色好看| 欧美日韩精品成人综合77777| 国内揄拍国产精品人妻在线| 69人妻影院| 中文在线观看免费www的网站| 欧美日韩黄片免| 国产精品三级大全| 天美传媒精品一区二区| 床上黄色一级片| 桃红色精品国产亚洲av| 老师上课跳d突然被开到最大视频| 熟女电影av网| 亚洲av中文av极速乱 | 精品欧美国产一区二区三| 成人一区二区视频在线观看| 成人精品一区二区免费| 五月伊人婷婷丁香| 床上黄色一级片| 国产精品一区www在线观看 | 久久久久免费精品人妻一区二区| 久久天躁狠狠躁夜夜2o2o| 老司机午夜福利在线观看视频| 欧美日韩瑟瑟在线播放| 久久精品夜夜夜夜夜久久蜜豆| a在线观看视频网站| av专区在线播放| 免费一级毛片在线播放高清视频| 精品久久久久久久久av| 欧美性猛交╳xxx乱大交人| 熟女人妻精品中文字幕| 99热这里只有精品一区| 伦理电影大哥的女人| 丰满的人妻完整版| av女优亚洲男人天堂| 国产色爽女视频免费观看| 一进一出抽搐gif免费好疼| 亚洲精品国产成人久久av| 别揉我奶头~嗯~啊~动态视频| 别揉我奶头 嗯啊视频| 免费在线观看影片大全网站| 人妻久久中文字幕网| 免费在线观看成人毛片| 十八禁网站免费在线| 国产人妻一区二区三区在| 韩国av在线不卡| 国产白丝娇喘喷水9色精品| 2021天堂中文幕一二区在线观| 亚洲av一区综合| 日本在线视频免费播放| 12—13女人毛片做爰片一| 美女高潮喷水抽搐中文字幕| 午夜激情福利司机影院| 伊人久久精品亚洲午夜| 亚州av有码| 亚洲人与动物交配视频| 少妇熟女aⅴ在线视频| 久久亚洲精品不卡| 午夜精品久久久久久毛片777| 伊人久久精品亚洲午夜| 美女高潮喷水抽搐中文字幕| 亚洲成a人片在线一区二区| 亚洲专区国产一区二区| 亚洲国产欧美人成| 91麻豆av在线| 夜夜看夜夜爽夜夜摸| 国产一区二区亚洲精品在线观看| 极品教师在线视频| 国内精品美女久久久久久| 亚洲av日韩精品久久久久久密| 亚洲国产精品久久男人天堂| 国产一级毛片七仙女欲春2| 91在线精品国自产拍蜜月| 中文字幕免费在线视频6| 悠悠久久av| 免费人成视频x8x8入口观看| 国产aⅴ精品一区二区三区波| 全区人妻精品视频| 综合色av麻豆| av黄色大香蕉| 又粗又爽又猛毛片免费看| 欧美激情在线99| 亚洲av中文字字幕乱码综合| 最新中文字幕久久久久| 国产精品久久视频播放| 欧美高清成人免费视频www| 国产蜜桃级精品一区二区三区| 国产麻豆成人av免费视频| 国产成人一区二区在线| 一进一出好大好爽视频| 乱系列少妇在线播放| 伊人久久精品亚洲午夜| 变态另类丝袜制服| 精品久久久久久成人av| 久久热精品热| 免费看美女性在线毛片视频| 午夜久久久久精精品| 国产一级毛片七仙女欲春2| 变态另类丝袜制服| 精品一区二区三区视频在线| 日本三级黄在线观看| 在线观看av片永久免费下载| 夜夜爽天天搞| 美女被艹到高潮喷水动态| 午夜福利欧美成人| 两个人视频免费观看高清| 一本精品99久久精品77| 精品人妻1区二区| www.www免费av| bbb黄色大片| 黄色日韩在线| 久久久精品欧美日韩精品| 尤物成人国产欧美一区二区三区| 91久久精品电影网| 毛片一级片免费看久久久久 | 亚洲成av人片在线播放无| 久久婷婷人人爽人人干人人爱| 午夜精品在线福利| 久久久久国产精品人妻aⅴ院| 女的被弄到高潮叫床怎么办 | 中亚洲国语对白在线视频| 精品久久久噜噜| 美女大奶头视频| 午夜爱爱视频在线播放| 老熟妇乱子伦视频在线观看| 噜噜噜噜噜久久久久久91| 亚洲精品456在线播放app | 男女视频在线观看网站免费| 欧美成人免费av一区二区三区| 欧美又色又爽又黄视频| 性色avwww在线观看| 黄色一级大片看看| 高清日韩中文字幕在线| 很黄的视频免费| 亚洲av五月六月丁香网| 久久精品国产99精品国产亚洲性色| 国产精品久久电影中文字幕| 看十八女毛片水多多多| 久久亚洲精品不卡| 久久久久久大精品| 男女之事视频高清在线观看| 一区二区三区免费毛片| 色综合婷婷激情| 99久久精品热视频| 老司机午夜福利在线观看视频| 中文字幕av成人在线电影| 日日干狠狠操夜夜爽| 在线播放国产精品三级| 国产 一区精品| 日本爱情动作片www.在线观看 | 九九爱精品视频在线观看| 日韩欧美国产一区二区入口| 久久久久九九精品影院| 乱人视频在线观看| 免费看av在线观看网站| 白带黄色成豆腐渣| 在线观看av片永久免费下载| 岛国在线免费视频观看| 国产精品一区二区免费欧美| 91av网一区二区| 亚洲欧美日韩高清专用| 变态另类成人亚洲欧美熟女| 狂野欧美激情性xxxx在线观看| 亚洲第一电影网av| 欧美日本亚洲视频在线播放| 久久天躁狠狠躁夜夜2o2o| 亚洲精华国产精华液的使用体验 | 欧美丝袜亚洲另类 | 亚洲av一区综合| 淫妇啪啪啪对白视频| 亚洲美女黄片视频| 久9热在线精品视频| 伦理电影大哥的女人| 99在线人妻在线中文字幕| 亚洲人与动物交配视频| 日本黄色视频三级网站网址| a级毛片免费高清观看在线播放| 老司机午夜福利在线观看视频| 国产高潮美女av| 岛国在线免费视频观看| 日日摸夜夜添夜夜添小说| 亚州av有码| 久久久国产成人精品二区| 老熟妇仑乱视频hdxx| 亚洲国产精品成人综合色| 国产v大片淫在线免费观看| av.在线天堂| 最近最新中文字幕大全电影3| 免费黄网站久久成人精品| 97碰自拍视频| 一a级毛片在线观看| 国产免费av片在线观看野外av| 男人的好看免费观看在线视频| 国产aⅴ精品一区二区三区波| 精华霜和精华液先用哪个| bbb黄色大片| 国产午夜福利久久久久久| 校园人妻丝袜中文字幕| 成人欧美大片| 国产色婷婷99| 国产女主播在线喷水免费视频网站 | 中文资源天堂在线| 亚洲av成人av| 九九爱精品视频在线观看| 国产乱人视频| 最近最新免费中文字幕在线| 中文资源天堂在线| 少妇的逼好多水| 内地一区二区视频在线| 日韩欧美精品v在线| 三级毛片av免费| 小说图片视频综合网站| 亚洲av熟女| 中文字幕高清在线视频| 热99在线观看视频| 国产成人aa在线观看| 观看免费一级毛片| 国产精品久久久久久亚洲av鲁大| 免费搜索国产男女视频| 国产精品乱码一区二三区的特点| 99热只有精品国产| 欧洲精品卡2卡3卡4卡5卡区| 国产女主播在线喷水免费视频网站 | 国产精品一区www在线观看 | 啪啪无遮挡十八禁网站| 色综合婷婷激情| АⅤ资源中文在线天堂| 亚洲av电影不卡..在线观看| 中文字幕av成人在线电影| 一个人观看的视频www高清免费观看| 日日夜夜操网爽| 国产高清激情床上av| 可以在线观看的亚洲视频| 亚洲欧美日韩高清在线视频| 亚洲精品日韩av片在线观看| 亚洲欧美清纯卡通| 午夜久久久久精精品| 欧美日韩精品成人综合77777| 一级黄片播放器| 人人妻人人澡欧美一区二区| 日本爱情动作片www.在线观看 | 性插视频无遮挡在线免费观看| 午夜精品一区二区三区免费看| 男女下面进入的视频免费午夜| 免费看日本二区| 老司机深夜福利视频在线观看| 日韩欧美一区二区三区在线观看| av女优亚洲男人天堂| 亚洲国产精品合色在线| 美女高潮喷水抽搐中文字幕| 日韩欧美国产在线观看| 一进一出好大好爽视频| 很黄的视频免费| 国产爱豆传媒在线观看| 亚洲成人久久爱视频| 男女啪啪激烈高潮av片| 熟女人妻精品中文字幕| 国产精品久久久久久久电影| 久久久久九九精品影院| 两性午夜刺激爽爽歪歪视频在线观看| 日本色播在线视频| 精品日产1卡2卡| 两个人的视频大全免费| 欧美性感艳星| 午夜久久久久精精品| 一级a爱片免费观看的视频| 毛片女人毛片| 国产亚洲精品av在线| 成人特级av手机在线观看| 黄色丝袜av网址大全| 日本色播在线视频| 欧美zozozo另类| 亚洲专区国产一区二区| 日本a在线网址| 亚洲国产日韩欧美精品在线观看| 俄罗斯特黄特色一大片| 亚洲四区av| 国产精品一区二区三区四区免费观看 | 看十八女毛片水多多多| АⅤ资源中文在线天堂| 久久精品国产鲁丝片午夜精品 | 欧美国产日韩亚洲一区| 久久中文看片网| 又黄又爽又刺激的免费视频.| 在线观看一区二区三区| 国产精品永久免费网站| 国产 一区精品| 午夜激情欧美在线| 黄色一级大片看看| 国产一区二区在线av高清观看| 亚洲成av人片在线播放无| 久久久久久大精品| 十八禁国产超污无遮挡网站| 国产黄a三级三级三级人| 伦理电影大哥的女人| bbb黄色大片| 久久精品国产自在天天线| 日日摸夜夜添夜夜添小说| 99国产精品一区二区蜜桃av| 中文字幕av成人在线电影| 内地一区二区视频在线| 狂野欧美激情性xxxx在线观看| 伦理电影大哥的女人| 三级毛片av免费| 日日撸夜夜添| 亚洲,欧美,日韩| 欧美成人一区二区免费高清观看| 人妻夜夜爽99麻豆av| 国产精品久久电影中文字幕| 亚洲精品粉嫩美女一区| 亚洲精品一卡2卡三卡4卡5卡| 久久久久性生活片| 国产精品,欧美在线| 精品久久久久久久人妻蜜臀av| 性插视频无遮挡在线免费观看| 国产精品1区2区在线观看.| 亚洲aⅴ乱码一区二区在线播放| 日韩一本色道免费dvd| 一卡2卡三卡四卡精品乱码亚洲| 有码 亚洲区| 国产乱人伦免费视频| 亚洲欧美激情综合另类| 亚洲av免费在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 在线免费十八禁| 搡老岳熟女国产| 国产探花在线观看一区二区| 亚洲欧美日韩卡通动漫| 成年女人看的毛片在线观看| 两人在一起打扑克的视频| 十八禁网站免费在线| 国产毛片a区久久久久| 在线免费观看不下载黄p国产 | 中文字幕久久专区| 日韩av在线大香蕉| 成人毛片a级毛片在线播放| 久久99热6这里只有精品| 在线观看一区二区三区| 免费在线观看日本一区| 亚洲av第一区精品v没综合| 大型黄色视频在线免费观看| 99热这里只有精品一区| 最近视频中文字幕2019在线8| 亚洲美女搞黄在线观看 | 毛片女人毛片| 欧美日韩精品成人综合77777| 免费一级毛片在线播放高清视频| 天天躁日日操中文字幕| av天堂在线播放| 大又大粗又爽又黄少妇毛片口| 中亚洲国语对白在线视频| 自拍偷自拍亚洲精品老妇| 啦啦啦韩国在线观看视频|