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

    南極特拉諾瓦灣下降風(fēng)特征

    2018-06-27 09:09:24孫啟振張占海丁卓銘沈輝張林孟上
    極地研究 2018年2期
    關(guān)鍵詞:難言諾瓦特拉

    孫啟振 張占海 丁卓銘 沈輝 張林 孟上

    (1物理海洋教育部重點(diǎn)實(shí)驗(yàn)室, 中國(guó)海洋大學(xué), 山東 青島 266100;2國(guó)家海洋環(huán)境預(yù)報(bào)中心, 極地環(huán)境研究預(yù)報(bào)室, 北京 100081;3國(guó)家海洋局極地科學(xué)重點(diǎn)實(shí)驗(yàn)室, 中國(guó)極地研究中心, 上海 200136)

    0 引言

    南極大陸近表層風(fēng)場(chǎng)的主要特征是: 由冰蓋表層冷卻作用產(chǎn)生的下降風(fēng)風(fēng)場(chǎng)在空間分布上極不均勻, 下降風(fēng)在大陸邊緣少數(shù)地區(qū)匯合, 形成風(fēng)向穩(wěn)定的強(qiáng)下降風(fēng)區(qū)。下降風(fēng)在這些地區(qū)的匯集現(xiàn)象是由當(dāng)?shù)鬲?dú)特的地形決定的[1-2]。南極下降風(fēng)在水平方向上綿延數(shù)千千米, 而在垂直方向上僅到達(dá)數(shù)百米的高度[3]。下降風(fēng)的風(fēng)場(chǎng)特征具有空間上的穩(wěn)定性和時(shí)間上的持續(xù)性, 早在20世紀(jì)60年代就引起了人們的研究興趣[4]; 20世紀(jì)80年代以來(lái), 對(duì)南極大尺度下降風(fēng)發(fā)生機(jī)制的研究越來(lái)越多。南極地區(qū)現(xiàn)場(chǎng)觀(guān)測(cè)資料少, 因此數(shù)值模式和衛(wèi)星遙感已經(jīng)成為研究下降風(fēng)動(dòng)力過(guò)程的重要手段[5-9],這些研究也涉及了下降風(fēng)的質(zhì)量輸送及其與中緯度大尺度環(huán)流之間的相互作用[10]。

    南極羅斯海西岸的特拉諾瓦灣(Terra Nova Bay,TNB)及其附近地區(qū)地形環(huán)境復(fù)雜, 氣候特征鮮明,是下降風(fēng)較強(qiáng)的地區(qū)之一[11]。來(lái)自南極大陸的下降風(fēng)從橫貫?zāi)蠘O山脈泄流而下, 流經(jīng)特拉諾瓦灣西岸的冰川, 其長(zhǎng)期的沖刷作用在該地區(qū)形成了面積廣闊的藍(lán)冰。強(qiáng)烈的下降風(fēng)將特拉諾瓦灣的海冰向東推動(dòng), 特拉諾瓦灣南側(cè)的Drygalski冰舌阻擋了更高緯度地區(qū)的海冰向北漂移, 使得該海灣即使在冬季也經(jīng)常處于無(wú)冰或少冰狀態(tài), 這是特拉諾瓦灣冰間湖在冬季形成和維持的重要原因[12-17]。這些地形特征從衛(wèi)星可見(jiàn)光圖像上清晰可辨(圖1)。

    圖1 南極羅斯海特拉諾瓦灣及其附近地區(qū)地形和衛(wèi)星可見(jiàn)光圖像.白色區(qū)域?yàn)槟蠘O冰蓋冰架, 曲線(xiàn)為間隔200 m的地形等高線(xiàn), 黑色區(qū)域?yàn)楸g湖, 藍(lán)色區(qū)域?yàn)榱_斯海Fig.1. The topography and visible satellite image of Terra Nova Bay and its surrounding areas, Antarctica. The ice sheet and ice shelf are indicated by white, polynyas by black, and the Ross Sea by blue. Black contours show topography with 200 m intervals

    自20世紀(jì)80年代起, 國(guó)際上就開(kāi)始了對(duì)南極特拉諾瓦灣下降風(fēng)的研究[18-20]。美國(guó)和意大利先后在特拉諾瓦灣西側(cè)的Priestley冰川、Reeves冰川、難言島(Inexpressible Island)和富蘭克林島(Franklin Island)等地安裝了自動(dòng)氣象站, 用以監(jiān)測(cè)這些特殊地形環(huán)境下降風(fēng)的分布情況[21-24]。這些研究工作歸納了特拉諾瓦灣西岸冰川的下降風(fēng)特點(diǎn), 發(fā)現(xiàn)下降風(fēng)在Reeves冰川東部匯合加速并流經(jīng)難言島地區(qū), 在空間和時(shí)間上具有很強(qiáng)的穩(wěn)定性[25]。從特拉諾瓦灣向東流向羅斯海的下降風(fēng)可以向東延續(xù)達(dá)200 km, 還可能引發(fā)中尺度氣旋[26]。

    除了利用現(xiàn)場(chǎng)觀(guān)測(cè)資料進(jìn)行局地氣候研究之外, 衛(wèi)星遙感、飛機(jī)觀(guān)測(cè)和數(shù)值模擬手段也被用于分析特拉諾瓦灣及其附近地區(qū)的大尺度下降風(fēng)特征。由于下降風(fēng)在垂直尺度上的混合作用不明顯,加之其下墊面海冰的熱力作用, 很難從衛(wèi)星圖像上對(duì)下降風(fēng)進(jìn)行辨識(shí)和分析[27]。飛機(jī)觀(guān)測(cè)[28]發(fā)現(xiàn), 在距離地面170 m高度的下降風(fēng)風(fēng)場(chǎng)仍然受到地形的強(qiáng)烈影響。對(duì)該地區(qū)下降風(fēng)的數(shù)值模擬研究[29-30]發(fā)現(xiàn), 特拉諾瓦灣地區(qū)的下降風(fēng)起源于距離沿岸約180 km處的南極內(nèi)陸, 下降風(fēng)在南極內(nèi)陸海拔大于1 500 m的地區(qū)風(fēng)速較小, 在大陸邊緣陡峭地區(qū)開(kāi)始加強(qiáng)。數(shù)值模擬研究還證實(shí)了下降風(fēng)在特拉諾瓦灣冰間湖的出現(xiàn)過(guò)程中發(fā)揮了重要作用[31]。

    特拉諾瓦灣由于其特殊的地理位置和豐富的科學(xué)研究資源, 吸引了國(guó)際上南極科學(xué)研究者的目光。繼意大利和韓國(guó)分別于1985年和2014年在特拉諾瓦灣建立考察站以來(lái), 我國(guó)也已經(jīng)開(kāi)始在該海灣沿岸的難言島建設(shè)新的考察站。我國(guó)南極考察隊(duì)于2012年年底和2016年年初在難言島安裝了兩套自動(dòng)氣象站, 用于長(zhǎng)期觀(guān)測(cè)該地區(qū)的天氣特征, 這些氣象觀(guān)測(cè)數(shù)據(jù)對(duì)于當(dāng)?shù)貧庀蟓h(huán)境研究具有重要作用[32-33]。然而, 難言島的自動(dòng)氣象站資料并不足以代表特拉諾瓦灣及其附近地區(qū)的下降風(fēng)特征, 因此對(duì)特拉諾瓦灣及其附近地區(qū)氣象要素場(chǎng)的分析顯得尤其重要。在南極地區(qū)采用高分辨率數(shù)值天氣預(yù)報(bào)模擬技術(shù), 有助于我們對(duì)特拉諾瓦灣地區(qū)的下降風(fēng)特征開(kāi)展深入研究。目前, 美國(guó)國(guó)家大氣研究中心和俄亥俄州立大學(xué)聯(lián)合開(kāi)發(fā)的南極中尺度預(yù)報(bào)系統(tǒng)(Antarctic Mesoscale Prediction System,AMPS)在南極天氣數(shù)值預(yù)報(bào)方面處于國(guó)際領(lǐng)先地位[34-36]。AMPS自2000年開(kāi)始為美國(guó)的南極科研計(jì)劃和國(guó)際南極科考隊(duì)提供天氣預(yù)報(bào)支持, 經(jīng)過(guò)十幾年的發(fā)展, AMPS可提供的南極地區(qū)短期天氣預(yù)報(bào)產(chǎn)品的水平分辨率不斷提高, 其多年存檔數(shù)據(jù)也曾被用于南極大尺度下降風(fēng)研究[37]。自2013年5月以來(lái), AMPS在特拉諾瓦灣地區(qū)的數(shù)值天氣預(yù)報(bào)水平分辨率提升為1.1 km, 這為分析該地區(qū)下降風(fēng)風(fēng)場(chǎng)提供了可能。

    我國(guó)即將在羅斯海西岸建設(shè)的常年考察站即位于特拉諾瓦灣沿岸的難言島上, 因此很有必要對(duì)該地區(qū)的下降風(fēng)時(shí)空分布特征做一些分析研究。本文采用南極中尺度預(yù)報(bào)系統(tǒng)高分辨率資料和特拉諾瓦灣難言島自動(dòng)氣象站實(shí)測(cè)數(shù)據(jù), 分析特拉諾瓦灣及其附近地區(qū)的下降風(fēng)特征。

    1 數(shù)據(jù)

    1.1 模式預(yù)報(bào)數(shù)據(jù)

    AMPS是一個(gè)針對(duì)南極地區(qū)的實(shí)時(shí)數(shù)值天氣預(yù)報(bào)系統(tǒng), 其核心是Polar WRF極地?cái)?shù)值天氣預(yù)報(bào)模式[38], 自2000年10月開(kāi)始為美國(guó)的南極科研計(jì)劃以及國(guó)際南極科考提供天氣預(yù)報(bào)支持。美國(guó)國(guó)家大氣研究中心負(fù)責(zé)每天兩次運(yùn)行AMPS系統(tǒng), 為國(guó)際預(yù)報(bào)團(tuán)體提供模式產(chǎn)品、存檔預(yù)報(bào)產(chǎn)品并維護(hù)模式代碼; 俄亥俄州立大學(xué)伯德極地研究中心負(fù)責(zé)繼續(xù)發(fā)展Polar WRF模式的物理過(guò)程并評(píng)估模式性能。根據(jù)目前的數(shù)值預(yù)報(bào)模式配置方案, AMPS在南極地區(qū)設(shè)置5個(gè)嵌套預(yù)報(bào)區(qū)域, 其中最外層(Domain 1)分辨率為30 km, 包含南極大陸的預(yù)報(bào)區(qū)域(Domain 2)分辨率為10 km, 分辨率最高的區(qū)域(Domain 5)以羅斯島為中心, 覆蓋了羅斯冰架西北部分、維多利亞地東岸, 并向北延伸至特拉諾瓦灣地區(qū)(圖2)。

    本文選取的數(shù)據(jù)為AMPS自2014年1月—2015年12月Domain 3和Domain 5的3 h預(yù)報(bào)數(shù)據(jù)和月平均數(shù)據(jù), 包含10 m高度風(fēng)、2 m高度氣溫和地面氣溫、地面氣壓等要素場(chǎng)數(shù)據(jù), 并選取2016年7月10日06 UTC Domain 5預(yù)報(bào)數(shù)據(jù)分析下降風(fēng)空間結(jié)構(gòu)。

    圖2 南極中尺度預(yù)報(bào)系統(tǒng)(AMPS)的部分區(qū)域設(shè)置示意圖. 本文采用的資料來(lái)自Domain 3和Domain 5, 分辨率分別為3.3 km和1.1 km; 填色表示地形高度, 南極內(nèi)陸海拔高度按照紅色、白色、棕色、黃色和綠色的順序逐漸減小;Fig.2. The domains of Antarctic Mesoscale Prediction System (AMPS). Archive data from Domain 3 and Domain 5 are used in this study, and the resolution of Domain 3 and Domain 5 are 3.3 km and 1.1 km, respectively. Terrain heights are indicated by shaded colors, decreasing with the order of red, white, brown, yellow and green

    1.2 現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)

    自20世紀(jì)50年代至今, 美國(guó)威斯康星大學(xué)在南極地區(qū)陸續(xù)建立了多套自動(dòng)氣象站, 形成了南極地面氣象自動(dòng)觀(guān)測(cè)網(wǎng)[39]。其中, 難言島南端的Manuela自動(dòng)氣象站(74.946°S, 163.687°E)建立于1984年。該自動(dòng)氣象站海拔高度78 m, 觀(guān)測(cè)要素有氣壓、氣溫、風(fēng)速、風(fēng)向和相對(duì)濕度等, 各傳感器距離地面高度為1.5 m, 采樣間隔為10 min。該站曾于1987年1月觀(guān)測(cè)到極端強(qiáng)風(fēng)風(fēng)速46.9 m·s–1,于1992年9月觀(guān)測(cè)到極端低溫–42.4℃。本文采用該站自2014年1月—2017年1月實(shí)測(cè)的1 h平均數(shù)據(jù)。

    2 結(jié)果分析

    2.1 難言島Manuela站氣候特征

    2.1.1 日循環(huán)特征

    難言島長(zhǎng)寬僅分別為14 km和4 km, 面積不足40 km2, 但其地形復(fù)雜, 單點(diǎn)測(cè)站的數(shù)據(jù)并不足以代表整個(gè)難言島的總體氣候狀況, 風(fēng)速的空間分布尤其如此。本文以2016年7月和2017年1月分別代表冬季和夏季, 以Manuela站的風(fēng)速、風(fēng)向、氣溫、相對(duì)濕度日平均曲線(xiàn)(圖3)來(lái)分析該站在這兩個(gè)季節(jié)的氣象要素日循環(huán)特征。該曲線(xiàn)圖中短豎線(xiàn)代表該月份31 d統(tǒng)計(jì)值的標(biāo)準(zhǔn)誤差。Manuela站在冬季7月份觀(guān)測(cè)到的風(fēng)向穩(wěn)定為西西南, 風(fēng)力較大, 約為7—8級(jí), 風(fēng)速平均為16.3 m·s-1,統(tǒng)計(jì)顯示7月份不同日期的風(fēng)速相差較大。夏季1月份的風(fēng)速日變化較大, 當(dāng)?shù)貢r(shí)間傍晚至夜間(世界時(shí)08:00—12:00)風(fēng)速較低, 僅有5 m·s–1左右, 中午前后風(fēng)速稍大, 風(fēng)力最高為4—5級(jí)。1月風(fēng)力較小時(shí), 不同日期的風(fēng)向有顯著不同,但基本集中在南至西南這個(gè)風(fēng)向范圍; 風(fēng)力較大時(shí), 風(fēng)向與冬季一致, 穩(wěn)定為西至西南。

    與風(fēng)速的季節(jié)特征類(lèi)似, 不同季節(jié)的氣溫日平均變化特征差異較大。冬季7月不同日期的氣溫相差較大, 但平均氣溫日變化較小, 日溫差約為1℃, 這主要是因?yàn)闃O夜期間沒(méi)有日照。夏季1月份極晝期間, 氣溫日變化顯著, 凌晨氣溫最低,傍晚氣溫最高, 日溫差可達(dá)4℃。從圖3可見(jiàn), 氣溫在冬季沒(méi)有顯著的日變化, 因?yàn)槎練鉁嘏c風(fēng)速?zèng)]有較強(qiáng)的相關(guān)性。夏季1月份Manuela站的氣溫與風(fēng)速均有顯著的日變化特點(diǎn), 風(fēng)速變化比氣溫變化滯后了約3 h, 例如日最高風(fēng)速出現(xiàn)時(shí)間比日最低氣溫滯后約3 h。不同季節(jié)的相對(duì)濕度日變化特征相似, 均表現(xiàn)出了微弱的日循環(huán)特征,當(dāng)?shù)貢r(shí)間中午前后的濕度較低, 夜間至凌晨濕度較高。夏季的濕度日變化幅度比冬季高, 夏季當(dāng)?shù)貢r(shí)間凌晨濕度較高, 比下午高約9%, 而冬季日濕度差約為5%, 濕度最高值出現(xiàn)在夜間。

    2.1.2 季節(jié)循環(huán)特征

    AMPS在覆蓋難言島的Domain 5區(qū)域分辨率高達(dá)1.1 km, 可以較好地模擬當(dāng)?shù)靥鞖夂蜌夂蛱卣?。此處取AMPS數(shù)據(jù)在Manuela站周?chē)?個(gè)格點(diǎn)的雙線(xiàn)性插值, 與自動(dòng)氣象站觀(guān)測(cè)數(shù)據(jù)進(jìn)行比較, 兩組數(shù)據(jù)時(shí)間范圍均為2014年1月—2015年12月。圖4為該時(shí)間段的月平均值。

    由圖4和表1可見(jiàn), AMPS對(duì)難言島地區(qū)的風(fēng)速和氣溫有較好的模擬能力??傮w而言, AMPS在難言島地區(qū)冬季的模擬風(fēng)速稍高于實(shí)測(cè)值, 而夏季模擬值稍低于實(shí)測(cè)值, 但平均偏差僅為1.79 m·s–1。冬季氣溫模擬值稍低于實(shí)測(cè)值, 夏季氣溫模擬接近實(shí)測(cè)值。難言島地區(qū)風(fēng)速?gòu)?月開(kāi)始迅速增大,至4月份達(dá)到全年最高值, 接著在長(zhǎng)達(dá)5個(gè)月的時(shí)間里, 風(fēng)速均處于較高的范圍, 此時(shí)段內(nèi)風(fēng)力平均為8級(jí)以上, 9月份之后, 風(fēng)速逐漸減弱。與之相對(duì)應(yīng), 氣溫從1月份開(kāi)始迅速降低, 4—9月的氣溫處于–25—–30℃且變化幅度較小, 9月之后氣溫開(kāi)始升高。

    難言島面積較小但地形復(fù)雜, 其東南部分是較為開(kāi)闊的平坦地形, 西北部分有海拔約400 m的山峰。AMPS在此地區(qū)雖然水平分辨率較高, 但仍然無(wú)法精確描述當(dāng)?shù)貜?fù)雜的地形,因而對(duì)風(fēng)向模擬的偏差不可避免。由圖4可見(jiàn),AMPS模擬得到的風(fēng)向比實(shí)測(cè)值偏西約30°,在夏季風(fēng)速相對(duì)較小時(shí), 風(fēng)向的模擬值偏差更大。相對(duì)濕度是數(shù)值預(yù)報(bào)的難點(diǎn)之一, 尤其是極端氣象環(huán)境下的沿海地區(qū)。雖然AMPS對(duì)濕度的模擬能力弱于風(fēng)速和氣溫, 但仍然處于可信賴(lài)的水平。

    圖4 2014年1月—2015年12月難言島Manuela自動(dòng)氣象站的各要素季節(jié)平均曲線(xiàn)圖. 實(shí)線(xiàn)為AMPS的1.1 km數(shù)據(jù)在Manuela站的插值結(jié)果; 虛線(xiàn)為Manuela自動(dòng)氣象站的觀(guān)測(cè)數(shù)據(jù)Fig.4. Seasonal mean time series of wind speed, direction, temperature and humidity of Manuela AWS at Inexpressible Island from Jan. 2014 to Dec. 2015. Solid: AMPS 1.1 km; Dashed:AWS

    表1 難言島Manuela自動(dòng)氣象站所在位置的AMPS 1.1 km高分辨率模擬結(jié)果統(tǒng)計(jì)值(2014年1月—2015年12月)Table 1. Statistics of AMPS 1.1km simulations at the position of Manuela AWS on Inexpressible Island(from January 2014 to December 2015)

    2.2 特拉諾瓦灣地區(qū)風(fēng)場(chǎng)特征

    上節(jié)以難言島的自動(dòng)氣象站觀(guān)測(cè)數(shù)據(jù)為主,分析了局地氣候特征。本節(jié)討論特拉諾瓦灣及其附近地區(qū)的大尺度下降風(fēng)特征。

    2.2.1 水平分布特征

    圖5為2014—2015年季節(jié)平均的羅斯海西岸近地層風(fēng)場(chǎng)特征??傮w來(lái)看, 該地區(qū)風(fēng)向的空間分布在不同季節(jié)有相似的特征, 但風(fēng)速有顯著的季節(jié)變化。近地層風(fēng)從橫貫?zāi)蠘O山脈和維多利亞地以西海拔較高的冰蓋流向沿岸, 冬季在特拉諾瓦灣西岸等幾處冰川地帶形成強(qiáng)風(fēng)匯集區(qū)。特拉諾瓦灣的北部和維多利亞地東岸的風(fēng)速相對(duì)較小。羅斯島地勢(shì)較高, 在地形阻擋作用下, 其以北的羅斯冰架地區(qū)形成了微弱的中尺度系統(tǒng), 但風(fēng)速較小。在夏季, 無(wú)論是冰蓋內(nèi)陸還是沿岸的強(qiáng)風(fēng)區(qū), 風(fēng)速都顯著減小, 但風(fēng)場(chǎng)的形態(tài)與冬季相比變化不大。強(qiáng)風(fēng)匯集區(qū)從特拉諾瓦灣向東延伸至羅斯海內(nèi)的范圍,在夏季明顯縮小。

    圖5 羅斯海西岸地區(qū)2014年1月—2015年12月季節(jié)性平均風(fēng)場(chǎng). 數(shù)據(jù)源于AMPS 3.3 km分辨率數(shù)據(jù), MAM: 3—5月, JJA: 6—8月, SON: 9—11月, DJF: 12—2月Fig.5. Seasonal mean wind pattern of the western coast of the Ross Sea from January 2014 to December 2015. Data from AMPS 3.3 km, MAM: from March to May, JJA: from June to August, SON: from September to November, DJF: from December to February

    圖6為特拉諾瓦灣及其附近地區(qū)的2014—2015年近地層月平均風(fēng)場(chǎng)。來(lái)自南極大陸冰蓋斜坡的強(qiáng)下降風(fēng)主要匯集區(qū)位于Reeves冰川,Reeves冰川的下降風(fēng)在到達(dá)難言島附近地區(qū)時(shí)風(fēng)速達(dá)到最大, 次要匯集區(qū)位于David冰川(各冰川位置見(jiàn)圖1)。這些冰川地區(qū)之所以成為下降風(fēng)匯集區(qū), 一是從圖1的地形等高線(xiàn)可見(jiàn), Reeves冰川、David冰川及Larsen冰川等地區(qū)西部的高度落差較大, 下降風(fēng)氣流從冰蓋內(nèi)陸可以加速流下;二是冰川內(nèi)部下墊面平坦, 且在冰川東部?jī)蓚?cè)的峽谷逐漸狹窄, 這種獨(dú)特的地形使得下降風(fēng)在前進(jìn)過(guò)程中不斷加速, 峽谷效應(yīng)也增強(qiáng)了近地面風(fēng)速。下降風(fēng)氣流長(zhǎng)期流經(jīng)這些冰川, 在雪粒的摩擦作用下, 形成了從衛(wèi)星圖像上可分辨的藍(lán)冰。從圖6還可以發(fā)現(xiàn), 下降風(fēng)氣流流經(jīng)上述冰川時(shí),其兩側(cè)邊緣界限非常明顯, 表明邊緣地區(qū)的風(fēng)速梯度較大。

    從季節(jié)變化特征來(lái)說(shuō), 強(qiáng)風(fēng)區(qū)的風(fēng)速在11月至次年1月較小, 2月和10月分別為風(fēng)速增大和減小的過(guò)渡期, 3—9月風(fēng)速較大。不同月份的風(fēng)場(chǎng)空間分布形態(tài)相似, 風(fēng)向穩(wěn)定, 可見(jiàn)該地區(qū)下降風(fēng)的形態(tài)特征主要來(lái)源于周邊地形的作用和下降風(fēng)自身的形成機(jī)制, 很少受南大洋天氣尺度系統(tǒng)的影響。圖6表明, 來(lái)自Reeves冰川和David冰川的強(qiáng)風(fēng)向東延伸至特拉諾瓦灣以東海域, 其延伸范圍隨季節(jié)而變化, 冬季可達(dá)100 km以上,夏季風(fēng)力較弱時(shí)僅能向東延伸30—50 km。Drygalski冰舌東南方向的羅斯海海域基本不受下降風(fēng)影響, 風(fēng)速常年較小。

    圖6 特拉諾瓦灣及其附近地區(qū)2014年1月—2015年12月的月平均風(fēng)場(chǎng). 數(shù)據(jù)源于AMPS 1.1 km分辨率數(shù)據(jù)Fig.6. Monthly mean wind pattern of the Terra Nova Bay(TNB) and its surrounding areas from January 2014 to December 2015. Data from AMPS 1.1 km

    為了分析特拉諾瓦灣及其附近地區(qū)風(fēng)向隨著時(shí)間變化的程度, 計(jì)算了無(wú)量綱的風(fēng)向穩(wěn)定度(wind directional constancy)DC, 即平均矢量風(fēng)速與標(biāo)量風(fēng)速的比值。該風(fēng)向穩(wěn)定度值為1時(shí)表示風(fēng)向完全穩(wěn)定。風(fēng)向穩(wěn)定度由如下公式計(jì)算得到:

    其中,

    根據(jù)圖3不同季節(jié)的風(fēng)速日變化曲線(xiàn), 選取2015年1月各日21:00 UTC、09:00 UTC和2015年7月各日06:00 UTC、15:00 UTC數(shù)據(jù)進(jìn)行計(jì)算,分別代表夏季風(fēng)速較大時(shí)和較小時(shí)、冬季風(fēng)速較大時(shí)和較小時(shí)的風(fēng)向穩(wěn)定度。圖7為4個(gè)時(shí)次的日平均風(fēng)向穩(wěn)定度。由該圖可見(jiàn), 夏季1月份風(fēng)速較大時(shí), 上述特拉諾瓦灣地區(qū)各冰川強(qiáng)風(fēng)匯集區(qū)及其西部大陸冰蓋上風(fēng)向較為穩(wěn)定, 1月份風(fēng)速較小時(shí), 風(fēng)向分布不集中, 即夏季風(fēng)速變化時(shí),各地區(qū)風(fēng)向有明顯變化。而在冬季7月份, 無(wú)論風(fēng)速相對(duì)較大還是較小, 圖中大部分地區(qū)的風(fēng)向穩(wěn)定度都接近于1。值得注意的是, 圖3和圖6已經(jīng)表明, 冬季風(fēng)速的最低值仍然高于夏季風(fēng)速的最高值。由此可見(jiàn), 特拉諾瓦灣及其附近地區(qū)的風(fēng)速較大時(shí), 風(fēng)向較為穩(wěn)定。

    為了更加精細(xì)地展現(xiàn)不同季節(jié)難言島的風(fēng)場(chǎng)特點(diǎn), 利用1.1 km水平分辨率的2014年1月—2015年12月AMPS數(shù)據(jù)繪制了不同季節(jié)難言島及其附近地區(qū)的平均風(fēng)場(chǎng)形態(tài)。圖8表明, 難言島處于來(lái)自Reeves冰川的強(qiáng)下降風(fēng)氣流的核心區(qū)。由于難言島最高海拔達(dá)400 m以上, 強(qiáng)下降風(fēng)氣流在到達(dá)難言島之前風(fēng)速有所減弱, 但越過(guò)難言島之后, 氣流再次加強(qiáng)。夏季(12—2月)該地區(qū)下降風(fēng)較弱時(shí), 這種地形阻擋作用不明顯。

    圖7 特拉諾瓦灣地區(qū)2015年1月和7月風(fēng)向穩(wěn)定度分布圖. 數(shù)據(jù)源于AMPS 1.1 km分辨率數(shù)據(jù); a),c)當(dāng)月風(fēng)速較大時(shí)的風(fēng)向穩(wěn)定度; b),d)當(dāng)月風(fēng)速較小時(shí)的風(fēng)向穩(wěn)定度Fig .7. The mean wind directional constancy (DC) of TNB and its surrounding areas of Jan. and July 2015. Data from AMPS 1.1 km; a) and c), DC of high wind apeed; b) and d), DC of low wind speed

    地面位勢(shì)溫度場(chǎng)的空間分布特征可以反映出近地層大氣動(dòng)力和熱力過(guò)程。圖9為各季節(jié)特拉諾瓦灣及其附近地區(qū)的地面位勢(shì)溫度場(chǎng)分布圖。在下降風(fēng)較強(qiáng)的冬季(6—8月), Reeves冰川、Larsen冰川和David冰川所在地區(qū)的地面位勢(shì)溫度在空間上有著近似均一的分布特征, 即下降風(fēng)從冰川西側(cè)開(kāi)始匯集到流經(jīng)冰川并加強(qiáng)的過(guò)程中, 地面位勢(shì)溫度變化幅度很小。這表明下降風(fēng)在從冰川西部海拔較高處到沿岸地區(qū)海拔較低處的流動(dòng)是干絕熱過(guò)程。位勢(shì)溫度場(chǎng)的這種分布形態(tài)還表明下降風(fēng)主要取決于其內(nèi)部的動(dòng)力過(guò)程。在特拉諾瓦灣存在位勢(shì)溫度較大區(qū), 是因?yàn)榇颂幱忻娣e較大的冰間湖。

    由圖6可知, 特拉諾瓦灣附近冰川地區(qū)在冬季下降風(fēng)的風(fēng)速較大、范圍較廣。選取2016年7月10日06:00 UTC的AMPS 1.1 km分辨率預(yù)報(bào)數(shù)據(jù)為例, 分析Reeves冰川冬季強(qiáng)下降風(fēng)匯集區(qū)的風(fēng)速剖面特征。圖10a為A點(diǎn)至B點(diǎn)連線(xiàn)的風(fēng)速剖面圖, 即沿著下降風(fēng)流向的剖面; 圖10b為C點(diǎn)和D點(diǎn)連線(xiàn)的風(fēng)速剖面圖, 即下降風(fēng)風(fēng)帶的橫截面, ABCD四點(diǎn)的位置見(jiàn)圖6g。

    圖8 難言島及其附近地區(qū)2014年1月—2015年12月近地面風(fēng)場(chǎng)季節(jié)平均圖. 數(shù)據(jù)源于AMPS 1.1 km分辨率數(shù)據(jù);MAM: 3—5月, JJA: 6—8月, SON: 9—11月, DJF: 12—2月. 難言島位于圖的中心位置Fig.8. Monthly mean wind pattern of the Inexpressible Island and its surrounding areas from January 2014 to December 2015.Data is extracted from AMPS 1.1 km dataset; MAM: from March to May, JJA: from June to August, SON: from September to November, DJF: from December to February. The Inexpressible Island locates in the center of the plots

    圖9 特拉諾瓦灣地區(qū)2014年和2015年6—8月平均的地表位勢(shì)溫度場(chǎng). 數(shù)據(jù)來(lái)自AMPS 1.1 km分辨率數(shù)據(jù)Fig.9. Mean potential temperature of TNB and its surrounding areas from June to August in 2014 and 2015.Data from AMPS 1.1 km

    2.2.2 垂直剖面特征

    由圖10可見(jiàn), 在空間結(jié)構(gòu)上, Reeves冰川下降風(fēng)匯集區(qū)的上邊界有清晰的分界線(xiàn), 風(fēng)速較強(qiáng)區(qū)的界線(xiàn)大約位于第11—12模式eta 層, 距離地面高度約為650—800 m, 其中, 圖10a中位于163.7°E附近的風(fēng)速減弱過(guò)程是由于難言島山峰的阻擋作用,這與圖8一致。沿著下降風(fēng)的流向(圖10a), 風(fēng)速逐漸增大, 至冰川最窄處風(fēng)速達(dá)到最大, 此后逐漸減弱。在風(fēng)速的垂向分布上, 下降風(fēng)在不同地區(qū)的風(fēng)速最大值高度變化不大, 幾乎都位于距離地面50—200 m高度上。從橫截面圖(圖10b)看,Reeves冰川下降風(fēng)的匯集區(qū)存在強(qiáng)風(fēng)速核心區(qū), 南北兩側(cè)均有明顯的界線(xiàn), 其中北側(cè)有陡峭地形阻擋,南側(cè)與Larsen冰川下降風(fēng)匯集區(qū)連接。

    圖10 2016年7月10日06:00 UTC特拉諾瓦灣西岸Reeves冰川下降風(fēng)匯集區(qū)空間剖面圖. a)A點(diǎn)至B點(diǎn)風(fēng)速剖面圖;b)C點(diǎn)和D點(diǎn)風(fēng)速剖面圖(ABCD四點(diǎn)位置見(jiàn)圖6)Fig.10. The profile of katabatic winds over Reeves glacier to the west of TNB at 06:00 UTC on 10 July 2016. a) Profile from A to B; b) Profile from C to D. The positions of A, B, C, and D are shown in Fig.6

    3 結(jié)論

    下降風(fēng)是南極特拉諾瓦灣及其附近地區(qū)的主要?dú)夂蛱卣髦?。?duì)該地區(qū)下降風(fēng)的研究將有助于為我國(guó)南極考察隊(duì)提供可靠的氣象保障。本文采用南極中尺度預(yù)報(bào)系統(tǒng)高分辨率資料和特拉諾瓦灣難言島自動(dòng)氣象站實(shí)測(cè)數(shù)據(jù), 分析了特拉諾瓦灣及其附近地區(qū)的下降風(fēng)特征。

    分析發(fā)現(xiàn), 難言島地區(qū)風(fēng)速?gòu)?月開(kāi)始迅速增大, 自4月份開(kāi)始的長(zhǎng)達(dá)6個(gè)月的時(shí)間里, 風(fēng)速均處于較高的范圍, 此時(shí)段內(nèi)風(fēng)力平均為8級(jí)以上。9月份之后, 風(fēng)速逐漸減弱。與之相對(duì)應(yīng), 4—9月的氣溫介于–30—–25℃之間, 變化幅度較小; 氣溫從1月開(kāi)始迅速降低, 9月開(kāi)始升高。夏季1月份Manuela站的氣溫與風(fēng)速均有顯著的日變化, 風(fēng)速的變化滯后于氣溫變化約3 h。Manuela站在冬季7月份觀(guān)測(cè)到的風(fēng)向穩(wěn)定為西至西南, 風(fēng)力較大, 約為7—8級(jí), 風(fēng)速平均為16.3 m·s-1。南極中尺度預(yù)報(bào)系統(tǒng)AMPS對(duì)難言島地區(qū)的風(fēng)速和氣溫有較好的模擬能力, 但模擬的風(fēng)向比實(shí)測(cè)值偏西約30°。

    羅斯海西岸地區(qū)風(fēng)向穩(wěn)定, 其空間分布特征基本不隨季節(jié)而變化, 風(fēng)向穩(wěn)定度的統(tǒng)計(jì)顯示,特拉諾瓦灣及其附近地區(qū)的風(fēng)速較大時(shí), 風(fēng)向更為穩(wěn)定。但風(fēng)速有顯著的季節(jié)性演變特征。近地層下降風(fēng)從橫貫?zāi)蠘O山脈和維多利亞地以西海拔較高的冰蓋流向沿岸, 冬季在特拉諾瓦灣西岸等幾處冰川地帶形成強(qiáng)風(fēng)匯集區(qū), 這些匯集區(qū)主要位于Reeves冰川、Larsen冰川和David冰川, 其中難言島處于Reeves冰川下降風(fēng)匯集區(qū)核心區(qū)。下降風(fēng)強(qiáng)風(fēng)區(qū)的風(fēng)速在11月至次年1月較小, 2月和10月分別為風(fēng)速增大和減小的過(guò)渡期, 3—9月風(fēng)速較大。

    從空間結(jié)構(gòu)來(lái)說(shuō), 特拉諾瓦灣及其附近地區(qū)下降風(fēng)匯集區(qū)的上邊界和南北兩側(cè)均有清晰的分界線(xiàn), 風(fēng)速較強(qiáng)區(qū)從地面延伸至650—800 m高處。沿著下降風(fēng)的前進(jìn)方向, 下降風(fēng)風(fēng)速逐漸增大, 至冰川最窄處風(fēng)速達(dá)到最大。在風(fēng)速的垂向分布上, 下降風(fēng)的匯集區(qū)存在強(qiáng)風(fēng)速核心區(qū), 風(fēng)速最大值位于距離地面50—200 m高度上。由于難言島最高海拔達(dá)400 m以上, 來(lái)自Reeves冰川的強(qiáng)下降風(fēng)氣流在到達(dá)難言島之前風(fēng)速有所減弱,但越過(guò)難言島之后, 氣流再次加強(qiáng)。下降風(fēng)從冰川西側(cè)開(kāi)始匯集到流經(jīng)冰川并加強(qiáng)的過(guò)程中, 地面位勢(shì)溫度變化幅度很小, 表明下降風(fēng)在從內(nèi)陸高原到沿岸地區(qū)的流動(dòng)是干絕熱過(guò)程。該地區(qū)下降風(fēng)的形態(tài)特征主要來(lái)源于周邊地形的作用和下降風(fēng)自身的形成機(jī)制, 而很少受南大洋天氣尺度系統(tǒng)的影響。

    1 PARISH T R, CASSANO J J. Diagnosis of the katabatic wind influence on the wintertime Antarctic surface wind field from numerical simulations[J]. Monthly Weather Review, 2003, 131(6): 1128—1139.

    2 孫啟振, 張林, 張占海, 等. 南極中山站夏季下降風(fēng)數(shù)值模擬個(gè)例研究[J]. 海洋學(xué)報(bào), 2016, 38(3): 71—81.

    3 BROMWICH D H, PARISH T R, PELLEGRINI A, et al. Spatial and temporal characteristics of the intense katabatic winds at Terra Nova Bay, Antarctica[M]//BROMWICH D H, STEARNS C R. Antarctic Meteorology and Climatology: Studies Based on Automatic Weather Stations. Washington DC: American Geophysical Union, 1993.

    4 STRETEN N A. Some observations of Antarctic katabatic winds[J]. Australian Meteorological Magazine, 1963, 42: 1—23.

    5 PARISH T R, BROMWICH D H. Continental-scale simulation of the Antarctic katabatic wind regime[J]. Journal of Climate, 1991, 4(2):135—146.

    6 PARISH T R. Surface airflow over east Antarctica[J]. Monthly Weather Review, 1982, 110(2): 84—90.

    7 PARISH T R. A numerical study of strong katabatic winds over Antarctica[J]. Monthly Weather Review, 1984, 112(3): 545—554.

    8 REMY F, BROSSIER C, MINSTER J F. Intensity of satellite radar-altimeter return power over continental ice: a potential measurement of katabatic wind intensity[J]. Journal of Glaciology, 1990, 36(123): 133—142.

    9 REMY F, LEDROIT M, MINSTER J F. Katabatic wind intensity and direction over Antarctica derived from scatterometer data[J].Geophysical Research Letters, 1992, 19(10): 1021—1024.

    10 PARISH T R, BROMWICH D H. A case study of Antarctic katabatic wind interaction with large-scale forcing[J]. Monthly Weather Review, 1998, 126(1): 199—209.

    11 BROMWICH D H, KURTZ D D. Katabatic wind forcing of the Terra Nova Bay polynya[J]. Journal of Geophysical Research: Oceans,1984, 89(C3): 3561—3572.

    12 MATHIOT P, BARNIER B, GALLéE H, et al. Introducing katabatic winds in global ERA40 fields to simulate their impacts on the Southern Ocean and sea-ice[J]. Ocean Modelling, 2010, 35(3): 146—160.

    13 MASSOM R A, HARRIS P T, MICHAEL K J, et al. The distribution and formative processes of latent-heat polynyas in East Antarctica[J]. Annals of Glaciology, 1998, 27: 420—426.

    14 PARMIGGIANI F. Is the katabatic wind the forcing factor of Terra Nova Bay polynya events?[C]//Proceedings of the SPIE Volume 8175, Remote Sensing of the Ocean, Sea Ice, Coastal Waters, and Large Water Regions 2011. Prague, Czech Republic: SPIE, 2011,8175: 817511.

    15 VAN WOERT M L. Wintertime expansion and contraction of the Terra Nova Bay polynya[M]//SPEZIE G, MANZELLA G M R.Oceanography of the Ross SeaAntarctica. Milano: Springer, 1999: 145—164.

    16 BROMWICH D H, KURTZ D D. Experiences of Scott's Northern Party: evidence for a relationship between winter katabatic winds and the Terra Nova Bay polynya[J]. Polar Record, 1982, 21(131): 137—146.

    17 PARMIGGIANI F. Multi-year measurement of Terra Nova Bay winter polynya extents[J]. The European Physical Journal Plus, 2011,126(4): 39.

    18 KURTZ D D, BROMWICH D H. A recurring, atmospherically forced polynya in Terra Nova Bay[M]//JACOBS S S. Oceanology of the Antarctic Continental Shelf. Washington DC: American Geophysical Union, 1985.

    19 BROMWICH D H, KURTZ D D. Katabatic wind forcing of the Terra Nova Bay polynya[J]. Journal of Geophysical Research: Oceans,1984, 89(C3): 3561—3672.

    20 KURTZ D D, BROMWICH D H. Satellite observed behavior of the Terra Nova Bay polynya[J]. Journal of Geophysical Research:Oceans, 1983, 88(C14): 9717—9722.

    21 BROMWICH D H. An extraordinary katabatic wind regime at Terra Nova Bay, Antarctica[J]. Monthly Weather Review, 1989, 117(3):688—695.

    22 BROMWICH D H, PARISH T R, PELLEGRINI A. The katabatic wind regime near Terra Nova Bay, Antarctica[J]. Antarctic Journal,1990, 25(5): 267—269.

    23 PARISH T R, BROMWICH D H. Automatic weather station observations of strong katabatic winds near Terra Nova Bay, Antarctica[J].Antarctic Journal, 1991, 26(5): 265—267.

    24 DAVOLIO S, BUZZI A. Mechanisms of Antarctic katabatic currents near Terra Nova Bay[J]. Tellus A: Dynamic Meteorology and Oceanography, 2002, 54(2): 187-204, doi: 10.3402/tellusa.v54i2.12132.

    25 BROMWICH D H, PARISH T R, PELLEGRINI A, et al. Spatial and temporal characteristics of the intense katabatic winds at Terra Nova Bay, Antarctica[M]//BROMWICH D H, STEARNS C R. Antarctic Meteorology and Climatology: Studies Based on Automatic Weather Stations. Washington DC: American Geophysical Union, 1993.

    26 BROMWICH D H. Mesoscale cyclogenesis over the southwestern Ross Sea linked to strong katabatic winds[J]. Monthly Weather Review, 1991, 119(7): 1736—1753.

    27 BROMWICH D H. Satellite analyses of Antarctic katabatic wind behavior[J]. Bulletin of the American Meteorological Society, 1989,70(7): 738—749.

    28 PARISH T R, BROMWICH D H. Instrumented aircraft observations of the katabatic wind regime near Terra Nova Bay[J]. Monthly Weather Review, 1989, 117(7): 1570—1585.

    29 PARISH T R. Numerical simulation of the Terra Nova Bay katabatic wind regime[J]. Antarctic Journal of the US, 1987, 22: 252—254.

    30 BROMWICH D H, PARISH T R, ZORMAN C A. The confluence zone of the intense katabatic winds at Terra Nova Bay, Antarctica, as derived from airborne sastrugi surveys and mesoscale numerical modeling[J]. Journal of Geophysical Research: Atmospheres, 1990,95(D5): 5495—5509.

    31 MORELLI S, PARMIGGIANI F. Wind over Terra Nova Bay (Antarctica) during a polynya event: Eta model simulations and satellite microwave observations[J]. The European Physical Journal Plus, 2013, 128(11): 135.

    32 丁明虎, 卞林根, 張林, 等. 中國(guó)南極維多利亞地新站氣象特征分析[J]. 極地研究, 2015, 27(4): 344—350.

    33 趙杰臣, 程凈凈, 孟上, 等. 南極羅斯海恩克斯堡島短期氣候特征分析[J]. 極地研究, 2015, 27(2): 140—149.

    34 POWERS J G, MANNING K W, BROMWICH D H, et al. A decade of Antarctic science support through AMPS[J]. Bulletin of the American Meteorological Society, 2012, 93(11): 1699—1712.

    35 BROMWICH D H, MONAGHAN A J, MANNING K W, et al. Real-time forecasting for the Antarctic: an evaluation of the Antarctic Mesoscale Prediction System (AMPS)[J]. Monthly Weather Review, 2005, 133(3): 579—603.

    36 孫啟振, 楊清華, 張林. 南極數(shù)值天氣預(yù)報(bào)應(yīng)用與研究進(jìn)展[J]. 極地研究, 2011, 23(2): 128—137.

    37 PARISH T R, BROMWICH D H. Reexamination of the near-surface airflow over the Antarctic continent and implications on atmospheric circulations at high southern latitudes[J]. Monthly Weather Review, 2007, 135(5): 1961—1973.

    38 BROMWICH D H, OTIENO F O, HINES K M, et al. Comprehensive evaluation of polar weather research and forecasting model performance in the Antarctic[J]. Journal of Geophysical Research: Atmospheres, 2013, 118(2): 274—292.

    39 LAZZARA M A, WEIDNER G A, KELLER L M, et al. Antarctic automatic weather station program: 30 years of polar observation[J].Bulletin of the American Meteorological Society, 2012, 93(10): 1519—1537.

    猜你喜歡
    難言諾瓦特拉
    特拉法爾加海戰(zhàn)
    女人的難言之痛
    特拉揚(yáng)漫游大海
    螞蟻與諾瓦泰羅樹(shù)
    螞蟻與諾瓦泰羅樹(shù)
    螞蟻與諾瓦泰羅樹(shù)
    絕經(jīng)前后的“難言之痛”
    女人的難言之痛
    特拉維夫
    世界博覽(2017年2期)2017-04-19 14:29:07
    難言
    777久久人妻少妇嫩草av网站| 久久久久久免费高清国产稀缺| 夜夜夜夜夜久久久久| 99精国产麻豆久久婷婷| 午夜福利,免费看| 国产xxxxx性猛交| 天天躁夜夜躁狠狠躁躁| 一进一出抽搐gif免费好疼 | 黑人猛操日本美女一级片| 老司机福利观看| 777米奇影视久久| 国产男女内射视频| 久久久久国产一级毛片高清牌| 18禁黄网站禁片午夜丰满| 99精品久久久久人妻精品| 久9热在线精品视频| 后天国语完整版免费观看| 天堂中文最新版在线下载| 一夜夜www| 一夜夜www| 丝袜美腿诱惑在线| 精品久久久久久久毛片微露脸| 人人妻人人爽人人添夜夜欢视频| 久久精品国产亚洲av香蕉五月 | 老司机午夜十八禁免费视频| 99精品久久久久人妻精品| 国产av又大| 国产极品粉嫩免费观看在线| 国产黄色免费在线视频| 亚洲五月天丁香| av网站免费在线观看视频| 交换朋友夫妻互换小说| 亚洲国产欧美一区二区综合| 麻豆av在线久日| 久久久久久久久久久久大奶| 日韩高清综合在线| 桃红色精品国产亚洲av| 亚洲最大成人手机在线| 久久性视频一级片| 久久久国产成人免费| 亚洲国产色片| 亚洲人与动物交配视频| 51国产日韩欧美| 成年人黄色毛片网站| 国产亚洲av嫩草精品影院| 男人和女人高潮做爰伦理| 成人无遮挡网站| 性色avwww在线观看| 给我免费播放毛片高清在线观看| 久久午夜亚洲精品久久| 精品欧美国产一区二区三| 男人和女人高潮做爰伦理| 18禁黄网站禁片午夜丰满| 久久久国产成人精品二区| 欧美成人免费av一区二区三区| 欧美中文综合在线视频| 午夜精品一区二区三区免费看| 国产亚洲欧美98| 亚洲精品456在线播放app | 久久久久亚洲av毛片大全| 欧美+亚洲+日韩+国产| 成人国产一区最新在线观看| 久久久久亚洲av毛片大全| 无遮挡黄片免费观看| av在线蜜桃| 欧美日韩瑟瑟在线播放| 久久久精品大字幕| 高潮久久久久久久久久久不卡| 波多野结衣巨乳人妻| 中文资源天堂在线| 久久精品国产亚洲av涩爱 | 国产精品爽爽va在线观看网站| 国产亚洲欧美在线一区二区| 国产欧美日韩精品亚洲av| 90打野战视频偷拍视频| 嫩草影院精品99| 香蕉丝袜av| 俺也久久电影网| 三级毛片av免费| 又爽又黄无遮挡网站| 亚洲第一欧美日韩一区二区三区| 亚洲av电影不卡..在线观看| 国产成人啪精品午夜网站| 脱女人内裤的视频| 好看av亚洲va欧美ⅴa在| 国产探花极品一区二区| 欧美三级亚洲精品| 亚洲专区中文字幕在线| 日韩成人在线观看一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 久久这里只有精品中国| 色吧在线观看| 九色成人免费人妻av| 偷拍熟女少妇极品色| 精品欧美国产一区二区三| www国产在线视频色| 中文在线观看免费www的网站| 亚洲欧美日韩高清专用| 又黄又粗又硬又大视频| 亚洲激情在线av| 丁香欧美五月| 国产精品久久电影中文字幕| 三级男女做爰猛烈吃奶摸视频| 一本综合久久免费| 99久久九九国产精品国产免费| 欧美区成人在线视频| 又粗又爽又猛毛片免费看| 亚洲精品日韩av片在线观看 | 国产黄色小视频在线观看| 欧美成人一区二区免费高清观看| 在线观看日韩欧美| 亚洲成人久久爱视频| 亚洲av电影不卡..在线观看| 人妻丰满熟妇av一区二区三区| 熟妇人妻久久中文字幕3abv| 国产精品久久久久久久电影 | 国产成人av教育| 亚洲av不卡在线观看| 精品久久久久久久人妻蜜臀av| 国产精品乱码一区二三区的特点| 亚洲中文字幕日韩| 日韩欧美精品v在线| 亚洲激情在线av| 90打野战视频偷拍视频| 亚洲aⅴ乱码一区二区在线播放| 久久午夜亚洲精品久久| 国产视频内射| 真实男女啪啪啪动态图| 亚洲国产精品合色在线| 日韩精品青青久久久久久| 国产精品av视频在线免费观看| 国产黄a三级三级三级人| 岛国在线免费视频观看| 99在线视频只有这里精品首页| 老司机午夜福利在线观看视频| 啦啦啦韩国在线观看视频| 欧美激情久久久久久爽电影| 老汉色∧v一级毛片| 3wmmmm亚洲av在线观看| 国产一区二区亚洲精品在线观看| 乱人视频在线观看| 高清毛片免费观看视频网站| 欧美丝袜亚洲另类 | www.色视频.com| 听说在线观看完整版免费高清| 日韩高清综合在线| 国产视频内射| 国产探花在线观看一区二区| 好看av亚洲va欧美ⅴa在| 午夜福利在线观看免费完整高清在 | 少妇人妻精品综合一区二区 | 国产伦精品一区二区三区视频9 | 国产视频内射| 国产在视频线在精品| 午夜免费男女啪啪视频观看 | 欧美午夜高清在线| 天美传媒精品一区二区| 午夜精品在线福利| 99热这里只有精品一区| 欧美乱色亚洲激情| 老汉色∧v一级毛片| 国产精品香港三级国产av潘金莲| 免费高清视频大片| 午夜福利在线观看吧| 亚洲成a人片在线一区二区| 怎么达到女性高潮| 婷婷精品国产亚洲av| 国产男靠女视频免费网站| 黄片小视频在线播放| 免费大片18禁| 久9热在线精品视频| 精品99又大又爽又粗少妇毛片 | 成人午夜高清在线视频| 热99re8久久精品国产| 一本精品99久久精品77| 国产精品99久久久久久久久| av在线蜜桃| 搡女人真爽免费视频火全软件 | 亚洲av电影在线进入| 男插女下体视频免费在线播放| 特大巨黑吊av在线直播| 淫妇啪啪啪对白视频| 久久精品国产亚洲av涩爱 | 亚洲精品乱码久久久v下载方式 | www.999成人在线观看| 91在线精品国自产拍蜜月 | 亚洲精品亚洲一区二区| av欧美777| 舔av片在线| 在线观看日韩欧美| 一本精品99久久精品77| 老熟妇乱子伦视频在线观看| 日本与韩国留学比较| 久久中文看片网| 欧美黑人巨大hd| 久久久国产成人精品二区| 亚洲av中文字字幕乱码综合| 黄色成人免费大全| 色在线成人网| 亚洲内射少妇av| 高清毛片免费观看视频网站| 91麻豆av在线| 亚洲无线观看免费| 97碰自拍视频| 国产精品一区二区三区四区久久| 国模一区二区三区四区视频| 欧美中文日本在线观看视频| 美女被艹到高潮喷水动态| 国产精品1区2区在线观看.| 国产午夜精品久久久久久一区二区三区 | 手机成人av网站| 日本与韩国留学比较| 九色成人免费人妻av| 欧美成人a在线观看| 免费看日本二区| 精品久久久久久久人妻蜜臀av| 久久精品91蜜桃| 97人妻精品一区二区三区麻豆| 噜噜噜噜噜久久久久久91| 久久性视频一级片| 国产成人啪精品午夜网站| 欧美激情久久久久久爽电影| 欧美又色又爽又黄视频| 深夜精品福利| 色精品久久人妻99蜜桃| 99久国产av精品| 亚洲成av人片在线播放无| 12—13女人毛片做爰片一| 亚洲,欧美精品.| 亚洲黑人精品在线| 国产97色在线日韩免费| 可以在线观看毛片的网站| 亚洲国产高清在线一区二区三| www.色视频.com| 母亲3免费完整高清在线观看| 欧美丝袜亚洲另类 | 欧美日韩乱码在线| 美女高潮的动态| 午夜日韩欧美国产| 国产乱人视频| 天堂动漫精品| 久久久久九九精品影院| 欧美另类亚洲清纯唯美| 熟妇人妻久久中文字幕3abv| 人人妻人人澡欧美一区二区| 亚洲av日韩精品久久久久久密| 亚洲最大成人手机在线| 久久久久久久亚洲中文字幕 | 日本黄大片高清| 美女黄网站色视频| 精品一区二区三区人妻视频| 国产主播在线观看一区二区| 一二三四社区在线视频社区8| 性色avwww在线观看| 最近最新中文字幕大全电影3| 欧美区成人在线视频| 亚洲片人在线观看| 欧美zozozo另类| 亚洲精品乱码久久久v下载方式 | 91在线精品国自产拍蜜月 | 美女大奶头视频| 黄色视频,在线免费观看| 免费在线观看成人毛片| 99久久99久久久精品蜜桃| 欧美又色又爽又黄视频| 色综合站精品国产| 青草久久国产| 男女那种视频在线观看| 一区二区三区高清视频在线| 国产精品一及| 日日干狠狠操夜夜爽| 美女高潮的动态| 少妇熟女aⅴ在线视频| 成人国产一区最新在线观看| 国产精品爽爽va在线观看网站| 精品一区二区三区视频在线 | 日韩欧美三级三区| 午夜免费激情av| 一卡2卡三卡四卡精品乱码亚洲| 两性午夜刺激爽爽歪歪视频在线观看| 欧美日本亚洲视频在线播放| 一个人免费在线观看电影| 又爽又黄无遮挡网站| 久久精品亚洲精品国产色婷小说| 嫩草影视91久久| 国产黄a三级三级三级人| tocl精华| 啦啦啦韩国在线观看视频| 日韩人妻高清精品专区| 露出奶头的视频| 久久久久性生活片| 91九色精品人成在线观看| 国产男靠女视频免费网站| or卡值多少钱| 热99re8久久精品国产| 国产亚洲精品av在线| 亚洲乱码一区二区免费版| 淫秽高清视频在线观看| 免费在线观看亚洲国产| 好看av亚洲va欧美ⅴa在| 亚洲精品成人久久久久久| 日本免费a在线| 特大巨黑吊av在线直播| 老熟妇乱子伦视频在线观看| 国产伦精品一区二区三区视频9 | 91麻豆av在线| 在线国产一区二区在线| 亚洲成人久久性| 午夜福利免费观看在线| 99久久精品国产亚洲精品| 特级一级黄色大片| 精品久久久久久,| 久久6这里有精品| 亚洲久久久久久中文字幕| 全区人妻精品视频| 有码 亚洲区| 两人在一起打扑克的视频| 无限看片的www在线观看| 亚洲国产精品sss在线观看| 一夜夜www| 内地一区二区视频在线| 久久国产精品影院| 男女床上黄色一级片免费看| 一级作爱视频免费观看| 精品久久久久久久人妻蜜臀av| 日韩欧美国产一区二区入口| 超碰av人人做人人爽久久 | 亚洲av五月六月丁香网| 一级a爱片免费观看的视频| 母亲3免费完整高清在线观看| 亚洲第一电影网av| 亚洲电影在线观看av| 国产精品三级大全| 97超级碰碰碰精品色视频在线观看| 亚洲国产欧美网| 国产av在哪里看| 男女之事视频高清在线观看| 99热这里只有精品一区| 国产视频一区二区在线看| 亚洲内射少妇av| av片东京热男人的天堂| 天天一区二区日本电影三级| 久久久成人免费电影| 亚洲精品影视一区二区三区av| 国产综合懂色| 亚洲成人久久性| 国产精品自产拍在线观看55亚洲| 亚洲av不卡在线观看| 国产午夜福利久久久久久| 亚洲狠狠婷婷综合久久图片| 免费看日本二区| 精品一区二区三区视频在线观看免费| 久久久久久国产a免费观看| 成年女人毛片免费观看观看9| 黄色日韩在线| 九色成人免费人妻av| 97超级碰碰碰精品色视频在线观看| 色老头精品视频在线观看| 国产高清视频在线观看网站| 一级黄片播放器| 99国产精品一区二区三区| 黄色女人牲交| 亚洲国产色片| 又黄又粗又硬又大视频| 午夜a级毛片| 午夜精品在线福利| 久久精品亚洲精品国产色婷小说| 亚洲国产精品999在线| 午夜福利视频1000在线观看| 脱女人内裤的视频| 亚洲精品亚洲一区二区| 69人妻影院| 国产 一区 欧美 日韩| 亚洲乱码一区二区免费版| 午夜老司机福利剧场| 99精品在免费线老司机午夜| 国产真实伦视频高清在线观看 | 国内精品美女久久久久久| 亚洲av电影在线进入| 精品99又大又爽又粗少妇毛片 | 高清毛片免费观看视频网站| 久久人妻av系列| 日日摸夜夜添夜夜添小说| 国产高清视频在线播放一区| 少妇的逼好多水| 中出人妻视频一区二区| av福利片在线观看| 九色成人免费人妻av| 老司机午夜福利在线观看视频| 黄色成人免费大全| 九色成人免费人妻av| 51国产日韩欧美| av专区在线播放| 国产精品综合久久久久久久免费| 日韩欧美三级三区| 免费看十八禁软件| 精品国内亚洲2022精品成人| 搡老岳熟女国产| 日韩亚洲欧美综合| 成人永久免费在线观看视频| 日本 av在线| 男人的好看免费观看在线视频| 欧美日本视频| 91字幕亚洲| 日韩av在线大香蕉| 欧美大码av| 日韩精品青青久久久久久| 欧美在线一区亚洲| 一级毛片高清免费大全| 人人妻人人澡欧美一区二区| 国产精品久久久久久人妻精品电影| 色尼玛亚洲综合影院| 亚洲真实伦在线观看| 午夜老司机福利剧场| 国产亚洲精品一区二区www| 久久婷婷人人爽人人干人人爱| 99在线视频只有这里精品首页| 伊人久久精品亚洲午夜| 日韩欧美三级三区| 亚洲一区高清亚洲精品| 午夜两性在线视频| 舔av片在线| 国产 一区 欧美 日韩| 亚洲精品亚洲一区二区| 亚洲国产精品成人综合色| 真实男女啪啪啪动态图| 女人高潮潮喷娇喘18禁视频| 啪啪无遮挡十八禁网站| 欧美最黄视频在线播放免费| 在线观看一区二区三区| 国产欧美日韩一区二区精品| 日韩国内少妇激情av| 国产国拍精品亚洲av在线观看 | 久久久久国产精品人妻aⅴ院| 少妇的逼好多水| 国产精品久久电影中文字幕| 亚洲成人久久性| 岛国在线免费视频观看| 99久久精品一区二区三区| 国产不卡一卡二| 18禁黄网站禁片午夜丰满| 国产老妇女一区| 观看美女的网站| 国产探花在线观看一区二区| 九九久久精品国产亚洲av麻豆| 97超级碰碰碰精品色视频在线观看| 亚洲av五月六月丁香网| 亚洲avbb在线观看| 久久久久久久亚洲中文字幕 | 好男人电影高清在线观看| 丝袜美腿在线中文| 18+在线观看网站| 中文字幕人妻丝袜一区二区| 成熟少妇高潮喷水视频| 亚洲精品亚洲一区二区| a在线观看视频网站| 熟妇人妻久久中文字幕3abv| 久久精品影院6| 亚洲男人的天堂狠狠| 婷婷精品国产亚洲av在线| 日韩欧美精品免费久久 | 国内精品一区二区在线观看| 麻豆国产97在线/欧美| 黄色丝袜av网址大全| 国产男靠女视频免费网站| 亚洲最大成人手机在线| 婷婷亚洲欧美| 欧美一区二区亚洲| 国产精品野战在线观看| 色吧在线观看| 婷婷丁香在线五月| 一区二区三区高清视频在线| 久久天躁狠狠躁夜夜2o2o| 国产伦一二天堂av在线观看| 国产一区二区在线av高清观看| 高清在线国产一区| 一级a爱片免费观看的视频| 欧美色欧美亚洲另类二区| 1000部很黄的大片| 精品久久久久久久末码| 9191精品国产免费久久| 国产黄a三级三级三级人| 观看免费一级毛片| 日韩 欧美 亚洲 中文字幕| 亚洲欧美激情综合另类| 国产精品综合久久久久久久免费| 亚洲五月婷婷丁香| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av一区综合| 亚洲成av人片在线播放无| 床上黄色一级片| 久久久久久久亚洲中文字幕 | 高清日韩中文字幕在线| 色av中文字幕| 男女午夜视频在线观看| 欧美日韩一级在线毛片| 国产老妇女一区| 91久久精品国产一区二区成人 | 91av网一区二区| 97超级碰碰碰精品色视频在线观看| 国产爱豆传媒在线观看| 一本久久中文字幕| 在线a可以看的网站| 色播亚洲综合网| 亚洲欧美日韩东京热| 露出奶头的视频| 午夜免费男女啪啪视频观看 | 日本 欧美在线| 特级一级黄色大片| 亚洲性夜色夜夜综合| 精品久久久久久久人妻蜜臀av| 五月玫瑰六月丁香| 丰满人妻熟妇乱又伦精品不卡| 男人舔女人下体高潮全视频| 狂野欧美激情性xxxx| 88av欧美| 国产高清videossex| 国产蜜桃级精品一区二区三区| 色av中文字幕| 丁香欧美五月| 99久久综合精品五月天人人| 村上凉子中文字幕在线| 最新在线观看一区二区三区| 国产乱人视频| 成人av在线播放网站| 身体一侧抽搐| 午夜精品在线福利| 亚洲成人免费电影在线观看| 免费人成视频x8x8入口观看| h日本视频在线播放| 亚洲欧美日韩卡通动漫| 国产免费男女视频| 国产精品女同一区二区软件 | 在线观看美女被高潮喷水网站 | 国产日本99.免费观看| 国产精品久久久久久亚洲av鲁大| 成人无遮挡网站| 小说图片视频综合网站| 天天躁日日操中文字幕| 在线观看舔阴道视频| a级毛片a级免费在线| 18+在线观看网站| 国产精品精品国产色婷婷| 成人精品一区二区免费| 国产一区二区激情短视频| 真人做人爱边吃奶动态| 久久久久久久午夜电影| 免费观看精品视频网站| 日本精品一区二区三区蜜桃| 九色国产91popny在线| 男女做爰动态图高潮gif福利片| 欧美日韩亚洲国产一区二区在线观看| 国产成人系列免费观看| 亚洲欧美精品综合久久99| 搡老岳熟女国产| 成人亚洲精品av一区二区| 国产一区在线观看成人免费| 亚洲第一欧美日韩一区二区三区| 国产午夜精品久久久久久一区二区三区 | 国产99白浆流出| 久久久色成人| 精品电影一区二区在线| 亚洲人与动物交配视频| 精品电影一区二区在线| 51国产日韩欧美| 久久久久亚洲av毛片大全| 在线看三级毛片| 精品久久久久久,| 婷婷丁香在线五月| 国产高清激情床上av| 国产在视频线在精品| 97超视频在线观看视频| 国产aⅴ精品一区二区三区波| 精品乱码久久久久久99久播| 最近最新中文字幕大全电影3| 长腿黑丝高跟| 国产精品电影一区二区三区| 久久精品综合一区二区三区| 精品久久久久久久人妻蜜臀av| 亚洲av电影在线进入| 熟妇人妻久久中文字幕3abv| 狠狠狠狠99中文字幕| 精品福利观看| www.www免费av| 我的老师免费观看完整版| 操出白浆在线播放| 精品人妻1区二区| 国产成人欧美在线观看| 少妇丰满av| 国产精品永久免费网站| 亚洲av美国av| 免费人成在线观看视频色| 久久久久久大精品| 国产伦在线观看视频一区| 久久久久久久精品吃奶| 麻豆成人av在线观看| 国产精品嫩草影院av在线观看 | 岛国视频午夜一区免费看| 啦啦啦韩国在线观看视频| 成人特级av手机在线观看| 中文资源天堂在线| 日日干狠狠操夜夜爽| 网址你懂的国产日韩在线| 日韩高清综合在线| 每晚都被弄得嗷嗷叫到高潮| 免费在线观看成人毛片| 亚洲乱码一区二区免费版| 午夜福利成人在线免费观看| 欧美午夜高清在线| 人妻丰满熟妇av一区二区三区| 久久香蕉国产精品| 亚洲一区二区三区不卡视频| 成人午夜高清在线视频| 在线天堂最新版资源| 亚洲成人久久爱视频|