• <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ù)年代際變化特征
    一级片免费观看大全| 欧美大码av| 国产伦在线观看视频一区| 精品国产超薄肉色丝袜足j| 国产三级在线视频| 国产精品精品国产色婷婷| a级毛片在线看网站| 久久久久性生活片| 看片在线看免费视频| 国语自产精品视频在线第100页| 男女午夜视频在线观看| 国产精品精品国产色婷婷| 麻豆国产av国片精品| 国产成人啪精品午夜网站| 两性午夜刺激爽爽歪歪视频在线观看 | 五月伊人婷婷丁香| 亚洲七黄色美女视频| 美女 人体艺术 gogo| 91大片在线观看| 亚洲va日本ⅴa欧美va伊人久久| 日韩精品青青久久久久久| 长腿黑丝高跟| www.www免费av| 午夜成年电影在线免费观看| 在线观看免费视频日本深夜| 国产精品久久久久久精品电影| 免费观看人在逋| 一级作爱视频免费观看| 真人做人爱边吃奶动态| 97碰自拍视频| 在线观看美女被高潮喷水网站 | 国产一级毛片七仙女欲春2| 久久亚洲精品不卡| 人成视频在线观看免费观看| 欧美日韩亚洲国产一区二区在线观看| 老司机在亚洲福利影院| 在线永久观看黄色视频| 成人手机av| 在线播放国产精品三级| 欧美成人午夜精品| 亚洲av美国av| 亚洲色图 男人天堂 中文字幕| 中文在线观看免费www的网站 | 99精品久久久久人妻精品| 人妻久久中文字幕网| 国产伦一二天堂av在线观看| 久久国产乱子伦精品免费另类| 国产av不卡久久| 国产伦人伦偷精品视频| 1024视频免费在线观看| 久久久精品国产亚洲av高清涩受| 亚洲午夜理论影院| 少妇人妻一区二区三区视频| av视频在线观看入口| 日韩有码中文字幕| 国产亚洲精品综合一区在线观看 | 欧美一级毛片孕妇| 757午夜福利合集在线观看| 99国产极品粉嫩在线观看| 变态另类丝袜制服| 麻豆国产av国片精品| 黄色a级毛片大全视频| 全区人妻精品视频| 欧美黄色片欧美黄色片| 久9热在线精品视频| 亚洲人成网站高清观看| 三级国产精品欧美在线观看 | 毛片女人毛片| 怎么达到女性高潮| 精品电影一区二区在线| 国产精品爽爽va在线观看网站| 99在线视频只有这里精品首页| 九色国产91popny在线| 国产亚洲av高清不卡| 美女午夜性视频免费| 亚洲国产精品久久男人天堂| 国产视频一区二区在线看| 国产高清激情床上av| 亚洲av熟女| 99热6这里只有精品| 日日爽夜夜爽网站| 久久久久久人人人人人| 久久久久久大精品| 桃色一区二区三区在线观看| 国产亚洲精品一区二区www| 2021天堂中文幕一二区在线观| 亚洲狠狠婷婷综合久久图片| 狂野欧美激情性xxxx| 波多野结衣巨乳人妻| 亚洲 国产 在线| 久久久久久久午夜电影| 久久久国产成人精品二区| 一级毛片高清免费大全| 手机成人av网站| 亚洲av成人不卡在线观看播放网| 久久中文字幕一级| 精品午夜福利视频在线观看一区| 国产黄片美女视频| 深夜精品福利| 非洲黑人性xxxx精品又粗又长| 午夜日韩欧美国产| 午夜福利欧美成人| 中国美女看黄片| 九色国产91popny在线| 免费在线观看影片大全网站| 琪琪午夜伦伦电影理论片6080| 中文字幕精品亚洲无线码一区| 在线观看舔阴道视频| 午夜福利18| 亚洲七黄色美女视频| 久久伊人香网站| 欧美成人午夜精品| 精品久久久久久久毛片微露脸| 两人在一起打扑克的视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲片人在线观看| 麻豆成人午夜福利视频| 国产激情偷乱视频一区二区| 久久人人精品亚洲av| 18禁美女被吸乳视频| 精品久久蜜臀av无| 亚洲国产高清在线一区二区三| 日韩三级视频一区二区三区| 国产单亲对白刺激| 中文亚洲av片在线观看爽| 亚洲18禁久久av| 国产真实乱freesex| 亚洲国产欧美网| 国产在线观看jvid| 亚洲五月婷婷丁香| 9191精品国产免费久久| 日韩欧美一区二区三区在线观看| 亚洲人成电影免费在线| 欧美日韩精品网址| 深夜精品福利| 精品欧美一区二区三区在线| 国产99白浆流出| av有码第一页| 国产精品久久电影中文字幕| a级毛片在线看网站| 久久伊人香网站| 国产精品1区2区在线观看.| 黄色毛片三级朝国网站| 亚洲 欧美 日韩 在线 免费| 亚洲国产欧美网| 91av网站免费观看| 91国产中文字幕| 亚洲成a人片在线一区二区| 欧美久久黑人一区二区| 女人爽到高潮嗷嗷叫在线视频| 一二三四在线观看免费中文在| 亚洲乱码一区二区免费版| 国产一区在线观看成人免费| 日本 av在线| 51午夜福利影视在线观看| 国产熟女午夜一区二区三区| 国产成人欧美在线观看| 亚洲一区二区三区不卡视频| 久久香蕉精品热| 亚洲美女黄片视频| 久久精品国产亚洲av香蕉五月| 在线永久观看黄色视频| 亚洲av电影在线进入| 99国产精品一区二区蜜桃av| 最近最新免费中文字幕在线| av天堂在线播放| 国产爱豆传媒在线观看 | 亚洲 国产 在线| 91麻豆精品激情在线观看国产| 欧美性猛交╳xxx乱大交人| 国产一区二区激情短视频| 亚洲一码二码三码区别大吗| 日本一区二区免费在线视频| 亚洲av电影不卡..在线观看| 久久 成人 亚洲| 久久精品国产综合久久久| 午夜福利高清视频| 午夜两性在线视频| 国产高清视频在线播放一区| 人妻夜夜爽99麻豆av| 国产亚洲精品第一综合不卡| av欧美777| 久久亚洲真实| 操出白浆在线播放| 色尼玛亚洲综合影院| 黄色毛片三级朝国网站| 久久婷婷成人综合色麻豆| 日韩精品免费视频一区二区三区| 日本精品一区二区三区蜜桃| 宅男免费午夜| 啦啦啦免费观看视频1| 久久精品夜夜夜夜夜久久蜜豆 | 久久精品aⅴ一区二区三区四区| 香蕉av资源在线| 99在线人妻在线中文字幕| 亚洲九九香蕉| 午夜免费观看网址| www日本黄色视频网| 日本五十路高清| 999精品在线视频| 国产91精品成人一区二区三区| 99精品欧美一区二区三区四区| 国产精品,欧美在线| 露出奶头的视频| 九九热线精品视视频播放| 大型av网站在线播放| 久久精品成人免费网站| aaaaa片日本免费| 亚洲第一电影网av| 日韩欧美在线二视频| 一级毛片女人18水好多| 老司机深夜福利视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 妹子高潮喷水视频| 久久人人精品亚洲av| 国产精品一区二区三区四区久久| 麻豆国产97在线/欧美 | 国产精品综合久久久久久久免费| av天堂在线播放| 黄色女人牲交| 免费看美女性在线毛片视频| 午夜激情福利司机影院| 欧美久久黑人一区二区| 成人永久免费在线观看视频| 婷婷亚洲欧美| 久久婷婷成人综合色麻豆| 亚洲乱码一区二区免费版| 久久天堂一区二区三区四区| 天堂√8在线中文| 亚洲国产欧洲综合997久久,| 国产三级黄色录像| 两个人视频免费观看高清| 国产成人影院久久av| 嫁个100分男人电影在线观看| 黄色片一级片一级黄色片| 欧美成人免费av一区二区三区| 国产精品 欧美亚洲| 啦啦啦韩国在线观看视频| 欧美乱色亚洲激情| 一级毛片精品| 欧美精品啪啪一区二区三区| 国产精品综合久久久久久久免费| 欧美又色又爽又黄视频| 999久久久国产精品视频| 免费在线观看影片大全网站| 97碰自拍视频| 精华霜和精华液先用哪个| 丰满的人妻完整版| 99在线人妻在线中文字幕| 大型av网站在线播放| 亚洲成av人片免费观看| 久久精品国产亚洲av高清一级| 18美女黄网站色大片免费观看| 嫩草影院精品99| 1024手机看黄色片| 美女扒开内裤让男人捅视频| 高潮久久久久久久久久久不卡| 久久精品国产综合久久久| 天天一区二区日本电影三级| ponron亚洲| 啦啦啦免费观看视频1| 在线免费观看的www视频| 欧美性猛交╳xxx乱大交人| 久久久久九九精品影院| 国产成人av教育| 亚洲国产精品久久男人天堂| 在线国产一区二区在线| 夜夜看夜夜爽夜夜摸| 午夜影院日韩av| 国产主播在线观看一区二区| 91麻豆精品激情在线观看国产| 琪琪午夜伦伦电影理论片6080| 这个男人来自地球电影免费观看| 欧美日韩一级在线毛片| 亚洲国产中文字幕在线视频| 国产高清videossex| 免费在线观看日本一区| 国产免费男女视频| 999精品在线视频| 国产蜜桃级精品一区二区三区| 老汉色av国产亚洲站长工具| 欧美性猛交╳xxx乱大交人| 丝袜人妻中文字幕| 老司机午夜福利在线观看视频| 免费一级毛片在线播放高清视频| 亚洲午夜精品一区,二区,三区| 午夜视频精品福利| 久久欧美精品欧美久久欧美| 99国产极品粉嫩在线观看| 亚洲精品粉嫩美女一区| 国产精品 国内视频| 免费在线观看亚洲国产| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产精品合色在线| 人妻夜夜爽99麻豆av| 日韩成人在线观看一区二区三区| 久久久久久久久久黄片| 男女视频在线观看网站免费 | 级片在线观看| 岛国在线观看网站| 成年版毛片免费区| 亚洲成人精品中文字幕电影| 亚洲av美国av| 香蕉丝袜av| 搡老岳熟女国产| 亚洲无线在线观看| 在线观看午夜福利视频| 成人18禁在线播放| 国产免费男女视频| 麻豆久久精品国产亚洲av| 午夜福利欧美成人| 琪琪午夜伦伦电影理论片6080| 1024视频免费在线观看| 久久久久性生活片| 麻豆一二三区av精品| 人人妻人人看人人澡| 午夜福利高清视频| 长腿黑丝高跟| 日本一区二区免费在线视频| 99热这里只有精品一区 | cao死你这个sao货| 亚洲成人久久爱视频| 国产成人av激情在线播放| 亚洲专区国产一区二区| 一a级毛片在线观看| 国产三级黄色录像| 十八禁人妻一区二区| 首页视频小说图片口味搜索| 国产亚洲精品av在线| 亚洲成人中文字幕在线播放| 舔av片在线| 亚洲色图 男人天堂 中文字幕| 99热这里只有精品一区 | 免费搜索国产男女视频| 99riav亚洲国产免费| 成人三级做爰电影| 亚洲成av人片在线播放无| 精品电影一区二区在线| av有码第一页| 国产成人影院久久av| 国产人伦9x9x在线观看| 天天一区二区日本电影三级| 女人高潮潮喷娇喘18禁视频| 中文字幕人妻丝袜一区二区| 欧美zozozo另类| 成人国语在线视频| 国产精品亚洲美女久久久| 欧美日本亚洲视频在线播放| 亚洲18禁久久av| 日日干狠狠操夜夜爽| 欧美色欧美亚洲另类二区| 亚洲av中文字字幕乱码综合| 变态另类成人亚洲欧美熟女| 精品久久久久久,| 亚洲狠狠婷婷综合久久图片| 国产成人aa在线观看| 看黄色毛片网站| 身体一侧抽搐| 日韩精品免费视频一区二区三区| 琪琪午夜伦伦电影理论片6080| 高清毛片免费观看视频网站| 免费看日本二区| 99久久精品国产亚洲精品| 精品久久久久久久人妻蜜臀av| 91字幕亚洲| av有码第一页| 亚洲成人久久爱视频| 99精品在免费线老司机午夜| 久久婷婷成人综合色麻豆| 五月伊人婷婷丁香| www日本黄色视频网| 国产高清有码在线观看视频 | 亚洲aⅴ乱码一区二区在线播放 | 婷婷亚洲欧美| 此物有八面人人有两片| 女生性感内裤真人,穿戴方法视频| 亚洲精品国产精品久久久不卡| 午夜免费激情av| 蜜桃久久精品国产亚洲av| 精品福利观看| 成人国产一区最新在线观看| 19禁男女啪啪无遮挡网站| 国产成人一区二区三区免费视频网站| 国产在线精品亚洲第一网站| 变态另类成人亚洲欧美熟女| 中文字幕人成人乱码亚洲影| 欧美午夜高清在线| 极品教师在线免费播放| 国产日本99.免费观看| 香蕉久久夜色| 黄色丝袜av网址大全| 欧美日韩国产亚洲二区| 丁香欧美五月| 毛片女人毛片| 一级毛片女人18水好多| 一进一出抽搐gif免费好疼| 国产一区二区三区在线臀色熟女| 日本三级黄在线观看| 精品国产超薄肉色丝袜足j| 国产黄a三级三级三级人| 国产成人av教育| 97超级碰碰碰精品色视频在线观看| 欧美日韩国产亚洲二区| 久久国产乱子伦精品免费另类| 午夜福利欧美成人| 国产成人精品久久二区二区91| 狂野欧美激情性xxxx| 搡老熟女国产l中国老女人| 婷婷精品国产亚洲av| 国产91精品成人一区二区三区| 国产成人影院久久av| 最新美女视频免费是黄的| 老司机在亚洲福利影院| 99re在线观看精品视频| 亚洲中文av在线| 制服诱惑二区| 国产探花在线观看一区二区| 亚洲专区中文字幕在线| 啦啦啦免费观看视频1| 女警被强在线播放| 50天的宝宝边吃奶边哭怎么回事| 国产黄色小视频在线观看| 热99re8久久精品国产| 欧美绝顶高潮抽搐喷水| 国内揄拍国产精品人妻在线| 午夜两性在线视频| 日韩欧美一区二区三区在线观看| 一边摸一边抽搐一进一小说| 国产精品 国内视频| 黄色丝袜av网址大全| 久久香蕉精品热| 国产三级中文精品| 舔av片在线| 日本一区二区免费在线视频| 天堂影院成人在线观看| 亚洲中文字幕日韩| 亚洲精品粉嫩美女一区| 亚洲精品中文字幕一二三四区| 99久久综合精品五月天人人| 国产精品一区二区三区四区免费观看 | 午夜福利欧美成人| www.www免费av| 黄色毛片三级朝国网站| 亚洲专区国产一区二区| 欧美午夜高清在线| 俺也久久电影网| 制服诱惑二区| 男男h啪啪无遮挡| 一a级毛片在线观看| 这个男人来自地球电影免费观看| 男女那种视频在线观看| 亚洲欧美日韩东京热| 午夜福利18| 久久天躁狠狠躁夜夜2o2o| 精品一区二区三区四区五区乱码| 99精品在免费线老司机午夜| 此物有八面人人有两片| 丝袜美腿诱惑在线| 久久久久九九精品影院| 一进一出抽搐动态| 午夜免费激情av| 欧洲精品卡2卡3卡4卡5卡区| 婷婷亚洲欧美| 亚洲全国av大片| 亚洲精品美女久久久久99蜜臀| aaaaa片日本免费| 成人精品一区二区免费| 女人高潮潮喷娇喘18禁视频| 十八禁人妻一区二区| 免费观看人在逋| 高清在线国产一区| 久久精品91蜜桃| 国产精品九九99| 国产v大片淫在线免费观看| 欧美精品啪啪一区二区三区| 制服人妻中文乱码| 国产精品一区二区三区四区免费观看 | 看片在线看免费视频| 99国产极品粉嫩在线观看| 日韩欧美免费精品| 国产成人啪精品午夜网站| 婷婷精品国产亚洲av| 久久精品夜夜夜夜夜久久蜜豆 | 久久精品综合一区二区三区| 一本久久中文字幕| 婷婷亚洲欧美| www.999成人在线观看| 国产成人啪精品午夜网站| 午夜视频精品福利| 九九热线精品视视频播放| 在线播放国产精品三级| 动漫黄色视频在线观看| 亚洲国产精品999在线| tocl精华| 国产精品一区二区三区四区免费观看 | 中出人妻视频一区二区| 欧美成人午夜精品| 免费人成视频x8x8入口观看| 91大片在线观看| 日韩精品青青久久久久久| a级毛片在线看网站| 高潮久久久久久久久久久不卡| 国产黄片美女视频| 久久精品综合一区二区三区| 91大片在线观看| 亚洲精品一区av在线观看| 毛片女人毛片| 国产探花在线观看一区二区| 日韩免费av在线播放| 亚洲狠狠婷婷综合久久图片| 夜夜夜夜夜久久久久| 日韩国内少妇激情av| 美女午夜性视频免费| 久久久国产成人免费| 91在线观看av| 国产亚洲精品久久久久5区| 好看av亚洲va欧美ⅴa在| 亚洲男人天堂网一区| 在线播放国产精品三级| 国产一区在线观看成人免费| 成人av在线播放网站| 亚洲中文字幕日韩| 国产精品爽爽va在线观看网站| 国产真实乱freesex| 久久亚洲真实| www日本在线高清视频| 日韩精品中文字幕看吧| 最近最新免费中文字幕在线| 欧美三级亚洲精品| 99精品久久久久人妻精品| 变态另类成人亚洲欧美熟女| 亚洲五月婷婷丁香| 999精品在线视频| 99久久国产精品久久久| 日本黄大片高清| tocl精华| 亚洲一区二区三区色噜噜| 18美女黄网站色大片免费观看| 黄色成人免费大全| 舔av片在线| 香蕉丝袜av| 午夜成年电影在线免费观看| 日日摸夜夜添夜夜添小说| 国产精品1区2区在线观看.| 国产一区二区激情短视频| 一卡2卡三卡四卡精品乱码亚洲| 两性午夜刺激爽爽歪歪视频在线观看 | 久久精品综合一区二区三区| bbb黄色大片| 真人做人爱边吃奶动态| 亚洲av五月六月丁香网| 亚洲av日韩精品久久久久久密| 精品久久久久久久人妻蜜臀av| 又粗又爽又猛毛片免费看| 国产一区二区激情短视频| 国产激情偷乱视频一区二区| 午夜福利18| 成人手机av| 中文资源天堂在线| 露出奶头的视频| 人妻丰满熟妇av一区二区三区| 18禁观看日本| 免费看十八禁软件| 淫妇啪啪啪对白视频| 91字幕亚洲| 国产亚洲欧美98| 久久婷婷成人综合色麻豆| 最近最新免费中文字幕在线| 女人爽到高潮嗷嗷叫在线视频| 一进一出好大好爽视频| 亚洲第一电影网av| 国产在线观看jvid| 操出白浆在线播放| 亚洲中文日韩欧美视频| 国产蜜桃级精品一区二区三区| 日本一二三区视频观看| 欧美乱码精品一区二区三区| 男女床上黄色一级片免费看| 免费在线观看亚洲国产| 日韩精品中文字幕看吧| АⅤ资源中文在线天堂| 精品第一国产精品| 91麻豆精品激情在线观看国产| 精品午夜福利视频在线观看一区| 精品高清国产在线一区| 亚洲电影在线观看av| av福利片在线观看| 久久中文字幕一级| 国产熟女xx| 丁香欧美五月| 国产欧美日韩一区二区精品| 欧美最黄视频在线播放免费| 欧美日韩瑟瑟在线播放| 国产高清视频在线观看网站| 无人区码免费观看不卡| 久久久久国产精品人妻aⅴ院| 精品久久久久久久毛片微露脸| 国产片内射在线| 不卡av一区二区三区| 国产蜜桃级精品一区二区三区| 国产精品亚洲av一区麻豆| 黄色 视频免费看| 久久久久国内视频| 成人国产综合亚洲| 亚洲av电影在线进入| 免费观看精品视频网站| 国产爱豆传媒在线观看 | 日本一二三区视频观看| 91九色精品人成在线观看| 国产精品久久久人人做人人爽| 欧美成人一区二区免费高清观看 | 精品国产乱码久久久久久男人| 黄色毛片三级朝国网站| 19禁男女啪啪无遮挡网站| 特大巨黑吊av在线直播|