• <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
    難言
    波多野结衣高清无吗| 人成视频在线观看免费观看| 99香蕉大伊视频| 欧洲精品卡2卡3卡4卡5卡区| 午夜激情av网站| 99久久人妻综合| 日本黄色日本黄色录像| 大香蕉久久成人网| 欧美久久黑人一区二区| 一边摸一边做爽爽视频免费| 国产精品免费一区二区三区在线| 精品免费久久久久久久清纯| 一级片免费观看大全| 免费高清视频大片| 日本五十路高清| 少妇粗大呻吟视频| 国产精品久久久久成人av| 高潮久久久久久久久久久不卡| 国产高清videossex| 国产片内射在线| 黄网站色视频无遮挡免费观看| 欧美黄色淫秽网站| 伊人久久大香线蕉亚洲五| 日韩精品中文字幕看吧| 久久 成人 亚洲| 一级毛片女人18水好多| 一区二区三区国产精品乱码| 老司机深夜福利视频在线观看| 美女国产高潮福利片在线看| 亚洲一区二区三区不卡视频| 又黄又粗又硬又大视频| 午夜福利影视在线免费观看| 一本大道久久a久久精品| 亚洲专区字幕在线| 一级a爱视频在线免费观看| 欧美+亚洲+日韩+国产| av在线播放免费不卡| 日本精品一区二区三区蜜桃| 久久久精品欧美日韩精品| 久久亚洲真实| 别揉我奶头~嗯~啊~动态视频| 少妇的丰满在线观看| 国产日韩一区二区三区精品不卡| 亚洲国产欧美网| 一边摸一边抽搐一进一出视频| 国产野战对白在线观看| 高清黄色对白视频在线免费看| 亚洲成国产人片在线观看| 国产熟女xx| 免费在线观看影片大全网站| 天天影视国产精品| 免费人成视频x8x8入口观看| 国产欧美日韩综合在线一区二区| 午夜免费鲁丝| av天堂在线播放| 黑丝袜美女国产一区| 精品无人区乱码1区二区| 久久久久久久久中文| 日韩免费av在线播放| 99在线人妻在线中文字幕| 久久午夜亚洲精品久久| 中文字幕高清在线视频| 色在线成人网| 久久中文字幕人妻熟女| 国产成人啪精品午夜网站| 激情在线观看视频在线高清| 亚洲,欧美精品.| 精品欧美一区二区三区在线| 50天的宝宝边吃奶边哭怎么回事| 老汉色∧v一级毛片| 午夜福利,免费看| 91国产中文字幕| 丰满迷人的少妇在线观看| 动漫黄色视频在线观看| 亚洲av电影在线进入| x7x7x7水蜜桃| 精品久久久精品久久久| 亚洲第一av免费看| 高清毛片免费观看视频网站 | 亚洲三区欧美一区| 岛国视频午夜一区免费看| 午夜精品在线福利| 色播在线永久视频| 涩涩av久久男人的天堂| 女性生殖器流出的白浆| 极品人妻少妇av视频| 好看av亚洲va欧美ⅴa在| 亚洲成人精品中文字幕电影 | 在线观看免费高清a一片| 亚洲 国产 在线| 嫩草影院精品99| 两人在一起打扑克的视频| 高清欧美精品videossex| 午夜福利欧美成人| 久久久久九九精品影院| 国产亚洲精品综合一区在线观看 | 亚洲欧洲精品一区二区精品久久久| 国产精品99久久99久久久不卡| 91麻豆精品激情在线观看国产 | 国产av在哪里看| 国产亚洲精品久久久久5区| 性欧美人与动物交配| 午夜日韩欧美国产| 精品人妻1区二区| 免费观看人在逋| 熟女少妇亚洲综合色aaa.| 97人妻天天添夜夜摸| 亚洲激情在线av| 成人手机av| 久久青草综合色| 好男人电影高清在线观看| 99国产综合亚洲精品| 狂野欧美激情性xxxx| 欧美精品啪啪一区二区三区| 日韩欧美三级三区| 欧美日韩福利视频一区二区| 国产精品久久久久成人av| avwww免费| 这个男人来自地球电影免费观看| 91精品国产国语对白视频| 99国产极品粉嫩在线观看| 日本免费a在线| 免费久久久久久久精品成人欧美视频| 悠悠久久av| 亚洲人成77777在线视频| 精品欧美一区二区三区在线| 亚洲精品国产区一区二| 免费女性裸体啪啪无遮挡网站| 精品久久久久久成人av| 国产精品综合久久久久久久免费 | 亚洲精华国产精华精| 久久国产精品人妻蜜桃| 国产亚洲精品第一综合不卡| 国产黄色免费在线视频| 亚洲av片天天在线观看| 国产免费av片在线观看野外av| 欧美乱妇无乱码| 美女福利国产在线| 国产不卡一卡二| 十八禁人妻一区二区| 亚洲 欧美一区二区三区| 麻豆成人av在线观看| 国产99白浆流出| 国产精品野战在线观看 | 91av网站免费观看| 天堂影院成人在线观看| 亚洲美女黄片视频| 亚洲色图 男人天堂 中文字幕| 99国产精品99久久久久| 国产蜜桃级精品一区二区三区| 成年人黄色毛片网站| 美女福利国产在线| 精品国产美女av久久久久小说| 真人一进一出gif抽搐免费| 国产成人免费无遮挡视频| 国产亚洲精品第一综合不卡| a在线观看视频网站| 黄色毛片三级朝国网站| 一区二区日韩欧美中文字幕| 99久久国产精品久久久| 国产精品免费视频内射| 叶爱在线成人免费视频播放| 我的亚洲天堂| 亚洲第一av免费看| 视频区图区小说| 99国产精品一区二区三区| 午夜免费观看网址| 少妇的丰满在线观看| 亚洲人成网站在线播放欧美日韩| 亚洲性夜色夜夜综合| 人妻久久中文字幕网| 亚洲 欧美 日韩 在线 免费| 欧美成人性av电影在线观看| 看片在线看免费视频| 国产在线观看jvid| 欧美日韩一级在线毛片| 久久精品亚洲熟妇少妇任你| 在线天堂中文资源库| 桃色一区二区三区在线观看| 热re99久久精品国产66热6| 淫秽高清视频在线观看| 久久影院123| 日韩精品青青久久久久久| av网站免费在线观看视频| 黄色怎么调成土黄色| 国产亚洲精品第一综合不卡| 久久午夜亚洲精品久久| 手机成人av网站| 亚洲男人天堂网一区| 无人区码免费观看不卡| 久久香蕉精品热| 亚洲精品国产精品久久久不卡| 久久久久久亚洲精品国产蜜桃av| 久久精品91蜜桃| 免费高清在线观看日韩| 久久影院123| av中文乱码字幕在线| 人成视频在线观看免费观看| 99国产精品99久久久久| 精品人妻在线不人妻| 桃红色精品国产亚洲av| 久久精品亚洲熟妇少妇任你| 亚洲成av片中文字幕在线观看| 在线av久久热| 久久精品91蜜桃| 国产亚洲精品久久久久5区| 久久热在线av| av欧美777| 别揉我奶头~嗯~啊~动态视频| 淫秽高清视频在线观看| 成人18禁在线播放| 日韩成人在线观看一区二区三区| 久99久视频精品免费| 亚洲一区中文字幕在线| 视频区欧美日本亚洲| 国产精品永久免费网站| 无遮挡黄片免费观看| 欧美成狂野欧美在线观看| 久久精品亚洲精品国产色婷小说| 亚洲熟妇熟女久久| 国产高清videossex| 女人高潮潮喷娇喘18禁视频| 波多野结衣一区麻豆| 中文字幕精品免费在线观看视频| 国产午夜精品久久久久久| 成人三级做爰电影| 大型黄色视频在线免费观看| 香蕉国产在线看| 不卡av一区二区三区| 国产欧美日韩一区二区三| 国产欧美日韩综合在线一区二区| 在线观看www视频免费| 国产伦人伦偷精品视频| 国产成人精品久久二区二区免费| 亚洲熟妇中文字幕五十中出 | 免费av毛片视频| 97超级碰碰碰精品色视频在线观看| 中文亚洲av片在线观看爽| 欧美激情高清一区二区三区| 亚洲欧美日韩高清在线视频| 日本 av在线| 99久久国产精品久久久| 久久婷婷成人综合色麻豆| 亚洲aⅴ乱码一区二区在线播放 | 日韩av在线大香蕉| 咕卡用的链子| 成年人免费黄色播放视频| 国产精品久久电影中文字幕| 免费看a级黄色片| 在线观看免费日韩欧美大片| 国产精品99久久99久久久不卡| 黑丝袜美女国产一区| 黄片播放在线免费| 国产精品国产av在线观看| 性少妇av在线| 老汉色∧v一级毛片| 国产成人精品久久二区二区91| 亚洲精品av麻豆狂野| 一区二区日韩欧美中文字幕| 国产区一区二久久| 一进一出抽搐gif免费好疼 | 亚洲熟女毛片儿| 国产有黄有色有爽视频| 91在线观看av| 露出奶头的视频| 最新美女视频免费是黄的| 久久国产精品男人的天堂亚洲| 在线观看免费午夜福利视频| 老汉色∧v一级毛片| 亚洲精华国产精华精| 一本综合久久免费| 嫩草影院精品99| 精品国产乱子伦一区二区三区| avwww免费| 国产亚洲欧美98| 每晚都被弄得嗷嗷叫到高潮| 久久午夜亚洲精品久久| 日本黄色视频三级网站网址| 中亚洲国语对白在线视频| 久久精品亚洲精品国产色婷小说| 国产精品自产拍在线观看55亚洲| 一级毛片女人18水好多| 99国产精品免费福利视频| 亚洲精品国产一区二区精华液| 男女午夜视频在线观看| 性少妇av在线| 人人妻,人人澡人人爽秒播| 在线观看免费日韩欧美大片| 变态另类成人亚洲欧美熟女 | av网站在线播放免费| 久久久久亚洲av毛片大全| 香蕉国产在线看| 在线观看66精品国产| 18禁美女被吸乳视频| 成人精品一区二区免费| 成人免费观看视频高清| av国产精品久久久久影院| 午夜免费观看网址| 正在播放国产对白刺激| 他把我摸到了高潮在线观看| 久久久久精品国产欧美久久久| 亚洲国产精品999在线| 一级片免费观看大全| 国产日韩一区二区三区精品不卡| 不卡av一区二区三区| 国产不卡一卡二| 夜夜爽天天搞| 久久热在线av| 亚洲 国产 在线| 亚洲国产精品一区二区三区在线| 国产高清videossex| 精品电影一区二区在线| 性色av乱码一区二区三区2| 欧洲精品卡2卡3卡4卡5卡区| 不卡一级毛片| 亚洲性夜色夜夜综合| 高清av免费在线| 亚洲中文字幕日韩| 午夜福利在线免费观看网站| 国产精品 欧美亚洲| 不卡一级毛片| 国产亚洲精品久久久久5区| 午夜免费鲁丝| 欧美日本中文国产一区发布| 日日夜夜操网爽| 国产成年人精品一区二区 | 天堂√8在线中文| 久久热在线av| 天堂动漫精品| 久久青草综合色| 男女床上黄色一级片免费看| 久久人妻av系列| 丰满的人妻完整版| 丰满人妻熟妇乱又伦精品不卡| av在线播放免费不卡| 久久久久国产一级毛片高清牌| 国产亚洲精品第一综合不卡| 成人国语在线视频| 无遮挡黄片免费观看| 操美女的视频在线观看| 亚洲成人精品中文字幕电影 | 999久久久精品免费观看国产| 亚洲一区二区三区欧美精品| 国产亚洲av高清不卡| 亚洲aⅴ乱码一区二区在线播放 | 国产极品粉嫩免费观看在线| 好男人电影高清在线观看| 欧美日韩福利视频一区二区| www.自偷自拍.com| 97超级碰碰碰精品色视频在线观看| 97碰自拍视频| 操出白浆在线播放| 老司机靠b影院| 亚洲人成电影观看| 国产无遮挡羞羞视频在线观看| 超色免费av| 久久香蕉精品热| 十八禁人妻一区二区| 国产精品亚洲一级av第二区| 一区二区三区精品91| 亚洲欧美日韩另类电影网站| 国产成人影院久久av| 久久国产精品影院| 亚洲国产欧美一区二区综合| 国产精品1区2区在线观看.| 老鸭窝网址在线观看| 精品福利观看| 嫁个100分男人电影在线观看| 一边摸一边做爽爽视频免费| 久久午夜亚洲精品久久| 日韩免费av在线播放| 99热国产这里只有精品6| 成人av一区二区三区在线看| 国产欧美日韩综合在线一区二区| av天堂在线播放| 亚洲成a人片在线一区二区| 人人妻,人人澡人人爽秒播| 亚洲人成伊人成综合网2020| 91成人精品电影| 国产一卡二卡三卡精品| 黄片播放在线免费| 国产三级黄色录像| 亚洲午夜理论影院| 日本免费a在线| 老司机亚洲免费影院| 午夜福利免费观看在线| 欧美乱色亚洲激情| 亚洲黑人精品在线| 老司机在亚洲福利影院| 亚洲一区二区三区欧美精品| 不卡一级毛片| 国产成人一区二区三区免费视频网站| 亚洲美女黄片视频| 免费久久久久久久精品成人欧美视频| 亚洲激情在线av| 国产精品一区二区免费欧美| 日本黄色日本黄色录像| 亚洲av片天天在线观看| 国产精品av久久久久免费| 好看av亚洲va欧美ⅴa在| 日本免费a在线| 在线十欧美十亚洲十日本专区| 亚洲七黄色美女视频| 美女国产高潮福利片在线看| 国产成人影院久久av| 久久久久久免费高清国产稀缺| 欧美色视频一区免费| 老熟妇仑乱视频hdxx| 国产成人av教育| 纯流量卡能插随身wifi吗| 午夜老司机福利片| 久久久国产欧美日韩av| 午夜成年电影在线免费观看| 国产成年人精品一区二区 | 日本a在线网址| 色老头精品视频在线观看| 最近最新免费中文字幕在线| 不卡一级毛片| 欧美日韩视频精品一区| 亚洲 国产 在线| 两人在一起打扑克的视频| 丁香欧美五月| 成人18禁在线播放| 91精品三级在线观看| 日日爽夜夜爽网站| 国产av精品麻豆| 757午夜福利合集在线观看| 一级片免费观看大全| 黄网站色视频无遮挡免费观看| 久9热在线精品视频| 波多野结衣高清无吗| 亚洲成国产人片在线观看| 一区福利在线观看| 丰满迷人的少妇在线观看| 啦啦啦 在线观看视频| netflix在线观看网站| 黄频高清免费视频| 亚洲精品成人av观看孕妇| 国产av又大| 国产深夜福利视频在线观看| 欧美一级毛片孕妇| 他把我摸到了高潮在线观看| 啦啦啦在线免费观看视频4| 国产精品爽爽va在线观看网站 | 欧美成人免费av一区二区三区| av欧美777| 欧美久久黑人一区二区| 成人亚洲精品一区在线观看| 在线看a的网站| 超碰成人久久| cao死你这个sao货| 一区二区三区激情视频| 日韩有码中文字幕| 亚洲成人免费av在线播放| 国产亚洲欧美98| 男女下面进入的视频免费午夜 | 欧美成人性av电影在线观看| 亚洲欧美激情在线| 热99re8久久精品国产| 在线观看一区二区三区激情| 黄色丝袜av网址大全| 欧美日本中文国产一区发布| 无限看片的www在线观看| 久久99一区二区三区| 国产亚洲精品久久久久久毛片| 成人黄色视频免费在线看| 成人av一区二区三区在线看| 黄色片一级片一级黄色片| 欧美日本亚洲视频在线播放| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品美女久久久久99蜜臀| 黄色怎么调成土黄色| 嫩草影院精品99| 欧美不卡视频在线免费观看 | 97碰自拍视频| 女人被狂操c到高潮| 少妇被粗大的猛进出69影院| 脱女人内裤的视频| 欧美日韩乱码在线| 日韩中文字幕欧美一区二区| 极品人妻少妇av视频| av欧美777| 天堂影院成人在线观看| 国产精品国产av在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产毛片av蜜桃av| 老司机午夜十八禁免费视频| 在线播放国产精品三级| 日韩精品免费视频一区二区三区| 久久久久国产一级毛片高清牌| 亚洲成人免费av在线播放| 日韩一卡2卡3卡4卡2021年| 91字幕亚洲| 一区在线观看完整版| 超碰97精品在线观看| 精品福利永久在线观看| 欧美 亚洲 国产 日韩一| 亚洲九九香蕉| 看片在线看免费视频| 老汉色av国产亚洲站长工具| 久久久久久久久免费视频了| 在线永久观看黄色视频| 午夜成年电影在线免费观看| 男女下面进入的视频免费午夜 | 涩涩av久久男人的天堂| 欧美黄色淫秽网站| av天堂在线播放| av片东京热男人的天堂| 人人妻,人人澡人人爽秒播| 国产精品久久电影中文字幕| 亚洲精品国产精品久久久不卡| 在线十欧美十亚洲十日本专区| 久久性视频一级片| 国产野战对白在线观看| 精品第一国产精品| 免费观看人在逋| 夜夜看夜夜爽夜夜摸 | 国产高清视频在线播放一区| av网站在线播放免费| 老汉色av国产亚洲站长工具| 欧美日韩瑟瑟在线播放| 国产成人影院久久av| 日本五十路高清| 国产人伦9x9x在线观看| 国产av精品麻豆| 欧美成人午夜精品| 中文字幕人妻熟女乱码| 18禁观看日本| aaaaa片日本免费| 美女午夜性视频免费| 国内久久婷婷六月综合欲色啪| 国产蜜桃级精品一区二区三区| 精品第一国产精品| 人人妻人人爽人人添夜夜欢视频| 9191精品国产免费久久| 日韩大码丰满熟妇| 免费看a级黄色片| 国产aⅴ精品一区二区三区波| 极品教师在线免费播放| 一区二区三区激情视频| www国产在线视频色| 中文字幕人妻丝袜一区二区| 免费日韩欧美在线观看| av福利片在线| 欧美在线黄色| 国产1区2区3区精品| 亚洲精华国产精华精| 日韩精品青青久久久久久| 日韩欧美国产在线观看| 听说在线观看完整版免费高清| 亚洲专区国产一区二区| 久久99热6这里只有精品| 在线免费观看的www视频| 欧美日韩国产亚洲二区| 欧美成狂野欧美在线观看| 精华霜和精华液先用哪个| 国产精品电影一区二区三区| 老鸭窝网址在线观看| 欧美又色又爽又黄视频| 久久国产精品人妻蜜桃| 69人妻影院| 日本免费一区二区三区高清不卡| 内地一区二区视频在线| 久久久精品欧美日韩精品| 在线天堂最新版资源| 女人被狂操c到高潮| 久久人人爽人人爽人人片va | 极品教师在线免费播放| 永久网站在线| 精品99又大又爽又粗少妇毛片 | 色视频www国产| 国产高清视频在线观看网站| 少妇高潮的动态图| 少妇人妻精品综合一区二区 | www.熟女人妻精品国产| 俄罗斯特黄特色一大片| 国内精品美女久久久久久| 亚洲美女黄片视频| 日韩大尺度精品在线看网址| 女人十人毛片免费观看3o分钟| 91午夜精品亚洲一区二区三区 | 日本黄大片高清| 精品久久久久久成人av| 九色成人免费人妻av| av福利片在线观看| 欧美3d第一页| 国产成人福利小说| 精品不卡国产一区二区三区| 一级av片app| 国产亚洲欧美在线一区二区| 国产精品亚洲一级av第二区| 老司机午夜福利在线观看视频| 长腿黑丝高跟| 亚洲专区国产一区二区| 又爽又黄无遮挡网站| 他把我摸到了高潮在线观看| 国产精品电影一区二区三区| 在线观看美女被高潮喷水网站 | 国产精品美女特级片免费视频播放器| 亚洲精品在线观看二区| 亚洲avbb在线观看| 天堂√8在线中文| 老司机深夜福利视频在线观看| 毛片一级片免费看久久久久 | 9191精品国产免费久久| 91在线精品国自产拍蜜月| h日本视频在线播放| 给我免费播放毛片高清在线观看| 国产乱人视频| 蜜桃亚洲精品一区二区三区| 久久这里只有精品中国| 观看免费一级毛片| 国产亚洲欧美在线一区二区| 欧美绝顶高潮抽搐喷水| 国产成人影院久久av|