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

    南極中山站風(fēng)的長期變化特征及其可能原因

    2022-10-15 04:15:34武維剛張文千武麥鳳馬靖凱凌新鋒
    極地研究 2022年3期
    關(guān)鍵詞:中山站緯向日數(shù)

    武維剛 張文千 武麥鳳 馬靖凱 凌新鋒

    (1華山氣象站, 陜西 華陰 714200;2中國氣象科學(xué)研究院, 北京 100081;3渭南市氣象局, 陜西 渭南 714000;4國家海洋環(huán)境預(yù)報(bào)中心, 北京 100081;5壽縣氣象局, 安徽 壽縣 232200)

    0 引言

    南極的氣候變化在全球氣候變化中起著重要作用。在全球變暖背景下, 很多研究都集中在分析和模擬南極地區(qū)溫度的變化趨勢。由于觀測資料缺乏, 對(duì)其他氣象要素的氣候變化研究較少。風(fēng)是研究大氣動(dòng)力學(xué)和氣候變化的一個(gè)重要參量,地面風(fēng)也是氣候?qū)W中最重要的研究對(duì)象之一。南極中山站是中國在南極地區(qū)建立的第二個(gè)常年科考站, 1989年2月建成于東南極拉斯曼丘陵、普里茲灣東南沿岸。胡勝利[1]指出中山站天氣惡劣,具有大風(fēng)多、風(fēng)速大和大風(fēng)時(shí)間長的特點(diǎn), 并受到多種天氣形勢影響。大風(fēng)天氣同時(shí)會(huì)造成吹雪和雪暴天氣[2]。大風(fēng)及其引發(fā)的次生災(zāi)害會(huì)對(duì)南極地區(qū)的科考站建筑造成破壞, 危及科研人員的人身安全, 影響航空飛行器的起降, 不利于船只航行。

    前人研究[3-6]指出中山站常年盛行東風(fēng)及偏東風(fēng); 平均風(fēng)速約為7 m·s-1, 具有冬強(qiáng)夏弱的年變化特征, 且風(fēng)速日變化在夏季明顯, 冬季不明顯; 年平均大風(fēng)日數(shù)為159 d。在全球變暖背景下,盡管中山站1989—2008年氣溫?zé)o明顯變化, 但風(fēng)速呈減小趨勢[4,6]。Turner等[7]對(duì)南極地區(qū)19個(gè)氣象站50年的溫度、氣壓和風(fēng)速資料進(jìn)行分析,指出大部分考察站風(fēng)速變化與南半球環(huán)狀模(Southern Hemisphere Annular Mode, SAM)的變化一致, 當(dāng)SAM指數(shù)為正位相時(shí), 環(huán)南極西風(fēng)增強(qiáng),使各站觀測到的風(fēng)速增加[8]。但受到地理位置與局地天氣系統(tǒng)影響, 相距不遠(yuǎn)的考察站之間風(fēng)速仍存在不同變化。例如, 1958—2000年, 戴維斯站風(fēng)速呈增加趨勢, 但和平站風(fēng)速呈減小趨勢[7]。1989—2008年, 戴維斯站風(fēng)速呈增加趨勢, 但中山站風(fēng)速呈減小趨勢[6,9]。

    南極中山站所處的伊麗莎白公主地是非常典型的冰蓋下降風(fēng)頻發(fā)地區(qū)。南極中山站下降風(fēng)風(fēng)速一般為5~6級(jí), 最大可達(dá)7~8級(jí), 風(fēng)向?yàn)槠珫|風(fēng)[6,9-12]。普里茲灣沿岸的下降風(fēng)強(qiáng)度受到南極大陸冷高壓強(qiáng)弱以及位置的影響[10]。此外, 受下降風(fēng)影響的強(qiáng)弱還與距冰蓋斜坡的遠(yuǎn)近有關(guān)。與中山站相距不遠(yuǎn)的戴維斯站, 由于距冰蓋斜坡較遠(yuǎn),受下降風(fēng)的影響較中山站要小[13]。

    本文利用1989—2020年中山站地面風(fēng)等氣象觀測數(shù)據(jù), 對(duì)中山站風(fēng)速、大風(fēng)日數(shù)以及下降風(fēng)的氣候變化特征進(jìn)行分析, 同時(shí)結(jié)合NCEP/NCAR再分析數(shù)據(jù)探討上述變化的成因,以期了解南極自由大氣底部的氣候變化并為深入理解區(qū)域氣候變化提供參考意義。

    1 資料和方法

    本文所用的資料均來自南極中山氣象臺(tái), 其地面氣象觀測儀器的安裝、運(yùn)行和觀測方法均按照中國氣象局頒布的《地面氣象觀測規(guī)范》執(zhí)行[14], 其觀測數(shù)據(jù)經(jīng)過嚴(yán)格的質(zhì)量審核。資料時(shí)間跨度為1989年12月—2020年12月, 資料類型包括逐年、月、日和小時(shí)的風(fēng)、氣溫和氣壓等地面氣象觀測資料。同時(shí), 為了分析大氣環(huán)流對(duì)風(fēng)速的影響, 本文還選取了1990—2020年NCEP/NCAR再分析資料, 包括逐月海平面氣壓場、表面氣溫場、高度場和風(fēng)場等資料, 水平分辨率為2.5°× 2.5°[15]。文中季節(jié)按12月至次年2月為夏季、3—5月為秋季、6—8月為冬季、9—11月為春季劃分。

    《地面氣象觀測規(guī)范》中規(guī)定, 瞬時(shí)風(fēng)速達(dá)到或超過17 m·s-1(或目測估計(jì)風(fēng)力達(dá)到或超過8級(jí))的風(fēng), 記為大風(fēng); 若某日中有大風(fēng)出現(xiàn), 一般稱該日為大風(fēng)日; 一日中出現(xiàn)多次大風(fēng)記為一個(gè)大風(fēng)日[14]。根據(jù)丁卓銘等[10]分析的普里茲灣沿岸下降風(fēng)特征, 本文將晴朗少云天氣風(fēng)速超過5級(jí)(≥8 m·s-1), 風(fēng)向介于東北至東南(33.75°~146.25°)的風(fēng)看作下降風(fēng)。

    將統(tǒng)計(jì)時(shí)段內(nèi)風(fēng)速按照下式進(jìn)行分解:

    其中,F為風(fēng)速, π為圓周率,φ為風(fēng)向。分解后u為東-西向風(fēng)速,v為南-北向風(fēng)速, 按照正、負(fù)值區(qū)分為4個(gè)分量, 然后各分量取絕對(duì)值后求和, 再除以總觀測次數(shù), 即得到相應(yīng)時(shí)段內(nèi)經(jīng)、緯向及東、西、南、北風(fēng)分量的平均風(fēng)速。計(jì)算年、季平均風(fēng)速傾向率時(shí)采用一元線性趨勢擬合, 可得到風(fēng)速線性變化率b和相關(guān)系數(shù)r[16]。文中對(duì)氣候變化趨勢的性質(zhì)和幅度采用氣候傾向率(b×10)、趨勢系數(shù)(無量綱)來表示,趨勢檢驗(yàn)采用t檢驗(yàn)方法[17]。采用 Mann-Kendall檢驗(yàn)法(以下簡稱M-K檢驗(yàn))對(duì)時(shí)間序列進(jìn)行突變分析[18]。

    2 中山站風(fēng)的變化特征

    2.1 風(fēng)的年際變化

    圖1a可以看出,1990—2020年中山站年平均風(fēng)速最大值為7.8 m·s-1, 出現(xiàn)在1991年和2007年; 最小值為5.7 m·s-1, 出現(xiàn)在2019年; 年平均風(fēng)速為6.9 m·s-1。圖2 中UF為按時(shí)間序列順序計(jì)算出的統(tǒng)計(jì)量序列, UB為按時(shí)間序列逆序計(jì)算出的統(tǒng)計(jì)量序列[18]。通過M-K檢驗(yàn)發(fā)現(xiàn), UF和UB兩條曲線相交, 交點(diǎn)位于臨界線內(nèi), 根據(jù)交點(diǎn)位置, 確定風(fēng)速在2002年附近發(fā)生突變。1990—2001年平均風(fēng)速為7.3 m·s-1, 2002—2020年平均風(fēng)速為6.6 m·s-1, 這說明中山站年平均風(fēng)速減小。通過計(jì)算線性趨勢也證實(shí), 年平均風(fēng)速有減少趨勢,傾向率為-0.4 m·s-1·(10 a)-1(通過0.001顯著性水平檢驗(yàn))。M-K檢驗(yàn)也得到一致的結(jié)果, 1992年以后年平均風(fēng)速呈下降趨勢(UF<0), 2003—2020年風(fēng)速下降趨勢顯著, UF值超過α=0.05臨界線甚至超過0.001顯著性水平(u0.001=2.56)臨界線。從圖1還可以看出, 中山站年平均氣溫和氣壓變化趨勢并不明顯, 氣候傾向率分別為0.015℃·(10 a)-1、-0.037 hPa·(10 a)-1, 均未通過0.05顯著性水平檢驗(yàn)。

    圖1 1990—2020年中山站風(fēng)年際變化。a)年平均風(fēng)速、氣溫;b)年平均氣壓Fig.1.Interannual variation of wind at Zhongshan Station from 1990 to 2020.a) annual mean wind speed and temperature;b) annual mean pressure

    圖 2 1990—2020年中山站年平均風(fēng)速時(shí)間序列Mann-Kendall突變檢驗(yàn)(虛線為α = 0.05顯著性水平臨界值)Fig.2.Mann-Kendall test of the time series of annual mean wind speed at Zhongshan Station from 1990 to 2020(dotted line is α = 0.05 critical value of significance level)

    中山站位于環(huán)極低壓帶的南部和南極大陸地面冷高壓北部之間, 常年盛行偏東氣流。年最多風(fēng)向以東北至東東南為主, 東風(fēng)出現(xiàn)頻率最高,1990—2020年間有22年最多風(fēng)向?yàn)闁|風(fēng)。

    2.2 風(fēng)的季節(jié)變化

    1990—2020年中山站月平均風(fēng)速隨時(shí)間呈倒“U”字型(圖3a), 中間高, 兩頭低。月平均風(fēng)速為6.9 m·s-1, 1月平均風(fēng)速最小為5.0 m·s-1,2月風(fēng)速迅速增大, 6月和8月風(fēng)速達(dá)到全年最大為7.9 m·s-1, 9月風(fēng)速開始減小, 3—9月平均風(fēng)速變化緩慢且高于平均值。中山站風(fēng)速存在明顯的季節(jié)差異, 冬季風(fēng)速最大, 秋季次之,春季較小, 夏季最小(圖3)。1990—2020年中山站春、夏、秋、冬季平均風(fēng)速分別為6.7、5.7、7.4、7.7 m·s-1。春季平均風(fēng)速最大值為8.3 m·s-1,出現(xiàn)在1994年; 最小值為5.7 m·s-1, 出現(xiàn)在2016年。秋季和冬季平均風(fēng)速最大值分別為8.9、9.9 m·s-1, 均出現(xiàn)在2018年, 最小值分別為5.8、5.6 m·s-1, 均出現(xiàn)在2019年。夏季平均風(fēng)速最大值為6.8 m·s-1, 出現(xiàn)在1992年; 最小值為4.4 m·s-1, 出現(xiàn)在2020年。中山站全年風(fēng)向變化較小, 2—11月東風(fēng)出現(xiàn)頻率最高, 東東北風(fēng)次之。1月和12月東東北風(fēng)出現(xiàn)頻率最高,東風(fēng)次之。中山站地面風(fēng)受南極大陸冷高壓和繞極低壓帶控制[19]。由于夏季南極大陸冷高壓和繞極低壓帶強(qiáng)度較弱, 且繞極低壓帶的平均位置明顯北移, 使得高、低壓之間氣壓梯度小,從而導(dǎo)致夏季中山站平均風(fēng)速最小。當(dāng)夏季轉(zhuǎn)向冬季時(shí), 南極大陸冷高壓逐漸加強(qiáng), 同時(shí)繞極低壓帶的平均位置南移, 使得高、低壓中心的氣壓梯度增大, 形成強(qiáng)烈且穩(wěn)定的東風(fēng), 從而使冬季中山站平均風(fēng)速最大。

    圖3 1990—2020年中山站不同季節(jié)的平均風(fēng)速變化。a)月平均風(fēng)速變化;b)春季;c)夏季;d)秋季;e)冬季Fig.3.Interannual variation of mean wind speed at Zhongshan Station from 1990 to 2020.a) monthly mean wind speed;b) spring; c) summer; d) autumn; e) winter

    中山站各季節(jié)平均風(fēng)速的變化趨勢(圖3)與年平均風(fēng)速變化相同, 均呈減小趨勢。對(duì)比各季節(jié)平均風(fēng)速變化趨勢系數(shù)發(fā)現(xiàn), 夏季和冬季減小趨勢最為顯著, 氣候傾向率、趨勢系數(shù)分別為-0.53 m·s-1·(10 a)-1、-0.09和-0.48 m·s-1·(10 a)-1、-0.05。春季和秋季風(fēng)速氣候傾向率、趨勢系數(shù)分別為-0.25 m·s-1·(10 a)-1、-0.04和-0.31 m·s-1·(10 a)-1、-0.04。顯著性檢驗(yàn)表明, 除秋季外其他季節(jié)均通過0.05顯著性水平檢驗(yàn), 表明秋季風(fēng)速變化趨勢不明顯。

    2.3 風(fēng)的日變化

    圖4為1990—2020年中山站年和四季定時(shí)10分鐘風(fēng)速的日變化??梢钥闯? 中山站風(fēng)速日變化呈單峰單谷型, 風(fēng)速在凌晨5時(shí)達(dá)到最大,峰值為8.0 m·s-1, 之后風(fēng)速開始減小, 下午15時(shí)風(fēng)速降至最小, 谷值為5.6 m·s-1, 隨后風(fēng)速又開始增加。日平均風(fēng)速為6.8 m·s-1, 11—21時(shí)風(fēng)速均小于日平均風(fēng)速。中山站定時(shí)平均風(fēng)速和氣溫呈負(fù)相關(guān),r為-0.97, 并通過0.001顯著性水平檢驗(yàn)。從季節(jié)變化來看, 春季和夏季受下降風(fēng)影響風(fēng)速具有明顯的日變化特征, 秋季和冬季風(fēng)速日變化不明顯。這是因?yàn)榍锛竞投景讜儠r(shí)間短,日照時(shí)間很少, 普利茲灣海面以及拉斯曼丘陵絕大部分被冰雪覆蓋, 整個(gè)地區(qū)下墊面單一, 受太陽輻射加熱比較平衡, 晝夜變化小, 因此風(fēng)速的日變化不明顯。春季和夏季白天受太陽輻射加熱的時(shí)間長, 容易產(chǎn)生海水與冰面的受熱不平衡, 夜間不同下墊面輻射散熱的差異大, 使海陸溫差加大。中山站位于岸邊緊鄰冰蓋, 海平面與冰蓋存在很大的高度差, 導(dǎo)致海風(fēng)和冰蓋下降風(fēng)的增強(qiáng)[3]。春季風(fēng)速在4時(shí)達(dá)到最大, 峰值為8.6 m·s-1; 15時(shí)風(fēng)速最小, 谷值為4.6 m·s-1。與春季不同, 夏季風(fēng)速在5時(shí)達(dá)到最大, 峰值為8.0 m·s-1; 16時(shí)風(fēng)速最小, 谷值為3.5 m·s-1, 夏季風(fēng)速峰值和谷值的出現(xiàn)時(shí)間均比春季晚1小時(shí)。總體來看, 中山站全年和四季晚上平均風(fēng)速大于白天。

    圖4 1990—2020年中山站年和四季定時(shí)風(fēng)速的日變化Fig.4.Diurnal variation of annual and seasonal wind speed at Zhongshan Station from 1990 to 2020

    2.4 經(jīng)、緯向風(fēng)分量的變化特征

    圖5給出了1990—2020年中山站經(jīng)、緯向平均風(fēng)速的年際變化。由圖可見, 中山站緯向風(fēng)速明顯大于經(jīng)向風(fēng)速。與年平均風(fēng)速變化趨勢相同, 緯向年平均風(fēng)速呈減小趨勢, 而經(jīng)向風(fēng)速變化趨勢不明顯。緯向平均風(fēng)速在1991年達(dá)到最大值7.2 m·s-1, 與年平均風(fēng)速最大值出現(xiàn)時(shí)間相同。2002—2020年緯向平均風(fēng)速為5.7 m·s-1, 較1990—2001年緯向平均風(fēng)速減小了1.0 m·s-1(表1), 2005年和2006年緯向年平均風(fēng)速達(dá)到最小值4.0 m·s-1。與緯向風(fēng)速變化相反, 經(jīng)向平均風(fēng)速在2002年以后明顯增大, 經(jīng)向平均風(fēng)速由1.8 m·s-1增加至2.6 m·s-1, 并在2007年達(dá)到最大值5.8 m·s-1。經(jīng)向平均風(fēng)速最小值為1.3 m·s-1, 出現(xiàn)在2019年, 與年平均風(fēng)速最小值出現(xiàn)時(shí)間相同。

    圖5 1990—2020年中山站經(jīng)、緯向平均風(fēng)速的年際變化Fig.5.Interannual variation of mean wind speed in meridional and zonal at Zhongshan Station from 1990 to 2020

    表1 2002年前后中山站經(jīng)、緯向及東、南、西、北風(fēng)分量平均風(fēng)速Table 1.Mean wind speed in meridional, zonal and components of east, south, west and north wind at Zhongshan Station around 2002單位: m·s-1

    對(duì)經(jīng)向和緯向風(fēng)進(jìn)行進(jìn)一步分解, 圖6為1990—2020年中山站東、西、南、北風(fēng)分量平均風(fēng)速年際變化。從圖中可以看出, 東風(fēng)分量平均風(fēng)速明顯大于其他方向風(fēng)分量, 東風(fēng)分量年際變化與緯向風(fēng)速年際變化基本一致。2002—2020年東風(fēng)分量平均風(fēng)速較1990—2001年減小了1.0 m·s-1。東風(fēng)分量平均風(fēng)速在1991年最大為7.7 m·s-1, 2006年最小為4.4 m·s-1。2002—2020年西風(fēng)分量平均風(fēng)速較1990—2001年減小了0.3 m·s-1, 年平均風(fēng)速最大值為2.4 m·s-1, 最小值為1.1 m·s-1, 極值出現(xiàn)時(shí)間與東風(fēng)分量極值出現(xiàn)時(shí)間相同。南風(fēng)分量平均風(fēng)速在1991年最大為3.6 m·s-1, 2010年最小為1.2 m·s-1。2002—2020年南風(fēng)分量平均風(fēng)速較1990—2001年減小了1.0 m·s-1。北風(fēng)分量平均風(fēng)速1994年最小為1.0 m·s-1, 之后風(fēng)速開始增大, 與經(jīng)向年平均風(fēng)速變化基本一致。2002—2020年北風(fēng)分量平均風(fēng)速較1990—2001年增加了1.4 m·s-1, 并在2007年風(fēng)速達(dá)到最大為6.2 m·s-1。

    圖6 1990—2020年中山站東、西、南、北風(fēng)分量平均風(fēng)速的年際變化Fig.6.Interannual variation of the mean wind speed of east, west, south and north wind components at Zhongshan Station from 1990 to 2020

    表2為中山站年與四季經(jīng)、緯向及東、南、西、北風(fēng)分量平均風(fēng)速氣候傾向率及趨勢系數(shù)。從經(jīng)、緯向風(fēng)速的氣候傾向率和趨勢系數(shù)來看,緯向風(fēng)明顯大于經(jīng)向風(fēng), 緯向風(fēng)的全年與四季風(fēng)速氣候傾向率和趨勢系數(shù)均為負(fù)值, 說明緯向風(fēng)速呈減小趨勢。緯向年平均風(fēng)速氣候傾向率為-0.34 m·s-1·(10 a)-1(通過0.05顯著性水平檢驗(yàn)),四季中只有夏季通過0.01顯著性水平檢驗(yàn), 氣候傾向率為-0.46 m·s-1·(10 a)-1。經(jīng)向風(fēng)的全年和四季風(fēng)速變化趨勢不明顯, 均未通過顯著性檢驗(yàn)。從風(fēng)分量氣候傾向率和趨勢系數(shù)來看, 南風(fēng)分量呈減小趨勢且變化速率遠(yuǎn)大于其他風(fēng)分量。南風(fēng)分量在夏季減小速率最大, 氣候傾向率為-0.65 m·s-1·(10 a)-1; 在春季減小速率最小, 氣候傾向率為-0.46 m·s-1·(10 a)-1, 均通過0.01顯著性水平檢驗(yàn)。東風(fēng)分量夏季氣候傾向率和趨勢系數(shù)明顯大于其他季節(jié)且與其變化趨勢相反, 風(fēng)速呈減小趨勢, 氣候傾向率為-0.48 m·s-1·(10 a)-1, 通過0.01顯著性水平檢驗(yàn)。西風(fēng)和北風(fēng)分量變化趨勢不顯著, 其全年和四季氣候傾向率和趨勢系數(shù)均未通過顯著性水平檢驗(yàn)。

    表2 中山站年與四季經(jīng)、緯向及東、南、西、北風(fēng)分量平均風(fēng)速氣候傾向率[m·s-1·(10 a)-1]/趨勢系數(shù)Table 2.The annual and seasonal climatic tendency rate [m·s-1·(10 a)-1]/trend coefficient of the mean wind speed in meridional, zonal and east, south, west and north components of Zhongshan Station

    以上分析表明中山站緯向平均風(fēng)速呈顯著減小趨勢, 經(jīng)向平均風(fēng)速變化趨勢不顯著; 南風(fēng)分量年平均風(fēng)速氣候傾向率和趨勢系數(shù)都顯著大于其他風(fēng)分量, 風(fēng)速呈顯著減小趨勢; 西風(fēng)和北風(fēng)分量變化趨勢不顯著; 東風(fēng)分量在夏季呈減小趨勢, 其他季節(jié)變化趨勢不顯著。

    3 下降風(fēng)的氣候變化特征

    圖7展示了1990—2020年中山站年和四季出現(xiàn)下降風(fēng)頻數(shù)的年際變化。由圖可見, 中山站年下降風(fēng)頻數(shù)年際變化明顯, 平均為1757.2次, 最大值為2211次(1994年), 最小值為1087次(2020年)。年下降風(fēng)頻數(shù)呈減少趨勢,氣候傾向率為-119.8次·(10 a)-1(通過0.05顯著性水平檢驗(yàn))。與風(fēng)速季節(jié)變化相同, 中山站冬季出現(xiàn)下降風(fēng)次數(shù)最多, 秋季次之, 春季較少, 夏季最少。冬、秋、春、夏季年平均下降風(fēng)頻數(shù)分別為534.6次、481.2次、318.5次和217.7次。冬季年平均下降風(fēng)頻數(shù)最大值為837次(2006年), 最小值為222次(2010年)。秋季最大值為659次(2007年), 最小值為255次(1999年)。春季最大值為442次(1991年), 最小值為154次(2010年)。夏季最大值為373次(2004年), 最小值為96次(2020年)。

    圖7 1990—2020年中山站年、夏季和冬季下降風(fēng)頻數(shù)的年際變化Fig.7.Interannual variation of annual, summer and winter katabatic wind frequency at Zhongshan Station from 1990 to 2020

    1990—2020年中山站年和四季下降風(fēng)頻數(shù)均呈減少趨勢, 其中夏季和冬季下降風(fēng)頻數(shù)減少趨勢最為顯著, 氣候傾向率、趨勢系數(shù)分別為-30.3次·(10 a)-1、-0.05和-66.0次·(10 a)-1、-0.04,均通過0.05顯著性水平檢驗(yàn)。春季和秋季下降風(fēng)頻數(shù)氣候傾向率、趨勢系數(shù)分別為-0.7次·(10 a)-1、0.0和-19.6次·(10 a)-1、-0.02, 未通過0.05顯著性水平檢驗(yàn), 表明春季和秋季下降風(fēng)變化趨勢不顯著。

    4 大風(fēng)日數(shù)的變化特征

    中山站大風(fēng)主要受極地氣旋、極地大陸冷高壓和東移鋒面擾動(dòng)的系統(tǒng)影響[1]。圖8可以看出,中山站年大風(fēng)日數(shù)范圍為81~196 d, 1999年出現(xiàn)的大風(fēng)日數(shù)最多, 2016年和2019年最少。年平均大風(fēng)日數(shù)為140 d, 占全年總?cè)諗?shù)的38%。年大風(fēng)日數(shù)在100 d以下的有3 a, 分別是2010年、2016年和2019年; 100~200 d有28 a, 占總年份的90%。1990—2002年期間只有1993年的大風(fēng)日數(shù)小于平均值, 2003—2020年期間僅有4 a(2004年、2006年、2007年和2018年)大風(fēng)日數(shù)大于平均值, 初步表明年大風(fēng)日數(shù)有減少趨勢。對(duì)1990—2020年大風(fēng)日數(shù)進(jìn)行趨勢分析, 結(jié)果表明雖然中山站大風(fēng)日數(shù)高(低)值年與平均風(fēng)速高(低)值年并不一一對(duì)應(yīng), 但是大風(fēng)日數(shù)的線性趨勢與平均風(fēng)速的變化趨勢一致, 均呈減小趨勢, 年大風(fēng)日數(shù)氣候傾向率為-25.8 d·(10 a)-1(通過0.001顯著性水平檢驗(yàn))。大風(fēng)日數(shù)與年平均風(fēng)速呈正相關(guān),r為0.86, 并通過0.001顯著性水平檢驗(yàn)。中山站大風(fēng)風(fēng)向主要為東風(fēng),頻率為42%, 東東南、東東北大風(fēng)次之, 頻率分別為27%和16%。

    圖8 1990—2020年中山站大風(fēng)日數(shù)的年際變化Fig.8.Interannual variation of gale days at Zhongshan Station from 1990 to 2020

    中山站月大風(fēng)日數(shù)隨時(shí)間呈“凸”字型, 月平均大風(fēng)日數(shù)為11.7 d。5—11月大風(fēng)日數(shù)均大于平均值, 8月是大風(fēng)出現(xiàn)最多的月份, 大風(fēng)日數(shù)達(dá)15.4 d, 平均每兩天出現(xiàn)一次大風(fēng)天氣。1月出現(xiàn)大風(fēng)天氣的次數(shù)最少, 大風(fēng)日數(shù)僅6.6 d。月大風(fēng)日數(shù)與月平均風(fēng)速呈正相關(guān),r為0.83, 并通過0.001顯著性水平檢驗(yàn)。中山站春、夏、秋、冬季大風(fēng)日數(shù)分別為36.8、25.8、34.4、43.6 d。大風(fēng)日數(shù)也具有明顯的季節(jié)變化特征, 冬季大風(fēng)日數(shù)最多, 春季次之, 秋季較少, 夏季最少。春季和冬季大風(fēng)日數(shù)最多分別為62 d和72 d, 均出現(xiàn)在1999年。夏季大風(fēng)日數(shù)最多為39 d, 出現(xiàn)在1996年。秋季大風(fēng)日數(shù)最多為48 d, 出現(xiàn)在1990年、1991和2000年。春、夏、秋、冬季大風(fēng)日數(shù)和風(fēng)速呈正相關(guān),r分別為0.65、0.75、0.81和0.85,均通過0.001顯著性水平檢驗(yàn)。

    1990—2020年中山站四季大風(fēng)日數(shù)與年大風(fēng)日數(shù)變化一致均呈減少趨勢, 其中夏季大風(fēng)日數(shù)減少趨勢最為顯著, 氣候傾向率和趨勢系數(shù)分別為-8.2 d·(10 a)-1、-0.09, 夏季是年平均風(fēng)速減小最顯著的季節(jié)。秋季和春季大風(fēng)日數(shù)氣候傾向率和趨勢系數(shù)分別為-6.3 d·(10 a)-1、-0.07和-5.0 d·(10 a)-1、-0.06。冬季大風(fēng)日數(shù)減小趨勢最弱, 氣候傾向率和趨勢系數(shù)為-6.1 d·(10 a)-1、-0.05。顯著性檢驗(yàn)表明, 四季大風(fēng)日數(shù)減小趨勢均通過0.01顯著性水平檢驗(yàn), 其中夏季和秋季變化趨勢通過0.001顯著性水平檢驗(yàn)。

    5 風(fēng)速的變化趨勢

    表3為中山站日平均風(fēng)速、日最大風(fēng)速在各風(fēng)速范圍內(nèi)的平均日數(shù)、氣候傾向率和趨勢系數(shù)。中山站日平均風(fēng)速主要出現(xiàn)在3.4~10.7 m·s-1之間, 即以3~5級(jí)風(fēng)為主, 年平均日數(shù)為245 d, 占全年總?cè)諗?shù)的67%, 其中3級(jí)、4級(jí)和5級(jí)風(fēng)分別占年總?cè)諗?shù)的21.3%、25.6%和20.1%。1990—2020年中山站日平均風(fēng)速≤7.9 m·s-1(≤4級(jí))日數(shù)呈增加趨勢, 其中日平均風(fēng)速1.6~3.3 m·s-1(2級(jí))日數(shù)增加趨勢最顯著, 氣候傾向率為6.6 d·(10 a)-1,通過0.01顯著性水平檢驗(yàn)。日平均風(fēng)速≥8.0 m·s-1(≥5級(jí))日數(shù)呈減小趨勢, 其中 5級(jí)風(fēng)(8.0~10.7 m·s-1)日數(shù)減少趨勢最為顯著, 氣候傾向率為-7.6 d·(10 a)-1, 通過0.01顯著性水平檢驗(yàn)。四季日平均風(fēng)速不同風(fēng)力等級(jí)日數(shù)的分布特征與年的分布特征基本一致, 但變化趨勢沒有年日平均風(fēng)速變化趨勢顯著。

    表3 中山站日平均風(fēng)速和日最大風(fēng)速在各風(fēng)速范圍內(nèi)的平均日數(shù)(d)及氣候傾向率[d·(10 a)-1]/趨勢系數(shù)Table 3.Average days (d) of daily wind speed and maximum wind speed in different speed range and climatic tendency rate[d·(10 a)-1]/trend coefficient of Zhongshan Station

    對(duì)于日最大風(fēng)速而言, 風(fēng)速≤5級(jí)日數(shù)呈增加趨勢, 其中日最大風(fēng)速4級(jí)日數(shù)增加趨勢最顯著,氣候傾向率為5.6 d·(10 a)-1。日最大風(fēng)速≥6級(jí)日數(shù)呈減小趨勢,6級(jí)風(fēng)是日最大風(fēng)速一年中出現(xiàn)日數(shù)最多的風(fēng)速段(年平均日數(shù)為97.4 d), 但變化趨勢未能通過0.05顯著性水平檢驗(yàn)。日最大風(fēng)速7級(jí)日數(shù)減小趨勢最顯著, 氣候傾向率為-6.7 d·(10 a)-1,通過0.01顯著性水平檢驗(yàn)。四季日最大風(fēng)速不同風(fēng)力等級(jí)日數(shù)的分布特征與年分布特征基本一致, 但變化趨勢沒有年日最大風(fēng)速變化趨勢顯著。

    通過對(duì)中山站日平均風(fēng)速和日最大風(fēng)速風(fēng)力等級(jí)日數(shù)變化趨勢分析, 初步認(rèn)為, 中山站中小風(fēng)日數(shù)(日平均風(fēng)速<8.0 m·s-1)增加, 中大風(fēng)或大風(fēng)日數(shù)(日平均風(fēng)速≥8.0 m·s-1)減少。年平均風(fēng)速減小主要是由中大風(fēng)和大風(fēng)日數(shù)減少造成的。

    6 風(fēng)速變化與大氣環(huán)流異常的關(guān)系

    圖1表明中山站年平均風(fēng)速在2002年發(fā)生突變性減小, 對(duì)影響風(fēng)速變化的大氣環(huán)流背景場進(jìn)行分析。選取海平面氣壓場、表面氣溫場以及高空高度場做相關(guān)分析(圖9), 并對(duì)2002年前后中山站附近的大氣環(huán)流場做合成分析(圖10)。圖9a為中山站年平均風(fēng)速與地面海平面氣壓的相關(guān)分布, 由圖可見, 在南極內(nèi)陸和普里茲灣海域北部的南印度洋海平面氣壓場存在顯著的相關(guān)區(qū), 中山站風(fēng)速與極地高壓呈顯著的正相關(guān),r超過0.5;與南印度洋海平面氣壓呈負(fù)相關(guān),r低于-0.35, 均通過0.05水平顯著性檢驗(yàn)。圖9b可以看出, 中山站年平均風(fēng)速與南極內(nèi)陸和普里茲灣海域北部的南印度洋表面氣溫存在顯著的負(fù)相關(guān), 與南極戴維斯海附近表面氣溫呈顯著的正相關(guān)。從海平面氣壓和風(fēng)差值場(圖10a)來看, 中山站南部極地高壓負(fù)距平, 北部南印度洋低壓為正距平, 極地冷高壓和繞極低壓帶強(qiáng)度減弱, 高、低壓之間氣壓梯度減小, 空氣流動(dòng)性減弱, 造成中山站平均風(fēng)速減小。影響中山站的氣旋活動(dòng)主要為西北路徑,中山站北部南印度洋的海平面氣壓正異常,使中山站不易受氣旋影響[20]。彼得森淺灘附近海平面氣壓呈負(fù)距平, 形成順時(shí)針環(huán)流風(fēng)場, 南極瑪麗皇后地附近偏南風(fēng)增強(qiáng),在戴維斯海附近形成負(fù)變溫中心。表面氣溫差值場上南極內(nèi)陸和普里茲灣海域北部的南印度洋表面氣溫為正變溫, 而且南極內(nèi)陸表面增溫幅度大于南印度洋表面。南半球環(huán)狀?;蚰蠘O濤動(dòng)(Antarctic Oscillation, AAO)與中山站風(fēng)速呈顯著負(fù)相關(guān), 南極濤動(dòng)的正位相會(huì)引起中山站西風(fēng)和北風(fēng)異常[21]。2002—2020年AAO指數(shù)正位相較1990—2001年數(shù)值增大, 中山站附近北風(fēng)增強(qiáng), 與觀測數(shù)據(jù)一致。

    圖9 中山站年平均風(fēng)速與同期海平面氣壓(a)、表面氣溫場(b)以及700 hPa高度場(c)的相關(guān)系數(shù)。圓點(diǎn)表示通過了0.05水平顯著性檢驗(yàn)Fig.9.Correlation coefficient of the annual mean wind speed at Zhongshan Station with sea level pressure (a), surface termperature filed (b) and 700 hPa height filed of the same period (c).Dot means that it passed the 0.05 reliability significance test

    圖9c為年平均風(fēng)速與700 hPa高度場的相關(guān)分布, 可以看出, 中山站風(fēng)速與普里茲灣海域北部的南印度洋700 hPa高度場呈顯著負(fù)相關(guān),r超過0.5, 與極地內(nèi)陸高度場雖呈正相關(guān), 但相關(guān)性不顯著。年平均風(fēng)速與500 hPa高度場的相關(guān)分布與700 hPa分布基本一致。從高空位勢高度和風(fēng)差值場來看, 700 hPa(圖10b)普里茲灣海域北部的南印度洋附近為正變高, 南極內(nèi)陸為負(fù)變高, 這樣的變高場配置, 使南印度洋附近低壓和極地高壓均減弱。中山站東部的瑪麗皇后地高度場為正變高, 形成弱高壓中心, 有利于中山站偏北氣流的加強(qiáng)和維持。500 hPa(圖10c)南極內(nèi)陸高原和南印度洋高度場均呈正變高, 在兩高之間, 形成順時(shí)針的風(fēng)場, 中山站附近經(jīng)向環(huán)流加強(qiáng)。

    圖10 2002—2020年與1990—2001年合成場之差。a)海平面氣壓(hPa)和風(fēng)場(m·s-1); b)700 hPa高度場(gpm)和風(fēng)場(m·s-1); c) 500 hPa高度場(gpm)和風(fēng)場(m·s-1)Fig.10.Difference of composite fields between 2002—2020 and 1990—2001.a) sea level pressure (hPa) and wind field(m·s-1); b) 700 hPa height field (gpm) and wind field (m·s-1); c) 500 hPa height field (gpm) and wind field (m·s-1)

    7 結(jié)論

    1.中山站月平均風(fēng)速隨時(shí)間呈倒“U”字型,1月平均風(fēng)速最小, 6月和8月風(fēng)速達(dá)到全年最大。風(fēng)速存在明顯的季節(jié)變化, 冬季風(fēng)速最大,秋季次之, 夏季最小。中山站風(fēng)速日變化呈單峰單谷型, 峰值出現(xiàn)在凌晨5時(shí), 谷值出現(xiàn)在15時(shí)。從季節(jié)變化來看, 春季和夏季風(fēng)速具有明顯日變化特征, 秋季和冬季風(fēng)速日變化不明顯。中山站全年風(fēng)向變化較小, 2—11月東風(fēng)出現(xiàn)頻率最高,1月和12月東東北風(fēng)出現(xiàn)頻率最大。

    2.中山站年最多風(fēng)向以東北至東東南為主,東風(fēng)出現(xiàn)頻率最高。近31年風(fēng)速呈減小趨勢, 四個(gè)季節(jié)中夏季風(fēng)速減小趨勢最大。中山站年下降風(fēng)頻數(shù)也呈減少趨勢, 其中夏季和冬季下降風(fēng)頻數(shù)減少趨勢最為顯著。中山站緯向風(fēng)速明顯大于經(jīng)向風(fēng)速, 緯向風(fēng)速呈減小趨勢, 而經(jīng)向風(fēng)速變化趨勢不明顯。南風(fēng)分量年平均風(fēng)速減小速率都顯著大于其他風(fēng)分量, 西風(fēng)和北風(fēng)變化趨勢不顯著, 東風(fēng)分量在夏季減小趨勢最顯著。

    3.近31年中山站大風(fēng)日數(shù)和各季節(jié)總大風(fēng)日數(shù)均呈明顯的減少趨勢。夏季大風(fēng)日數(shù)減少趨勢最大, 冬季是大風(fēng)日數(shù)最多的季節(jié), 但大風(fēng)日數(shù)減少趨勢最弱。中山站中小風(fēng)日數(shù)增加了, 中大風(fēng)或大風(fēng)日數(shù)減少了, 年平均風(fēng)速減小主要由中大風(fēng)和大風(fēng)日數(shù)減少造成的。

    4.中山站年平均風(fēng)速在2002年發(fā)生突變,突變后風(fēng)速減小。從大氣環(huán)流的變化來看, 2002年以后中山站南部極地冷高壓和北部南印度洋繞極低壓帶強(qiáng)度減弱, 導(dǎo)致高、低壓之間氣壓梯度減小, 造成中山站平均風(fēng)速減小。高空經(jīng)向環(huán)流加強(qiáng), 造成中山站附近緯向風(fēng)減弱, 經(jīng)向風(fēng)增強(qiáng)。

    猜你喜歡
    中山站緯向日數(shù)
    南極中山站和長城站日照時(shí)數(shù)的氣候變化特征及其影響因素
    漢江上游漢中區(qū)域不同等級(jí)降水日數(shù)的氣候變化特征分析
    綠色科技(2022年16期)2022-09-15 03:04:46
    紗線強(qiáng)力對(duì)純棉平紋面料強(qiáng)力的影響
    利用掩星溫度數(shù)據(jù)推算大氣月平均緯向風(fēng)場
    天津市濱海新區(qū)塘沽地域雷暴日數(shù)變化規(guī)律及特征分析
    天津科技(2020年2期)2020-03-03 05:09:48
    溫度對(duì)絲綢面料粘襯熱縮率的影響
    絲綢(2018年10期)2018-10-15 09:54:16
    南極洲——中山站
    “雪龍”號(hào)搭載南極科考隊(duì)抵達(dá)中山站
    今日華人(2018年2期)2018-01-16 18:23:33
    柞蠶絲面料在粘襯過程中的熱縮率變化分析
    海南省雷暴日數(shù)年代際變化特征
    18禁在线播放成人免费| 乱人视频在线观看| 在现免费观看毛片| 精品熟女少妇av免费看| aaaaa片日本免费| 国产精品,欧美在线| 婷婷精品国产亚洲av| 欧美成人免费av一区二区三区| 日韩av在线大香蕉| 久久欧美精品欧美久久欧美| 中文字幕av在线有码专区| 久久99热6这里只有精品| 亚洲综合色惰| 午夜免费男女啪啪视频观看 | av专区在线播放| 亚洲婷婷狠狠爱综合网| 免费看a级黄色片| 精品久久久噜噜| 亚洲欧美清纯卡通| 91av网一区二区| 一级毛片aaaaaa免费看小| 69av精品久久久久久| 久久久久免费精品人妻一区二区| 亚洲精品国产av成人精品 | 一级毛片aaaaaa免费看小| 欧美成人a在线观看| 日韩中字成人| 少妇的逼水好多| 亚洲av一区综合| 男女啪啪激烈高潮av片| 国产成人freesex在线 | 六月丁香七月| 寂寞人妻少妇视频99o| 精品久久久久久久久久久久久| 偷拍熟女少妇极品色| 国产成人精品久久久久久| 国产又黄又爽又无遮挡在线| 欧美成人一区二区免费高清观看| 亚洲精品在线观看二区| 真实男女啪啪啪动态图| 久久韩国三级中文字幕| 久久国产乱子免费精品| 欧美日韩国产亚洲二区| 日本a在线网址| 亚洲av免费在线观看| 国产黄片美女视频| 欧美极品一区二区三区四区| 51国产日韩欧美| 国国产精品蜜臀av免费| 国产精品久久久久久久电影| 特级一级黄色大片| 亚洲人成网站高清观看| 给我免费播放毛片高清在线观看| 黑人高潮一二区| 不卡视频在线观看欧美| 免费人成视频x8x8入口观看| 日韩中字成人| 欧美不卡视频在线免费观看| 尤物成人国产欧美一区二区三区| 俺也久久电影网| 免费高清视频大片| 老熟妇乱子伦视频在线观看| 夜夜夜夜夜久久久久| 久久久久久国产a免费观看| 不卡一级毛片| 亚洲精华国产精华液的使用体验 | 欧美潮喷喷水| 国产精品美女特级片免费视频播放器| 九九热线精品视视频播放| 国产精品综合久久久久久久免费| 久久久久久久亚洲中文字幕| 国产人妻一区二区三区在| 99热只有精品国产| 亚洲第一区二区三区不卡| 精品人妻一区二区三区麻豆 | 99热这里只有是精品50| 18禁在线无遮挡免费观看视频 | 国产成人一区二区在线| 国产精品久久久久久久久免| 九九热线精品视视频播放| 91午夜精品亚洲一区二区三区| 麻豆精品久久久久久蜜桃| 亚洲性夜色夜夜综合| 国产av在哪里看| 尾随美女入室| 精品日产1卡2卡| 欧美一区二区亚洲| 成人三级黄色视频| 日韩国内少妇激情av| 久久久a久久爽久久v久久| 18禁裸乳无遮挡免费网站照片| 日日撸夜夜添| 成年版毛片免费区| 一卡2卡三卡四卡精品乱码亚洲| 国产欧美日韩精品亚洲av| 人人妻人人澡欧美一区二区| 亚洲欧美成人精品一区二区| 国产淫片久久久久久久久| 国产一区亚洲一区在线观看| 亚洲成av人片在线播放无| 国产精品久久久久久久久免| 欧美在线一区亚洲| 国产精品一区二区三区四区久久| 欧美一区二区精品小视频在线| 欧美人与善性xxx| 亚洲人成网站在线观看播放| 亚洲不卡免费看| 中出人妻视频一区二区| 99在线视频只有这里精品首页| 免费看美女性在线毛片视频| 黄色一级大片看看| 精品少妇黑人巨大在线播放 | 我要看日韩黄色一级片| 亚洲av.av天堂| 久久久精品欧美日韩精品| 亚洲av一区综合| 嫩草影视91久久| av在线观看视频网站免费| 国产免费一级a男人的天堂| 国产一区二区三区在线臀色熟女| 亚洲av免费在线观看| 日本精品一区二区三区蜜桃| 久久久久久国产a免费观看| 校园人妻丝袜中文字幕| 一级毛片久久久久久久久女| 国内精品宾馆在线| 少妇熟女欧美另类| 亚洲精品乱码久久久v下载方式| 精品久久久久久成人av| 人妻丰满熟妇av一区二区三区| 淫妇啪啪啪对白视频| 在线播放无遮挡| 亚洲人成网站在线观看播放| 欧美人与善性xxx| 黄色日韩在线| 国产午夜福利久久久久久| 精品人妻视频免费看| 亚洲国产欧洲综合997久久,| 久久久久久久午夜电影| 男女视频在线观看网站免费| 久久精品国产亚洲网站| 高清毛片免费看| 精品国内亚洲2022精品成人| 尤物成人国产欧美一区二区三区| 国产在线男女| 在线观看66精品国产| 国产伦精品一区二区三区四那| 免费人成在线观看视频色| 99久久精品一区二区三区| 欧美日韩国产亚洲二区| 夜夜看夜夜爽夜夜摸| 精品人妻视频免费看| 亚洲自偷自拍三级| 国产 一区精品| 99久久久亚洲精品蜜臀av| 国产乱人视频| 99热网站在线观看| 男女边吃奶边做爰视频| 国产欧美日韩精品亚洲av| 国产男人的电影天堂91| 日韩精品有码人妻一区| 亚洲自偷自拍三级| 听说在线观看完整版免费高清| 自拍偷自拍亚洲精品老妇| 97人妻精品一区二区三区麻豆| 婷婷色综合大香蕉| 一进一出好大好爽视频| 国产高清视频在线观看网站| 国产亚洲91精品色在线| 国产成人91sexporn| 免费不卡的大黄色大毛片视频在线观看 | 免费av毛片视频| 免费观看在线日韩| 久久久a久久爽久久v久久| 成年av动漫网址| 欧美成人免费av一区二区三区| 国产黄色视频一区二区在线观看 | 一个人看的www免费观看视频| 成人鲁丝片一二三区免费| 乱系列少妇在线播放| 欧美成人a在线观看| 中国国产av一级| 午夜a级毛片| 色在线成人网| 日本黄大片高清| 乱码一卡2卡4卡精品| 亚洲av电影不卡..在线观看| 99久久精品国产国产毛片| 欧美成人一区二区免费高清观看| 免费大片18禁| 97在线视频观看| 亚洲乱码一区二区免费版| 久久热精品热| 在线观看av片永久免费下载| 一卡2卡三卡四卡精品乱码亚洲| 久久天躁狠狠躁夜夜2o2o| 中文字幕熟女人妻在线| 国产午夜精品久久久久久一区二区三区 | 永久网站在线| 免费一级毛片在线播放高清视频| 国产日本99.免费观看| 一级毛片我不卡| 日日撸夜夜添| 搡老岳熟女国产| 嫩草影视91久久| 日韩,欧美,国产一区二区三区 | 看十八女毛片水多多多| 大又大粗又爽又黄少妇毛片口| 成年免费大片在线观看| 精品久久国产蜜桃| 我的女老师完整版在线观看| 2021天堂中文幕一二区在线观| 老师上课跳d突然被开到最大视频| 亚洲av成人精品一区久久| 中文资源天堂在线| 亚洲第一区二区三区不卡| 一区二区三区高清视频在线| www日本黄色视频网| 神马国产精品三级电影在线观看| 精品熟女少妇av免费看| 人人妻人人澡人人爽人人夜夜 | 国产黄a三级三级三级人| 久久热精品热| 色哟哟哟哟哟哟| 日本成人三级电影网站| 国产黄片美女视频| h日本视频在线播放| 国产精品乱码一区二三区的特点| 午夜福利高清视频| 老熟妇仑乱视频hdxx| av在线蜜桃| 日韩三级伦理在线观看| 一本久久中文字幕| 秋霞在线观看毛片| 亚洲国产精品久久男人天堂| 97超碰精品成人国产| 熟妇人妻久久中文字幕3abv| 亚洲中文字幕一区二区三区有码在线看| 白带黄色成豆腐渣| 成人av一区二区三区在线看| 国产在线男女| 赤兔流量卡办理| 尾随美女入室| 人妻少妇偷人精品九色| 99久国产av精品国产电影| 亚洲欧美日韩卡通动漫| 麻豆国产97在线/欧美| 欧美+亚洲+日韩+国产| 99视频精品全部免费 在线| 在线观看免费视频日本深夜| 你懂的网址亚洲精品在线观看 | 欧美日韩一区二区视频在线观看视频在线 | 变态另类成人亚洲欧美熟女| 午夜福利高清视频| 久久亚洲精品不卡| 99在线人妻在线中文字幕| 午夜亚洲福利在线播放| 亚洲av免费高清在线观看| 久久99热6这里只有精品| av在线老鸭窝| 欧美激情国产日韩精品一区| 在线观看免费视频日本深夜| av中文乱码字幕在线| 国产精品国产三级国产av玫瑰| 成人漫画全彩无遮挡| 99久国产av精品国产电影| 国产v大片淫在线免费观看| 免费观看精品视频网站| 国产黄片美女视频| 婷婷精品国产亚洲av在线| 国产一区二区亚洲精品在线观看| 一个人观看的视频www高清免费观看| 午夜免费激情av| 如何舔出高潮| 小蜜桃在线观看免费完整版高清| 黄色日韩在线| 国产男靠女视频免费网站| 97超视频在线观看视频| 超碰av人人做人人爽久久| 最近手机中文字幕大全| 久久99热6这里只有精品| 国产色爽女视频免费观看| 亚洲,欧美,日韩| 22中文网久久字幕| 丰满的人妻完整版| 国产一区二区亚洲精品在线观看| 精品久久久久久久末码| 男人狂女人下面高潮的视频| 国产精品野战在线观看| 三级男女做爰猛烈吃奶摸视频| 99国产极品粉嫩在线观看| 国产欧美日韩精品亚洲av| 亚洲欧美日韩高清专用| 成人精品一区二区免费| 亚洲乱码一区二区免费版| 99国产精品一区二区蜜桃av| 中文字幕精品亚洲无线码一区| 色尼玛亚洲综合影院| 国产精品亚洲美女久久久| 国产乱人视频| 在线免费观看的www视频| 国产精品美女特级片免费视频播放器| 国产高清三级在线| 久久鲁丝午夜福利片| 久久久久久九九精品二区国产| 久久精品人妻少妇| 黑人高潮一二区| 性插视频无遮挡在线免费观看| 中文字幕精品亚洲无线码一区| 亚洲在线自拍视频| 国产免费男女视频| 久久久午夜欧美精品| 天堂√8在线中文| 午夜免费激情av| 精品99又大又爽又粗少妇毛片| 特级一级黄色大片| 大型黄色视频在线免费观看| 日本免费a在线| 狂野欧美白嫩少妇大欣赏| 真实男女啪啪啪动态图| 成年版毛片免费区| 日本与韩国留学比较| 免费高清视频大片| 亚洲精品久久国产高清桃花| 91精品国产九色| 亚洲性久久影院| 国产精品久久久久久亚洲av鲁大| 日日干狠狠操夜夜爽| 成人一区二区视频在线观看| 成人特级黄色片久久久久久久| 国产一区二区在线观看日韩| 18禁在线播放成人免费| 精品无人区乱码1区二区| 99精品在免费线老司机午夜| 国产亚洲精品久久久久久毛片| 亚洲精品成人久久久久久| 亚洲av一区综合| 美女 人体艺术 gogo| 美女被艹到高潮喷水动态| av.在线天堂| 国产精品av视频在线免费观看| 我要搜黄色片| 国产精品一二三区在线看| 久久久a久久爽久久v久久| 一个人看视频在线观看www免费| 99国产极品粉嫩在线观看| 欧美一区二区国产精品久久精品| 在现免费观看毛片| 成年av动漫网址| 日日摸夜夜添夜夜添av毛片| 精品一区二区三区人妻视频| 亚洲18禁久久av| 男人和女人高潮做爰伦理| 男女那种视频在线观看| 国产精品1区2区在线观看.| 久久久久久伊人网av| 国模一区二区三区四区视频| 亚洲精品乱码久久久v下载方式| 欧美+日韩+精品| 亚洲自拍偷在线| 在线观看美女被高潮喷水网站| 性欧美人与动物交配| 欧美绝顶高潮抽搐喷水| eeuss影院久久| 一区二区三区四区激情视频 | 午夜激情福利司机影院| 九九在线视频观看精品| 长腿黑丝高跟| 97超级碰碰碰精品色视频在线观看| 国产精品人妻久久久久久| 国产伦精品一区二区三区视频9| 六月丁香七月| 国产探花在线观看一区二区| 一区二区三区高清视频在线| 一级毛片久久久久久久久女| 亚洲第一电影网av| 国产精品亚洲一级av第二区| 色吧在线观看| 亚洲国产精品成人久久小说 | 淫妇啪啪啪对白视频| 免费在线观看影片大全网站| 成人漫画全彩无遮挡| 成熟少妇高潮喷水视频| 欧美区成人在线视频| 国产一区二区在线av高清观看| 少妇熟女欧美另类| 婷婷精品国产亚洲av| 国产国拍精品亚洲av在线观看| 一级av片app| av在线蜜桃| 欧美性感艳星| 99热这里只有精品一区| 熟女电影av网| 激情 狠狠 欧美| 色av中文字幕| 欧美又色又爽又黄视频| 免费观看的影片在线观看| 久久久欧美国产精品| 99热网站在线观看| 成年免费大片在线观看| 人妻久久中文字幕网| 精品一区二区三区视频在线观看免费| 一a级毛片在线观看| 99久国产av精品国产电影| 成人美女网站在线观看视频| 日产精品乱码卡一卡2卡三| 毛片一级片免费看久久久久| 黄色欧美视频在线观看| 人人妻人人澡欧美一区二区| 久久6这里有精品| av在线老鸭窝| 亚洲精品在线观看二区| 性插视频无遮挡在线免费观看| 男女那种视频在线观看| 蜜臀久久99精品久久宅男| 天天躁日日操中文字幕| 久久人妻av系列| 亚洲av中文av极速乱| 日本a在线网址| 十八禁国产超污无遮挡网站| 精品人妻视频免费看| 久久久久国内视频| 天堂网av新在线| 午夜福利18| 少妇猛男粗大的猛烈进出视频 | 长腿黑丝高跟| 亚洲一区高清亚洲精品| 一级黄色大片毛片| 噜噜噜噜噜久久久久久91| 精品99又大又爽又粗少妇毛片| 国产一级毛片七仙女欲春2| 国产精品精品国产色婷婷| 久久精品国产鲁丝片午夜精品| 国产精品综合久久久久久久免费| 国产成人影院久久av| 99久久精品一区二区三区| avwww免费| 少妇熟女欧美另类| 一级黄色大片毛片| АⅤ资源中文在线天堂| 亚洲av不卡在线观看| 我的老师免费观看完整版| 小蜜桃在线观看免费完整版高清| 在线天堂最新版资源| 激情 狠狠 欧美| 日韩欧美国产在线观看| 免费看a级黄色片| 六月丁香七月| 国产成人freesex在线 | 久久精品国产99精品国产亚洲性色| 国产精品亚洲一级av第二区| 亚洲国产欧美人成| 色av中文字幕| 国产爱豆传媒在线观看| 免费搜索国产男女视频| 中文字幕人妻熟人妻熟丝袜美| av天堂在线播放| 一本一本综合久久| 大香蕉久久网| 国产精品日韩av在线免费观看| 女的被弄到高潮叫床怎么办| 国产乱人视频| 亚洲va在线va天堂va国产| 色综合色国产| 亚洲av第一区精品v没综合| 国产成人一区二区在线| 久久精品综合一区二区三区| 亚洲第一电影网av| 一a级毛片在线观看| 非洲黑人性xxxx精品又粗又长| 变态另类成人亚洲欧美熟女| 午夜亚洲福利在线播放| 91在线观看av| 亚洲精品色激情综合| 欧美高清性xxxxhd video| 亚洲乱码一区二区免费版| 熟女人妻精品中文字幕| 午夜福利高清视频| 国内精品宾馆在线| 国产老妇女一区| 欧美在线一区亚洲| 久久亚洲国产成人精品v| 亚洲自偷自拍三级| 直男gayav资源| 欧美日本亚洲视频在线播放| 性插视频无遮挡在线免费观看| 在线播放无遮挡| 床上黄色一级片| 日本一本二区三区精品| 欧美日韩精品成人综合77777| 精品欧美国产一区二区三| 亚洲在线自拍视频| 亚洲精品成人久久久久久| 欧美xxxx性猛交bbbb| 国产伦一二天堂av在线观看| 麻豆国产97在线/欧美| 亚洲久久久久久中文字幕| 国产老妇女一区| 两个人视频免费观看高清| 成人高潮视频无遮挡免费网站| 亚洲av熟女| 国产精品不卡视频一区二区| 色哟哟·www| 国内久久婷婷六月综合欲色啪| 深夜a级毛片| a级毛色黄片| 中文资源天堂在线| 男女边吃奶边做爰视频| 看黄色毛片网站| 国产av不卡久久| 级片在线观看| 国产真实伦视频高清在线观看| 精品欧美国产一区二区三| 在线天堂最新版资源| av视频在线观看入口| 一级毛片aaaaaa免费看小| 亚洲av免费在线观看| 蜜桃久久精品国产亚洲av| av中文乱码字幕在线| 国产精品美女特级片免费视频播放器| 亚洲中文字幕日韩| 18禁在线无遮挡免费观看视频 | 十八禁网站免费在线| 国产大屁股一区二区在线视频| 亚洲成人中文字幕在线播放| 亚洲成人久久爱视频| 老熟妇仑乱视频hdxx| 最近在线观看免费完整版| 你懂的网址亚洲精品在线观看 | 在线天堂最新版资源| 99国产精品一区二区蜜桃av| 亚洲性久久影院| 一级毛片电影观看 | 国产极品精品免费视频能看的| 69av精品久久久久久| 欧美日本视频| 女生性感内裤真人,穿戴方法视频| 亚洲在线自拍视频| 免费电影在线观看免费观看| 观看免费一级毛片| 久久这里只有精品中国| 蜜桃亚洲精品一区二区三区| 国产精品福利在线免费观看| 99热6这里只有精品| av在线观看视频网站免费| 99在线视频只有这里精品首页| 寂寞人妻少妇视频99o| 最近的中文字幕免费完整| 久久久午夜欧美精品| 免费看av在线观看网站| 五月伊人婷婷丁香| 国产大屁股一区二区在线视频| 免费不卡的大黄色大毛片视频在线观看 | 最近在线观看免费完整版| 亚洲精品亚洲一区二区| 午夜激情欧美在线| 日本与韩国留学比较| 又爽又黄无遮挡网站| av在线观看视频网站免费| 干丝袜人妻中文字幕| 国产午夜精品久久久久久一区二区三区 | 欧美zozozo另类| 国产精品99久久久久久久久| 免费看美女性在线毛片视频| 少妇的逼水好多| 国产黄a三级三级三级人| 国产成人a区在线观看| 国产v大片淫在线免费观看| 日韩精品有码人妻一区| 国产综合懂色| 我的女老师完整版在线观看| 成人美女网站在线观看视频| 亚洲国产欧美人成| 91狼人影院| 久久6这里有精品| 国产精品亚洲一级av第二区| 国产精品嫩草影院av在线观看| 老熟妇仑乱视频hdxx| 最新在线观看一区二区三区| 国内揄拍国产精品人妻在线| 一本精品99久久精品77| 精品久久久久久久久av| 久久久欧美国产精品| 丝袜美腿在线中文| 永久网站在线| 热99re8久久精品国产| 色吧在线观看| 天天一区二区日本电影三级| 麻豆成人午夜福利视频| 久久精品夜夜夜夜夜久久蜜豆| 国产亚洲精品综合一区在线观看| 国产成人影院久久av| 亚洲成人中文字幕在线播放| 精品久久久久久久人妻蜜臀av| 亚洲va在线va天堂va国产| 18禁黄网站禁片免费观看直播| 国产伦精品一区二区三区四那| 最后的刺客免费高清国语| 国产 一区精品| 99热精品在线国产| 97热精品久久久久久| 成人三级黄色视频| 真人做人爱边吃奶动态| 国产精华一区二区三区| 久久精品久久久久久噜噜老黄 | 99热这里只有是精品50| 免费搜索国产男女视频| 日本与韩国留学比较| 欧美潮喷喷水| 高清毛片免费观看视频网站| 亚洲自拍偷在线| 乱码一卡2卡4卡精品| 日韩强制内射视频| 国产视频一区二区在线看| 女的被弄到高潮叫床怎么办|