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

    1970—2020 年昌吉地區(qū)日照時數(shù)時空變化特征分析

    2021-11-08 01:46:20葉爾克江霍依哈孜阿帕爾肉孜
    沙漠與綠洲氣象 2021年5期
    關(guān)鍵詞:昌吉日照時數(shù)日數(shù)

    葉爾克江·霍依哈孜,黃 健,王 森,阿帕爾·肉孜*

    (1.昌吉州氣象局,新疆 昌吉 831100;2.中國氣象局烏魯木齊沙漠氣象研究所,新疆 烏魯木齊 830002;3.新疆農(nóng)業(yè)氣象臺,新疆 烏魯木齊 830002)

    近些年來,氣候變化越來越受人類的重視,同時區(qū)域性氣候已引起人們普遍關(guān)注[1]。 據(jù)科學(xué)研究表明,平均地表溫度一直在增高,近50 a 來平均地表氣溫增幅約為1.1 ℃,每10 a 增溫速率約為0.22 ℃,21 世紀(jì)以來增溫趨勢更加明顯[2]。氣候變暖背景下,風(fēng)速、氣溫、降水、光照等氣象要素都有不同程度的變化[3]。 日照時數(shù)是太陽直射光線照射地面的時間,是太陽輻射變化的直接反映, 日照時數(shù)的長短與大氣熱力狀況、 動植物生長發(fā)育、 農(nóng)牧業(yè)生產(chǎn)密切相關(guān),是植物光合作用必不可少的條件,是重要的農(nóng)業(yè)氣候資源,也是影響氣候變化的主要要素之一[4-5]。隨著氣候的變暖、 城市的發(fā)展和人為污染物排放的加劇,各地日照時間的長短發(fā)生了巨大的變化。日照時數(shù)時空變化研究對氣候資源開發(fā)利用、 農(nóng)業(yè)結(jié)構(gòu)調(diào)整、農(nóng)業(yè)生產(chǎn)合理布局具有重要意義[6-7]。 眾多學(xué)者對各地區(qū)日照時數(shù)的變化進(jìn)行了研究。 李慧群等[8]指出中國絕大部分地區(qū)近50 a 來日照時數(shù)總量呈減少的趨勢,黃小燕等[9]研究表明中國西北地區(qū)大部分地區(qū)日照時數(shù)呈明顯的減少趨勢,苗運玲[10]、楊艷玲等[11]研究了新疆哈密地區(qū),劉強(qiáng)吉等[12]研究了新疆巴州地區(qū),禹樸家等[13]研究新疆天山北坡,鄭玉萍等[14]研究烏魯木齊市北部農(nóng)區(qū),趙勇等[15]研究新疆天山地區(qū)日照時數(shù)的氣候特征, 結(jié)果顯示日照時數(shù)均呈減少趨勢。 昌吉地區(qū)是新疆主要的農(nóng)牧業(yè)生產(chǎn)地,也是太陽能豐富的地區(qū),研究其日照時數(shù)的變化對指導(dǎo)農(nóng)業(yè)經(jīng)濟(jì)發(fā)展具有重要意義。目前,昌吉地區(qū)日照時數(shù)變化趨勢和規(guī)律及其影響因子的研究較少,因此,本文選取1970—2020 年昌吉地區(qū)10 個氣象站逐月日照時數(shù)資料,對昌吉地區(qū)近51 a 日照時數(shù)時空變化及其影響因素進(jìn)行分析。

    1 研究區(qū)氣候特征及研究方法

    1.1 研究區(qū)域概況

    新疆昌吉回族自治州地處天山北麓、 準(zhǔn)噶爾盆地東南緣,43°20′~45°00′N,85°17′~91°32′E。地勢南高北低,由東南向西北傾斜,南部是天山山地,中部為廣袤的沖積平原,北部為浩瀚的沙漠盆地。昌吉回族自治州轄有5 縣2 市(圖1)。 昌吉位于歐亞大陸腹地、四周環(huán)山的盆地中,降水少,干燥度大,蒸發(fā)強(qiáng)烈,氣溫變化劇烈,為典型的大陸性干旱氣候。 具有冬季寒冷漫長,夏季炎熱,光照充足,熱量豐富,氣溫年較差、日較差大,春、秋季氣溫變化劇烈,春、夏季多風(fēng)的氣候特點。 各地常年平均氣溫在2.3~8.0 ℃,無霜期為145~172 d,≥0 ℃的活動積溫為2 446.3~4 161.8 ℃,≥10 ℃的活動積溫為1 849.6~3 671.8 ℃,年降水量為153.8~343.5 mm,年蒸發(fā)量為1 606.8~2 163.3 mm,年日照時數(shù)為2 701.8~3 075.8 h,年平均大風(fēng)日數(shù)為4~22 d。

    圖1 研究區(qū)域地形

    1.2 研究方法

    1.2.1 數(shù)據(jù)與分析方法

    選用1970—2020 年新疆昌吉回族自治州(下稱昌吉地區(qū))各縣市10 個國家氣象觀測站的日照時數(shù)觀測資料; 同時參照地貌和氣候特點將昌吉地區(qū)劃分為西部、中部、東部和山區(qū)4 個區(qū)域,其中,瑪納斯縣和呼圖壁縣氣象站代表西部,昌吉市、蔡家湖和阜康市氣象站代表中部,吉木薩爾縣、奇臺縣和木壘縣氣象站代表東部,天池和北塔山氣象站代表山區(qū)。季節(jié)的劃分根據(jù)北半球典型的劃分方法,3—5 月為春季,6—8 月為夏季,9—11 月為秋季,12 月—次年2月為冬季。

    本文采用反距離加權(quán)(IDW)插值法、突變檢驗、氣候傾向率和趨勢變化等分析昌吉地區(qū)日照時數(shù)空間分布和趨勢變化特征, 運用Pearson 相關(guān)性分析昌吉地區(qū)10 個氣象站近51 a 總云量、揚沙日數(shù)、水汽壓、 大風(fēng)日數(shù)和霧日數(shù)等5 個氣象因子與日照時數(shù)的相關(guān)性。

    1.2.2 氣候傾向率

    氣候傾向率可以用來顯示氣候要素隨時間的變化趨勢[16],用yi表示樣本量為n 某一氣候變量,用表示xi所對應(yīng)的時間序列,建立yi與xi之間的一元線性回歸方程[17-21]:

    式中,i 為年序列,a 為回歸系數(shù)(線性方程的斜率),也就是氣候要素的線性變化趨勢和速率, 也稱為傾向率。a>0 表示氣候要素隨時間x 的增加氣候變量y呈增加趨勢;若a<0,則說明隨時間x 的增加y 呈下降趨勢。b 為常數(shù),可通過最小二乘法求取[16-21]。本文運用趨勢分析法探究近51 a 昌吉地區(qū)日照趨勢變化,采用一元線性回歸進(jìn)行氣候趨勢傾向估計。

    1.2.3 突變檢驗

    借助Matlab R 2018b 和Origin 2018 等軟件,主要從最常見的Mann-Kendall(簡稱M-K)檢驗[16-22]著手, 再考慮日照時數(shù)演變趨勢和累積距平值的變化情況, 進(jìn)一步檢驗昌吉地區(qū)日照變化趨勢和突變特征。

    1.2.4 反距離權(quán)重插值方法

    運用ArcGIS10.7 中的反距離權(quán)重插值IDW(Inverse Distance Weighted)方法將昌吉地區(qū)日照數(shù)據(jù)進(jìn)行空間插值,得到相應(yīng)的空間分布。 IDW 是一種常用而簡便的空間插值方法, 它以插值點與樣本點間的距離為權(quán)重進(jìn)行加權(quán)平均, 是基于相近相似的原理,每個采樣點都對插值點有一定的影響,即權(quán)重。 權(quán)重隨著采樣點和插值點之前距離的增加而減小,離插值點越近的樣本點賦予的權(quán)重越大,而且當(dāng)采樣點在距離插值點一定范圍以外時, 權(quán)重可以忽略不計[23-24]。 IDW 通過對鄰近區(qū)域的每個采樣點值平均運算獲得內(nèi)插單元值。 IDW 是一個均分過程,這一方法要求離散點均勻分布, 并且密集程度足以滿足在分析中反映局部表面變化。

    設(shè)平面上分布一系列離散點, 任一離散點的值是各采樣點權(quán)重之和, 已知其坐標(biāo)和值為xi,yi,zi,(i=1,2,3,…,n),根據(jù)周圍離散點的值,通過距離權(quán)值求目標(biāo)插值Zp,表示為:

    式中,zi為第i 個采樣點數(shù)值,di為第i 個采樣點到插值點距離,di2=(x-xi)2(y-yi)2,d-μ為距離衰減函數(shù),權(quán)指數(shù)μ 表示距離越大而影響越小。 當(dāng)μ=0 時,距離沒有影響;當(dāng)μ=1 時,距離的影響為線性;當(dāng)μ>1時, 遙遠(yuǎn)位置的影響快速減少。 μ 通常取值為1 或2,但多數(shù)學(xué)者認(rèn)為μ=2 得到的試驗效果更好,本文中μ 值取2。

    1.2.5 Pearson 相關(guān)系數(shù)分析

    采用Pearson 相關(guān)系數(shù)分析法對每一個像元相應(yīng)的年日照時數(shù)與相應(yīng)的各年總云量、揚沙日數(shù)、水汽壓、 大風(fēng)日數(shù)和霧日數(shù)等氣象因子進(jìn)行相關(guān)性分析,以此分析日照時數(shù)與氣候因子之間的響應(yīng)關(guān)系。Pearson 相關(guān)分析是一種準(zhǔn)確度量2 個變量之間的關(guān)系密切程度的統(tǒng)計學(xué)方法, 系數(shù)大小可以反映2個變量相關(guān)關(guān)系的密切程度, 將兩組氣象要素記為(xi,yi)(i=1,2,…,n),則相關(guān)系數(shù)的數(shù)學(xué)表達(dá)式為[25-26]:

    2 結(jié)果分析

    2.1 日照時數(shù)年際和年代變化

    統(tǒng)計昌吉地區(qū)各縣市10 個國家氣象觀測站近51 a 各年累計日照時數(shù), 分析其變化特征和趨勢(圖2a)。昌吉地區(qū)近51 a 日照時數(shù)呈明顯的減少趨勢,氣候傾向率為-62.8 h/10 a,并通過0.001 的顯著性檢驗,51 a 來減少320 h, 并且1987 年以來減少速率有加快趨勢。昌吉地區(qū)近51 a 日照時數(shù)年際間波動較大, 多年平均為2 840.3 h, 最小值出現(xiàn)在2016 年,為2 531.4 h;最大值出現(xiàn)在1974 年,為3 044.5 h; 有24 a 的日照時數(shù)超過多年平均值,有27 a 的日照時數(shù)低于多年平均值;最大與最小年日照時數(shù)之差為513.1 h。 從5 a 滑動平均曲線來看,近51 a 日照時數(shù)逐年下降, 其中1972—1987 年、1993—2002 年和2011—2015 年3 個階段日照時數(shù)呈現(xiàn)較大的減少趨勢, 其中, 第1 階段減少幅度最大, 其次是第3 階段。 從不同年代日照時數(shù)總量來看,20 世紀(jì)70—90 年代的日照時數(shù)偏多,均高于多年平均值;21 世紀(jì)以來,各年代日照時數(shù)明顯偏少,遠(yuǎn)低于多年平均值(表1)。 在M-K 檢驗中(圖2b),日照時數(shù)M-K 統(tǒng)計曲線和年際變化趨勢線基本吻合, 整體呈減少趨勢,UFk和UBk曲線在1987 年出現(xiàn)交點, 是突變開始的年份,UFk曲線在2002 年超出0.05 顯著性水平臨界線, 在2016 年超出0.01 顯著性水平臨界線,說明昌吉地區(qū)日照時數(shù)逐年減少趨勢十分顯著,1987 年開始突變減少并在2002 年開始減少比較明顯, 突變前多年平均日照時數(shù)比突變后多180.41 h。

    圖2 昌吉地區(qū)日照時數(shù)年際、年代變化(a)及M-K 檢驗(b)

    表1 昌吉地區(qū)年和四季平均日照時數(shù)變化 h

    2.2 日照時數(shù)季節(jié)變化

    對昌吉地區(qū)10 個站近51 a 各年各季平均日照時數(shù)的分析表明, 各季日照時數(shù)均呈現(xiàn)不同程度的下降趨勢(圖3)。 日照時數(shù)春季平均為803.0 h,下降速率為-7.2 h/10 a;夏季平均為912.5 h,下降速率為-11.9 h/10 a; 秋季平均為672.8 h, 下降速率為-12.9 h/10 a; 冬季平均為452.4 h, 下降速率最大,為-31.7 h/10 a;其中,春季下降趨勢未通過顯著性檢驗,變化趨勢不顯著,其它季度通過0.001 的顯著性檢驗。 日照時數(shù)夏季最多,其次是春季和秋季,冬季最少。各季度日照時數(shù)均有下降趨勢,冬季日照時數(shù)下降趨勢最快,其次是秋季和夏季,春季最慢,說明冬季日照時數(shù)的逐年減少對全年日照時數(shù)逐年減少作用最大,是年日照時數(shù)減少的主要貢獻(xiàn)者;各季年代際變化與全年年代際變化基本一致(表1)。

    圖3 昌吉地區(qū)春、夏、秋、冬季日照時數(shù)逐年變化

    2.3 日照時數(shù)月變化

    昌吉地區(qū)10 個站近51 a 各月多年平均日照時數(shù)年內(nèi)變化基本一致,均呈現(xiàn)雙峰型(圖4)。最大值出現(xiàn)在5、7 和8 月,分別為312.1、308.8、302.9 h,最小值出現(xiàn)在12 月, 為128.7 h;1—5 月日照時數(shù)持續(xù)上升,5 月達(dá)到最大值,8 月以后不斷下降, 在12月達(dá)到最小值(圖5)。1—12 月各月的日照時數(shù)均呈現(xiàn)不同速率的減少趨勢, 減少速率為-1.37~-0.15 h/10 a(表2),其中,1 月日照時數(shù)減少最為明顯,近51 a共減少了68.5 h,12 月次之。 各月各站最高日照時數(shù)依次為: 呼圖壁>蔡家湖>阜康>瑪納斯>昌吉>北塔山>木壘>奇臺>吉木薩爾>天池。

    圖4 昌吉地區(qū)各站日照時數(shù)逐月變化

    圖5 昌吉地區(qū)日照時數(shù)逐月變化

    表2 昌吉地區(qū)日照時數(shù)各月變化趨勢

    2.4 日照時數(shù)空間變化特征

    昌吉地區(qū)日照時數(shù)年平均值為2 840.3 h,空間分布呈現(xiàn)“東部多、中西部少,山區(qū)多、平原少”的格局(表3)。 木壘縣和北塔山日照時數(shù)在3 000 h以上,呼圖壁和奇臺縣日照時數(shù)在2 900 h 以上,蔡家湖和吉木薩爾縣日照時數(shù)在2 800 h 以上, 瑪納斯、昌吉市和阜康市日照時數(shù)在2 700 h 以上,天池日照時數(shù)最小,為2 603 h(圖6a)。 從空間變化趨勢來看, 阜康市近51 a 日照時數(shù)減少趨勢最為明顯(圖6b),減少率為-113.4 h/10 a;其次是呼圖壁縣和蔡家湖,減少速率分別為-77.6、-75.9 h/10 a;吉木薩爾縣、奇臺縣和北塔山,減少速率分別為-69.9、-71.4、-70.9 h/10 a; 其余縣市減少速率較慢,在-2.54~-59.8 h/10 a。 中西部地區(qū)的日照時數(shù)減少速率明顯大于東部地區(qū)和山區(qū)。

    圖6 昌吉地區(qū)日照時數(shù)空間分布(a)和變化趨勢空間分布(b)

    表3 昌吉地區(qū)不同區(qū)域日照時數(shù)變化

    2.5 日照時數(shù)的主要影響因子分析

    氣象因子是造成日照時數(shù)變化的主要原因,在全球變暖的氣候背景下, 日照時數(shù)的年際變化與氣象因子變化有密切的關(guān)系, 為深入分析影響日照時數(shù)變化的主要氣象因子, 選取與日照時數(shù)變化密切相關(guān)的總云量、揚沙日數(shù)、水汽壓、大風(fēng)日數(shù)和霧日數(shù)等氣象因子,并對選取的氣象因子進(jìn)行趨勢分析,對氣象因子與日照時數(shù)進(jìn)行Pearson 相關(guān)性分析。結(jié)果顯示(表4),各氣象因子變化趨勢有一定差異,總云量中西部呈增加趨勢,東部和山區(qū)呈減少趨勢;揚沙日數(shù)和大風(fēng)日數(shù)各地均呈減少趨勢; 水汽壓和霧日數(shù)各地均呈增加趨勢。昌吉地區(qū)總云量、水汽壓和霧日數(shù)呈增加趨勢,增加趨勢顯著,其中霧日數(shù)增加速率最大,為11.72 d/10 a;揚沙日數(shù)和大風(fēng)日數(shù)呈減少趨勢,減少趨勢顯著。

    表4 昌吉地區(qū)各站總云量、揚沙日數(shù)、水汽壓、大風(fēng)日數(shù)和霧日數(shù)的變化

    各站日照時數(shù)與氣象因子相關(guān)分析結(jié)果表明(表5),昌吉地區(qū)日照時數(shù)與總云量、水汽壓和霧日數(shù)呈顯著負(fù)相關(guān),其中水汽壓和霧日數(shù)通過0.01 的顯著性檢驗,與大風(fēng)日數(shù)和揚沙日數(shù)呈顯著正相關(guān),并通過0.01 的顯著性檢驗。

    表5 昌吉地區(qū)各站日照時數(shù)與有關(guān)因子之間的相關(guān)性

    日照時數(shù)與水汽壓和霧日數(shù)的負(fù)相關(guān)性最強(qiáng),這表明水汽壓逐年增加和霧日數(shù)的增多是日照時數(shù)減少的主要原因。 昌吉地區(qū)各地日照時數(shù)減少趨勢明顯,其中阜康較突出,從瑪納斯到奇臺是水汽壓和霧日數(shù)分布較多區(qū)域,近51 a 昌吉地區(qū)各地水汽壓顯著增加趨勢, 水汽壓的增加表明大氣中的水汽含量增多,水汽對太陽輻射的散射、折射及吸收等直接影響著日照時數(shù), 大氣中水汽含量的高低也直接影響霧日數(shù)的多或少(尤其是冬季),這表明水汽壓和霧日數(shù)顯著增加是影響昌吉地區(qū)日照時數(shù)減少的主要因素。另外日照時數(shù)與大風(fēng)日數(shù)呈顯著正相關(guān),大風(fēng)日數(shù)變少引起云層易集聚、污染物擴(kuò)散慢、不同程度地降低大氣透明度,削弱太陽輻射,進(jìn)而降低日照時數(shù), 因此大風(fēng)日數(shù)也是影響昌吉地區(qū)日照時數(shù)的因素之一。

    日照時數(shù)與水汽壓和霧日數(shù)呈顯著負(fù)相關(guān),并通過0.01 的顯著性檢驗,與大風(fēng)日數(shù)和揚沙日數(shù)呈顯著正相關(guān),并通過0.01 的顯著性檢驗。 水汽壓的增大、霧日數(shù)的增多、大風(fēng)日數(shù)的減少是造成昌吉地區(qū)日照時數(shù)減少的主要原因。

    3 討論

    昌吉地區(qū)近51 a 日照時數(shù)呈減少趨勢,這與李慧群、陳少勇、張山清[27-29]等研究結(jié)論一致。 昌吉地區(qū)51 a 水汽壓的增大、霧日數(shù)的增多、大風(fēng)日數(shù)的減少是造成昌吉地區(qū)日照時數(shù)減少的主要原因,這與肖風(fēng)勁[30]、虞海燕[31]、黃小燕[32]、張立波[33]、武鵬飛等[34]研究日照時數(shù)減少與降水、空氣濕度、風(fēng)、能見度有關(guān)的結(jié)論基本一致。本研究發(fā)現(xiàn)水汽壓的增大、霧日數(shù)增多、大風(fēng)日數(shù)減少對日照時數(shù)減少作用大;冬季日照時數(shù)減少的最快, 且對全年日照時數(shù)減少的影響最大;這很可能是由于當(dāng)?shù)毓I(yè)活動增加,大氣污染排放增加, 環(huán)境中氣溶膠等污染物增加造成霧日數(shù)增加(尤其是冬季)而導(dǎo)致能見度降低所致,這與孔鋒[35]等研究結(jié)論一致。 本文分析了昌吉地區(qū)日照時數(shù)時空變化規(guī)律特征并初步分析了日照時數(shù)與氣象因子的關(guān)系, 研究結(jié)果對明確區(qū)域氣候變化特征及其驅(qū)動因子具有一定的科學(xué)指導(dǎo)意義。 影響日照時數(shù)變化的因素是多方位和復(fù)雜的, 隨著階段和區(qū)域的不同,其影響因素也存在階段性和區(qū)域性。本文只對各氣象因素對日照時數(shù)的影響進(jìn)行了初步研究, 并未具體分析大氣污染等其它環(huán)境因素的影響,后期結(jié)合氣象因素進(jìn)一步考慮大氣環(huán)流、大氣污染特征、城市發(fā)展、人類活動和地形等因素對日照時數(shù)的影響,日照時數(shù)對農(nóng)業(yè)生產(chǎn)的影響方面的分析。

    4 結(jié)論

    (1)昌吉地區(qū)日照時數(shù)年際波動較大,且呈明顯的減少趨勢,氣候傾向率為-62.8 h/10 a,昌吉地區(qū)各季日照時數(shù)變化不一, 不同季的日照時數(shù)均呈下降趨勢, 其中冬季減少趨勢最為顯著, 氣候傾向率為-31.7 h/10 a, 且對全年日照時數(shù)的減少影響最大。

    (2)昌吉地區(qū)日照時數(shù)年內(nèi)變化呈雙峰型,最大值出現(xiàn)在5、7、8 月,分別為312.1、308.8、302.9 h,最小值出現(xiàn)在12 月,為128.7 h,日照時數(shù)最高月份是最低月份的2.4 倍;各月日照分布不均,波動較大,1—12 月日照時數(shù)均呈減少趨勢,減少速率為-1.37~-0.15 h/10 a,其中,1 月減少最為明顯,近51 a 共減少了68.5 h,12 月次之。

    (3) 昌吉地區(qū)日照時數(shù)多年平均為2 840.3 h,東部地區(qū)的日照時數(shù)明顯高于中西部地區(qū),呈現(xiàn)“東部多、中西部少,山區(qū)多、平原少”的分布特征。 昌吉地區(qū)各縣市日照時數(shù)均呈減少趨勢,其中,阜康市日照時數(shù)減少最為明顯,減少速率為-113.4 h/10 a;其余縣市減少速率為-2.54~-77.6 h/10 a。 中西部地區(qū)的日照時數(shù)減少速率明顯大于東部地區(qū)和山區(qū)。

    (4)昌吉地區(qū)日照時數(shù)在1987 年發(fā)生突變,在2002 年超出0.05 顯著性水平臨界線,在2016 年超出0.01 顯著性水平臨界線, 說明昌吉地區(qū)日照時數(shù)減少十分明顯。水汽壓的增大、霧日數(shù)的增多、大風(fēng)日數(shù)的減少是造成昌吉地區(qū)日照時數(shù)減少的主要原因。

    猜你喜歡
    昌吉日照時數(shù)日數(shù)
    適宜在昌吉春麥區(qū)種植的早熟高產(chǎn)春小麥品種篩選
    福州市近70年日照變化趨勢分析
    漢江上游漢中區(qū)域不同等級降水日數(shù)的氣候變化特征分析
    綠色科技(2022年16期)2022-09-15 03:04:46
    西昌近60年日照時數(shù)的變化特征分析
    1961~2020年曲麻萊縣日照時數(shù)變化特征
    天津市濱海新區(qū)塘沽地域雷暴日數(shù)變化規(guī)律及特征分析
    天津科技(2020年2期)2020-03-03 05:09:48
    以十九大精神為指引 展現(xiàn)新作為新氣象,開創(chuàng)昌吉學(xué)院發(fā)展新局面
    1980年~2017年大冶市日照時數(shù)變化特征分析
    在昌吉,我們品嘗到了豐收的味道——新疆昌吉漢和7S店無人機(jī)飛防作業(yè)小記
    海南省雷暴日數(shù)年代際變化特征
    最好的美女福利视频网| 午夜福利在线观看吧| 色哟哟·www| 身体一侧抽搐| 高清午夜精品一区二区三区 | 女的被弄到高潮叫床怎么办| 国产欧美日韩一区二区精品| 男插女下体视频免费在线播放| 日韩av在线大香蕉| 亚洲五月天丁香| 久久午夜亚洲精品久久| 成人无遮挡网站| 一区二区三区免费毛片| 狂野欧美激情性xxxx在线观看| 99热只有精品国产| 免费不卡的大黄色大毛片视频在线观看 | 精品久久久久久久久久久久久| 亚洲图色成人| 如何舔出高潮| 日本撒尿小便嘘嘘汇集6| 久久99热6这里只有精品| 少妇熟女欧美另类| 简卡轻食公司| 简卡轻食公司| 亚洲av中文字字幕乱码综合| 老司机福利观看| 高清毛片免费看| 日本黄色片子视频| 亚洲成a人片在线一区二区| 亚洲国产精品成人久久小说 | 国内精品宾馆在线| 99久久无色码亚洲精品果冻| 国产成人a∨麻豆精品| 天美传媒精品一区二区| 亚洲性夜色夜夜综合| 欧美激情在线99| 成人午夜高清在线视频| 欧美在线一区亚洲| 嫩草影院新地址| 最近视频中文字幕2019在线8| 国产精品三级大全| 在线免费观看的www视频| 国产伦精品一区二区三区视频9| 99国产极品粉嫩在线观看| 午夜影院日韩av| 国产精品久久久久久亚洲av鲁大| 午夜福利在线在线| 国产免费一级a男人的天堂| 伦精品一区二区三区| 99久久精品一区二区三区| 联通29元200g的流量卡| 一本精品99久久精品77| 最近最新中文字幕大全电影3| 国产在视频线在精品| 精品午夜福利视频在线观看一区| 欧美成人精品欧美一级黄| 国产黄色视频一区二区在线观看 | 我的老师免费观看完整版| 国产一级毛片七仙女欲春2| 亚洲欧美精品自产自拍| h日本视频在线播放| 欧美人与善性xxx| 免费一级毛片在线播放高清视频| 精品一区二区三区av网在线观看| 国产色婷婷99| 日韩欧美国产在线观看| 最近最新中文字幕大全电影3| 日韩欧美国产在线观看| 国产成人a区在线观看| 成熟少妇高潮喷水视频| 成人亚洲精品av一区二区| 51国产日韩欧美| 国产成人影院久久av| 一本久久中文字幕| 给我免费播放毛片高清在线观看| 禁无遮挡网站| 麻豆成人午夜福利视频| 午夜爱爱视频在线播放| 欧美成人一区二区免费高清观看| 毛片一级片免费看久久久久| 我要搜黄色片| 亚洲av第一区精品v没综合| 丝袜美腿在线中文| 丝袜美腿在线中文| 岛国在线免费视频观看| 久久6这里有精品| 淫秽高清视频在线观看| 久久6这里有精品| 精品免费久久久久久久清纯| 日韩成人av中文字幕在线观看 | a级毛色黄片| 久久久国产成人免费| 日韩 亚洲 欧美在线| 国产淫片久久久久久久久| 好男人在线观看高清免费视频| 精品一区二区三区视频在线| 国产片特级美女逼逼视频| 日韩高清综合在线| 亚洲va在线va天堂va国产| 亚洲精品久久国产高清桃花| 亚洲四区av| 在线免费十八禁| 国产精品美女特级片免费视频播放器| 麻豆国产av国片精品| 韩国av在线不卡| 亚洲成a人片在线一区二区| 欧美另类亚洲清纯唯美| 在线播放国产精品三级| 欧美潮喷喷水| 日韩成人伦理影院| 国产伦在线观看视频一区| 日韩精品中文字幕看吧| 高清午夜精品一区二区三区 | 国产精品亚洲一级av第二区| 1000部很黄的大片| 嫩草影视91久久| 亚洲精品日韩av片在线观看| 99国产极品粉嫩在线观看| 一区福利在线观看| 亚洲熟妇熟女久久| 日韩精品中文字幕看吧| 看免费成人av毛片| 两性午夜刺激爽爽歪歪视频在线观看| 寂寞人妻少妇视频99o| 91午夜精品亚洲一区二区三区| 免费看a级黄色片| 日韩欧美精品v在线| 看非洲黑人一级黄片| 亚洲综合色惰| 99久久中文字幕三级久久日本| 99热网站在线观看| 精品午夜福利在线看| 亚洲av美国av| 亚洲av一区综合| 免费人成视频x8x8入口观看| 日韩高清综合在线| 99热6这里只有精品| 亚洲精品乱码久久久v下载方式| 亚洲图色成人| 亚洲在线自拍视频| 伦精品一区二区三区| 免费av不卡在线播放| 国产精品久久久久久亚洲av鲁大| 搡女人真爽免费视频火全软件 | 久久韩国三级中文字幕| 深爱激情五月婷婷| 一个人免费在线观看电影| 欧美成人一区二区免费高清观看| 久久久色成人| 蜜桃亚洲精品一区二区三区| 日本三级黄在线观看| 国产综合懂色| 国产亚洲精品av在线| 欧美最新免费一区二区三区| 黄色日韩在线| 男女啪啪激烈高潮av片| 最好的美女福利视频网| 久久亚洲精品不卡| 久久久久久国产a免费观看| 日日摸夜夜添夜夜爱| 国产黄色小视频在线观看| 男插女下体视频免费在线播放| 麻豆av噜噜一区二区三区| 精品久久久久久成人av| 日日摸夜夜添夜夜添小说| 两性午夜刺激爽爽歪歪视频在线观看| 永久网站在线| 22中文网久久字幕| 日韩中字成人| 国产精品人妻久久久影院| 国产一区二区在线av高清观看| 欧美在线一区亚洲| 成人av在线播放网站| 国产91av在线免费观看| 久久久精品免费免费高清| 熟女av电影| 亚洲中文av在线| 色网站视频免费| 亚洲精品久久久久久婷婷小说| 黄片无遮挡物在线观看| 成人国产av品久久久| 久久久国产一区二区| 亚洲内射少妇av| 成人无遮挡网站| 国产视频内射| 91成人精品电影| 欧美97在线视频| 午夜免费男女啪啪视频观看| 高清午夜精品一区二区三区| 亚洲欧美精品专区久久| 激情五月婷婷亚洲| 啦啦啦中文免费视频观看日本| 亚洲精品国产av成人精品| 青春草国产在线视频| 十八禁高潮呻吟视频 | 黑人猛操日本美女一级片| 日日摸夜夜添夜夜爱| 日韩亚洲欧美综合| 免费看日本二区| 国产精品国产av在线观看| 黄色一级大片看看| 国产成人精品久久久久久| 亚洲精品视频女| 最新中文字幕久久久久| 色哟哟·www| 亚洲高清免费不卡视频| 午夜激情福利司机影院| 又爽又黄a免费视频| 午夜福利,免费看| 人体艺术视频欧美日本| 久久99蜜桃精品久久| 国产一区亚洲一区在线观看| 精品少妇久久久久久888优播| 99久国产av精品国产电影| 丝袜在线中文字幕| 亚洲国产精品国产精品| 人人妻人人爽人人添夜夜欢视频 | 久久精品国产自在天天线| 亚洲综合色惰| 18禁在线播放成人免费| 成年人午夜在线观看视频| 精品人妻熟女av久视频| 香蕉精品网在线| 亚洲精品一区蜜桃| 99热这里只有是精品在线观看| 如何舔出高潮| 大又大粗又爽又黄少妇毛片口| 精品少妇黑人巨大在线播放| 久久鲁丝午夜福利片| 中文天堂在线官网| 日韩中文字幕视频在线看片| 免费av中文字幕在线| 91久久精品国产一区二区成人| 久久久久国产精品人妻一区二区| 国产精品久久久久久精品古装| tube8黄色片| 精品人妻熟女毛片av久久网站| 一级毛片黄色毛片免费观看视频| 精品国产露脸久久av麻豆| 99热6这里只有精品| 亚洲精品日韩av片在线观看| 香蕉精品网在线| 汤姆久久久久久久影院中文字幕| 日韩视频在线欧美| 91精品一卡2卡3卡4卡| 欧美三级亚洲精品| 少妇人妻 视频| 偷拍熟女少妇极品色| 免费黄色在线免费观看| 日韩电影二区| 国产精品久久久久久久电影| 亚洲欧洲国产日韩| 欧美性感艳星| 欧美日韩av久久| 极品少妇高潮喷水抽搐| www.av在线官网国产| 欧美精品一区二区免费开放| 777米奇影视久久| 日本av手机在线免费观看| 亚洲av综合色区一区| 日韩视频在线欧美| 夜夜看夜夜爽夜夜摸| av免费在线看不卡| 午夜福利影视在线免费观看| 国产精品伦人一区二区| a级一级毛片免费在线观看| 午夜福利网站1000一区二区三区| 久久久久久伊人网av| 啦啦啦视频在线资源免费观看| 亚洲欧美一区二区三区国产| 免费黄频网站在线观看国产| 99国产精品免费福利视频| 亚洲精品乱码久久久久久按摩| 国产精品免费大片| 王馨瑶露胸无遮挡在线观看| 韩国av在线不卡| 又黄又爽又刺激的免费视频.| 精品国产一区二区三区久久久樱花| 久久99蜜桃精品久久| 天美传媒精品一区二区| 香蕉精品网在线| 男女无遮挡免费网站观看| 波野结衣二区三区在线| 美女主播在线视频| 久久ye,这里只有精品| 午夜老司机福利剧场| 欧美精品高潮呻吟av久久| 看免费成人av毛片| 少妇人妻一区二区三区视频| 噜噜噜噜噜久久久久久91| 国产精品一区www在线观看| 精品一品国产午夜福利视频| 亚洲电影在线观看av| 搡女人真爽免费视频火全软件| 久久6这里有精品| 中文字幕制服av| 人人妻人人爽人人添夜夜欢视频 | 搡女人真爽免费视频火全软件| 青春草亚洲视频在线观看| 一级黄片播放器| 国产伦理片在线播放av一区| 日韩一本色道免费dvd| 男女国产视频网站| 麻豆精品久久久久久蜜桃| 日韩在线高清观看一区二区三区| 蜜臀久久99精品久久宅男| √禁漫天堂资源中文www| 一区二区三区免费毛片| 内地一区二区视频在线| 久久久久久久久久久免费av| 美女国产视频在线观看| 大香蕉97超碰在线| 99热这里只有是精品50| 大话2 男鬼变身卡| 成人特级av手机在线观看| 狂野欧美激情性bbbbbb| 精品一品国产午夜福利视频| 边亲边吃奶的免费视频| 777米奇影视久久| 99久久精品国产国产毛片| 夜夜爽夜夜爽视频| 国产成人精品无人区| 国产一区亚洲一区在线观看| 在现免费观看毛片| 秋霞伦理黄片| 只有这里有精品99| 国精品久久久久久国模美| 久久精品久久久久久久性| 久久av网站| 婷婷色麻豆天堂久久| 亚洲欧洲精品一区二区精品久久久 | 久久人人爽av亚洲精品天堂| 精品一区二区三区视频在线| 国产淫语在线视频| 天天操日日干夜夜撸| 大码成人一级视频| 精品久久久精品久久久| 晚上一个人看的免费电影| 熟女电影av网| 一级片'在线观看视频| 2021少妇久久久久久久久久久| 日本av手机在线免费观看| 另类精品久久| 欧美变态另类bdsm刘玥| av在线观看视频网站免费| 最近中文字幕高清免费大全6| 国产无遮挡羞羞视频在线观看| 国产精品人妻久久久久久| 国产午夜精品久久久久久一区二区三区| 2018国产大陆天天弄谢| 欧美日韩国产mv在线观看视频| 久热久热在线精品观看| 伊人久久精品亚洲午夜| 婷婷色综合大香蕉| 久久 成人 亚洲| 国产亚洲一区二区精品| 亚洲人成网站在线观看播放| 天天躁夜夜躁狠狠久久av| 人人澡人人妻人| 日韩一区二区三区影片| 国产一区二区三区综合在线观看 | 在线免费观看不下载黄p国产| 欧美变态另类bdsm刘玥| a级毛片在线看网站| 国产一区亚洲一区在线观看| 五月天丁香电影| 大又大粗又爽又黄少妇毛片口| 亚洲精品乱码久久久久久按摩| 色视频在线一区二区三区| 韩国高清视频一区二区三区| 亚洲欧美日韩东京热| 七月丁香在线播放| 国产亚洲精品久久久com| 一级毛片aaaaaa免费看小| 欧美亚洲 丝袜 人妻 在线| 我要看日韩黄色一级片| 亚洲精品456在线播放app| 久久久国产精品麻豆| 日韩大片免费观看网站| 日韩亚洲欧美综合| 少妇的逼水好多| 91aial.com中文字幕在线观看| 我要看黄色一级片免费的| 日韩电影二区| 在线观看av片永久免费下载| 日本av手机在线免费观看| 中文在线观看免费www的网站| 日韩一本色道免费dvd| 啦啦啦中文免费视频观看日本| 免费看av在线观看网站| 十八禁网站网址无遮挡 | 亚洲精品中文字幕在线视频 | 九九久久精品国产亚洲av麻豆| 春色校园在线视频观看| a级一级毛片免费在线观看| 免费观看a级毛片全部| 99九九在线精品视频 | 色婷婷av一区二区三区视频| 国精品久久久久久国模美| 中文字幕制服av| 熟女电影av网| 搡老乐熟女国产| 在线亚洲精品国产二区图片欧美 | 精品人妻熟女毛片av久久网站| 蜜臀久久99精品久久宅男| 精品少妇久久久久久888优播| 蜜桃久久精品国产亚洲av| 男女免费视频国产| 激情五月婷婷亚洲| a级毛色黄片| 精品酒店卫生间| 成人国产av品久久久| 国产精品偷伦视频观看了| 亚洲av电影在线观看一区二区三区| 国产精品伦人一区二区| 你懂的网址亚洲精品在线观看| 少妇的逼好多水| 欧美3d第一页| 女性生殖器流出的白浆| 日韩制服骚丝袜av| 国产亚洲精品久久久com| 亚洲,欧美,日韩| 色吧在线观看| 国产乱来视频区| 老司机影院毛片| 91久久精品国产一区二区三区| 纯流量卡能插随身wifi吗| 大香蕉久久网| 国产成人精品婷婷| 婷婷色av中文字幕| 肉色欧美久久久久久久蜜桃| 精品久久国产蜜桃| 国产欧美日韩一区二区三区在线 | 亚洲第一区二区三区不卡| 日韩强制内射视频| 老女人水多毛片| av在线观看视频网站免费| 亚洲欧美清纯卡通| 中文资源天堂在线| 最近最新中文字幕免费大全7| 秋霞在线观看毛片| 亚洲成人av在线免费| 伦理电影免费视频| 特大巨黑吊av在线直播| 亚洲精品aⅴ在线观看| 卡戴珊不雅视频在线播放| 国产伦精品一区二区三区四那| 人人妻人人澡人人看| 亚洲va在线va天堂va国产| 日韩成人伦理影院| 国产一区二区在线观看日韩| 色婷婷久久久亚洲欧美| 少妇的逼水好多| 亚洲av.av天堂| 亚洲成人av在线免费| 欧美区成人在线视频| 久久久久人妻精品一区果冻| 99热6这里只有精品| 欧美三级亚洲精品| 高清不卡的av网站| 国产精品麻豆人妻色哟哟久久| 汤姆久久久久久久影院中文字幕| 热re99久久国产66热| 青春草视频在线免费观看| 成人国产麻豆网| 建设人人有责人人尽责人人享有的| 老司机影院毛片| 99九九线精品视频在线观看视频| av在线老鸭窝| 久久精品国产自在天天线| 99热全是精品| 亚洲国产精品一区二区三区在线| 青春草亚洲视频在线观看| 国产成人精品一,二区| 久久久久精品久久久久真实原创| 国产老妇伦熟女老妇高清| 精品酒店卫生间| 成人影院久久| 老熟女久久久| 日韩一本色道免费dvd| 最后的刺客免费高清国语| 一个人看视频在线观看www免费| √禁漫天堂资源中文www| 蜜桃久久精品国产亚洲av| 免费人妻精品一区二区三区视频| 精品久久久精品久久久| 日韩一区二区视频免费看| 亚洲av.av天堂| 国产精品久久久久久久久免| 午夜福利网站1000一区二区三区| 大片免费播放器 马上看| 尾随美女入室| 在线观看三级黄色| 中文字幕亚洲精品专区| 日产精品乱码卡一卡2卡三| 国产淫片久久久久久久久| 国产av码专区亚洲av| av国产精品久久久久影院| 大又大粗又爽又黄少妇毛片口| 99热国产这里只有精品6| 国内精品宾馆在线| 九九在线视频观看精品| 久久韩国三级中文字幕| 亚洲国产精品999| 国产精品人妻久久久影院| 亚洲精品一区蜜桃| 五月伊人婷婷丁香| 99热这里只有是精品在线观看| 伊人亚洲综合成人网| 好男人视频免费观看在线| 亚洲久久久国产精品| 国产亚洲欧美精品永久| 高清毛片免费看| 婷婷色综合大香蕉| 久久久久国产精品人妻一区二区| 亚洲av二区三区四区| 黄色怎么调成土黄色| 日韩视频在线欧美| 亚洲av男天堂| 在线亚洲精品国产二区图片欧美 | 日韩一本色道免费dvd| 色婷婷久久久亚洲欧美| 成人免费观看视频高清| 国产亚洲最大av| av视频免费观看在线观看| 亚洲av成人精品一二三区| 香蕉精品网在线| 男女啪啪激烈高潮av片| 国产美女午夜福利| 欧美精品一区二区大全| 亚洲国产色片| 人体艺术视频欧美日本| 精品国产国语对白av| 伊人久久精品亚洲午夜| 青青草视频在线视频观看| 国产av一区二区精品久久| 丰满人妻一区二区三区视频av| 亚洲国产精品国产精品| 精品国产乱码久久久久久小说| 久久精品夜色国产| 成人毛片a级毛片在线播放| 春色校园在线视频观看| 亚洲成人av在线免费| 美女内射精品一级片tv| 成人免费观看视频高清| 国产精品一区二区在线不卡| 美女中出高潮动态图| 美女国产视频在线观看| 狠狠精品人妻久久久久久综合| 国产午夜精品久久久久久一区二区三区| 秋霞在线观看毛片| 一级,二级,三级黄色视频| 亚洲精品国产av成人精品| 免费av不卡在线播放| 能在线免费看毛片的网站| 99精国产麻豆久久婷婷| 高清不卡的av网站| 欧美人与善性xxx| 少妇裸体淫交视频免费看高清| 99国产精品免费福利视频| av一本久久久久| 精品亚洲成a人片在线观看| 免费播放大片免费观看视频在线观看| 国产视频内射| 欧美日韩在线观看h| 国产一级毛片在线| 欧美bdsm另类| 九九爱精品视频在线观看| 五月伊人婷婷丁香| 久久久久久久大尺度免费视频| 欧美日韩在线观看h| 妹子高潮喷水视频| 久久99精品国语久久久| 国产精品伦人一区二区| 观看av在线不卡| 男人狂女人下面高潮的视频| 国产 一区精品| 嘟嘟电影网在线观看| 久久精品久久久久久噜噜老黄| 亚洲精品456在线播放app| 深夜a级毛片| 另类精品久久| 久久免费观看电影| 成年人午夜在线观看视频| 国产熟女欧美一区二区| 欧美 亚洲 国产 日韩一| 夜夜看夜夜爽夜夜摸| 老司机亚洲免费影院| 成年美女黄网站色视频大全免费 | av国产久精品久网站免费入址| 夫妻性生交免费视频一级片| 国语对白做爰xxxⅹ性视频网站| 免费观看在线日韩| 亚洲国产精品成人久久小说| 日韩一本色道免费dvd| 人妻系列 视频| 国产午夜精品久久久久久一区二区三区| 亚洲av综合色区一区| 亚洲av中文av极速乱| 麻豆精品久久久久久蜜桃| 欧美日韩在线观看h| 女的被弄到高潮叫床怎么办| 少妇丰满av| av黄色大香蕉| 日韩av免费高清视频| 国产日韩欧美在线精品| 激情五月婷婷亚洲| 97超碰精品成人国产| 在线观看美女被高潮喷水网站| 九草在线视频观看| 丰满少妇做爰视频| 久久99一区二区三区| 亚洲国产色片| 只有这里有精品99| 亚洲国产成人一精品久久久|