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

    1964-2017年山西省潛在蒸發(fā)量時空變化及其影響因素分析

    2019-09-05 08:47:16孫從建鄭振婧侯慧新
    水土保持研究 2019年5期
    關(guān)鍵詞:蒸發(fā)量站點山西省

    孫從建, 鄭振婧, 李 偉, 楊 偉, 侯慧新

    (1.山西師范大學(xué) 地理科學(xué)學(xué)院, 山西 臨汾 041000;2.中國科學(xué)院 新疆生態(tài)與地理研究所 荒漠與綠洲生態(tài)國家重點實驗室, 烏魯木齊 830011)

    潛在蒸發(fā)量是指水分充足條件下蒸散發(fā)到空氣中的水量,又稱可能蒸發(fā)量或蒸發(fā)能力[1]。作為水循環(huán)過程的一個重要參數(shù),對地表降水及能量輻射的再次分配產(chǎn)生極大的影響[1-2]。因而,近幾十年來有關(guān)潛在蒸發(fā)量研究也在不斷的增多。潛在蒸發(fā)量的計算模型有基于氣溫的Thornthwaite方法、基于太陽輻射的Makkink方法、Penman-Monteith模型及蒸發(fā)皿蒸發(fā)等[3],但是Penman-Monteith模型因其理論基礎(chǔ)較為完整且準(zhǔn)確性較高,在世界各地被廣泛應(yīng)用,如國外的佛羅里達(dá)州南部[4]、伊朗[5]、加拿大[6]、澳大利亞[7]、朝鮮[8]等,我國的黃河上游瑪曲地區(qū)[9]、天山北坡[10]、祁連山及河西走廊[11]、涇惠渠灌區(qū)[12]、三江源區(qū)[13]等均運用Penman-Monteith模型計算潛在蒸發(fā)量并進行研究。

    由于潛在蒸發(fā)量的影響因子眾多,其變化成因十分復(fù)雜,不同地區(qū)潛在蒸發(fā)量的時空分布特征及其影響因子也不盡相同。劉普幸等[14]利用甘肅省27個氣象站點的逐日氣象資料,分析了潛在蒸發(fā)量的時空變化及影響因子,研究得出甘肅省潛在蒸發(fā)量除河西走廊外均呈上升趨勢,在空間分布上自東南向西北遞減,風(fēng)速與太陽輻射是其主導(dǎo)因素;劉聞等[15]研究得出渭河流域關(guān)中段潛在蒸發(fā)量隨海拔的降低而增高,夏季潛在蒸發(fā)量的變化趨勢與年變化趨勢一致,潛在蒸發(fā)量與日較差、平均氣溫、平均風(fēng)速、日照時數(shù)呈正相關(guān),與相對濕度和水汽壓呈負(fù)相關(guān);朱國鋒等[16]利用橫斷山區(qū)20個氣象站點的資料分析得出橫斷山區(qū)潛在蒸發(fā)量在研究時段呈先減少后上升的趨勢,氣溫、風(fēng)速及日照時數(shù)是影響潛在蒸發(fā)量的主導(dǎo)因素。自20世紀(jì)80年代以來,有很多學(xué)者做了關(guān)于潛在蒸發(fā)量的研究[17],但是目前有關(guān)山西省潛在蒸發(fā)量及其影響因素的研究較少,對于影響潛在蒸發(fā)量的影響因子分析更是較為匱乏。

    山西省地處黃土高原東緣,地貌類型復(fù)雜[18],不同地形條件對氣象因子影響較大,并且同時存在水資源緊缺等問題[19],因而研究山西省潛在蒸發(fā)量的變化及其影響因素,有助于合理規(guī)劃水資源,為農(nóng)業(yè)發(fā)展提供科學(xué)依據(jù)。鑒于此,本文基于山西省17個氣象站點1964—2017年逐日氣象數(shù)據(jù),采用Penman-Monteith模型計算出山西省潛在蒸發(fā)量,并且運用小波分析、Mann-Kendall非參數(shù)檢驗和地統(tǒng)計等方法,分析山西省近53 a潛在蒸發(fā)量時空分布特征,運用主成分分析法分析不同氣象因素及地理因素對潛在蒸發(fā)量的影響,為山西省水資源永續(xù)利用及農(nóng)業(yè)健康發(fā)展提供重要的科學(xué)依據(jù)。

    1 研究區(qū)概況

    本文選取黃土高原東部的山西省作為典型研究區(qū),涉及大同、右玉、河曲、五臺、五寨、興縣、原平、太原、離石、介休、榆社、隰縣、臨汾、長治、侯馬、陽城及運城17個氣象站點(圖1)。山西省位于34°34′—40°43′N,110°14′—114°34′E。平均海拔高度1 000 m左右,最低海拔202 m,最高海拔達(dá)3 059 m。研究區(qū)地勢東北部高西南部低,地貌類型山地、丘陵、臺地、高原及盆地均有分布[20]。氣候?qū)儆跍貛Т箨懶约撅L(fēng)氣候,夏季高溫多雨,冬季寒冷干燥。年均降水量400~600 mm,年均溫3~14℃。河流較多,擁有黃河水系和海河水系。土壤類型主要為褐土。植被類型多樣,主要以草地、林地和耕地為主[18]。

    圖1 山西省各氣象站點分布

    2 數(shù)據(jù)與方法

    2.1 數(shù)據(jù)來源

    選取山西省空間分布均勻、數(shù)據(jù)長度連續(xù)的17個氣象站點,收集各站點1964—2017年的風(fēng)速(m/s),日均溫(℃),日最低氣溫(℃),日最高氣溫(℃),日照時數(shù)(h),相對濕度(%)及平均水氣壓(KPa)數(shù)據(jù),計算月、季、年潛在蒸發(fā)量,并通過算術(shù)平均計算每個站點多年平均潛在蒸發(fā)量。此外,還需收集17個站點的經(jīng)緯度和高程數(shù)據(jù)。本文季節(jié)劃分運用氣象季節(jié),即春季3—5月份,夏季6—8月份,秋季9—11月份,冬季12月份—次年2月份。氣象數(shù)據(jù)收錄于中國氣象局國家氣象信息中心(http:∥data.cma.cn/site/index.html)。

    2.2 研究方法

    2.2.1 Penman-Monteith模型 采用1998年世界糧農(nóng)組織(FAO)修正后的Penman-Monteith模型計算潛在蒸發(fā)量。其模型如下:

    (1)

    式中:ET0表示潛在蒸發(fā)量(mm/d);Δ代表飽和水汽壓曲線斜率(kPa/℃) ;Rn是太陽凈輻射[MJ/(m2·d)];G為土壤熱通量[MJ/(m2·d)];γ為干濕常數(shù)(kPa/℃);u2表示離地面2 m高處的風(fēng)速(m/s);T為平均氣溫(℃);es是飽和水汽壓平均值;ea代表實際水汽壓;(es-ea)是飽和水汽壓差值(kPa)。

    2.2.2 時空特征分析方法 潛在蒸發(fā)量時間變化主要采用趨勢分析及周期分析方法,空間變化主要采用地統(tǒng)計插值。由于Mann-Kendall檢驗不受少數(shù)異常值的干擾,樣本分布不受一定條件的限制,在水文、氣象等非正態(tài)分布的數(shù)據(jù)資料分析中得到了廣泛的應(yīng)用[21]。通過Mann-Kendall統(tǒng)計量Z值判斷序列數(shù)據(jù)的變化趨勢,當(dāng)Z>0時表示呈上升趨勢,Z<0時表示下降趨勢,Z的絕對值越大,說明該序列的變化趨勢越顯著,當(dāng)|Z|>1.96,說明通過0.05顯著性檢驗;當(dāng)|Z|>2.58,說明通過0.01顯著性檢驗。采用小波分析方法(Wavelet Analysis)[22]分析潛在蒸發(fā)量的周期性變化規(guī)律。采用ArcGIS地統(tǒng)計模塊的反距離權(quán)重插值法[23],分析年潛在蒸發(fā)量的空間分布。

    2.2.3 主成分分析法 采用主成分分析法對山西省潛在蒸發(fā)量的影響因素進行分析。主成分分析是把多個指標(biāo)化為少數(shù)幾個綜合指標(biāo)的一種統(tǒng)計分析方法。主要是把給定的一組相關(guān)變量通過線性變換轉(zhuǎn)換成另一組不相關(guān)的變量,新變量按方差依次遞減排列,數(shù)學(xué)變換過程保持變量的總方差不變,第1個變量方差最大,稱為第一主成分,第2個變量方差其次,且和第1個變量不相關(guān),稱為第二主成分,依次類推[24]。降維產(chǎn)生的新變量在不損失原有信息的情況下,使原有變量所代表的信息更集中、更典型地體現(xiàn)出來[25]。

    3 結(jié)果與分析

    3.1 潛在蒸發(fā)量動態(tài)變化特征

    3.1.1 潛在蒸發(fā)量的年際變化 山西省1964—2017年平均潛在蒸發(fā)量為1 148.6 mm,年際傾向斜率為0.3 mm/10 a,總體呈現(xiàn)微弱下降趨勢(圖2)。1964—2017年潛在蒸發(fā)量表現(xiàn)為波動變化趨勢,其最值分別出現(xiàn)在1965年和2016年。2016年潛在蒸發(fā)量達(dá)到最小值,最小蒸發(fā)量為997.7 mm。1965年潛在蒸發(fā)量達(dá)到1 234.1 mm,為研究時段最大值。1965—1996年潛在蒸發(fā)量處于波動下降趨勢,且保持在較低水平。1996—2001年潛在蒸發(fā)量出現(xiàn)大幅度提高,但是2001—2003年又出現(xiàn)大幅下降,3 a間潛在蒸發(fā)量下降了117.94 mm,隨后進入小幅波動上升期,2013年后潛在蒸發(fā)量急劇下降。

    圖2 1964-2017年山西省年潛在蒸發(fā)量的變化

    采用小波周期分析法對山西省53 a潛在蒸發(fā)量序列進行分析(圖3),發(fā)現(xiàn)山西省潛在蒸發(fā)量存在27 a左右的主周期變化以及9 a,45 a左右的次周期變化。在27 a左右主周期變化上,有2次大的變動,經(jīng)歷了“-,+”相位1.5次循環(huán)。9 a及45 a左右為小周期變化,其周期分布比較明顯,在整個研究時段均有表現(xiàn)。

    3.1.2 潛在蒸發(fā)量季節(jié)變化 山西省潛在蒸發(fā)量季節(jié)變化趨勢如圖4所示。春季及冬季呈現(xiàn)增長趨勢,夏季及秋季呈現(xiàn)下降趨勢。山西省春季平均潛在蒸發(fā)量為321.6 mm,占全年潛在蒸發(fā)量的27.9%,1964—2017年春季潛在蒸發(fā)量的年際傾向斜率為4.2 mm/10 a,呈現(xiàn)顯著增長趨勢。夏季潛在蒸發(fā)量呈現(xiàn)下降趨勢,每10 a下降6.4 mm,平均潛在蒸發(fā)量為483.2 mm,占全年潛在蒸發(fā)量的42.1%。秋季平均潛在蒸發(fā)量為242.5 mm,占全年潛在蒸發(fā)量的21.1%,年際傾向斜率為0.8 mm/10 a,在研究時段小幅下降。冬季潛在蒸發(fā)量整體呈現(xiàn)上升趨勢,每10 a增加2.4 mm,平均潛在蒸發(fā)量為102.1 mm,占全年潛在蒸發(fā)量的8.9%。

    3.2 潛在蒸發(fā)量空間變化特征

    3.2.1 潛在蒸發(fā)量空間分布規(guī)律 根據(jù)Penman-Monteith模型計算出山西省17個氣象站點1964—2017年53 a的潛在蒸發(fā)量數(shù)據(jù)進行插值得到研究區(qū)潛在蒸發(fā)量的空間分布情況(圖5)。山西省多年平均潛在蒸發(fā)量空間分布特征表現(xiàn)為由中部向南北兩邊遞減,北部地區(qū)潛在蒸發(fā)量低于南部地區(qū)。潛在蒸發(fā)量由中部太原盆地地區(qū)的3 759.1 mm減少至北部五臺山地區(qū)的535.4 mm。

    圖3 基于小波分析年潛在蒸發(fā)量動態(tài)變化

    圖4 1964-2017年山西省不同季節(jié)潛在蒸發(fā)量變化

    3.2.2 潛在蒸發(fā)量空間變化趨勢分析 通過統(tǒng)計每一個氣象站1964—2017年潛在蒸發(fā)量Mann-Kendall統(tǒng)計量Z值,進一步分析山西省潛在蒸發(fā)量的變化趨勢,基于ArcGIS繪制出潛在蒸發(fā)量Z值的空間分布圖(圖6)。

    山西省潛在蒸發(fā)量的Mann-KendallZ值呈現(xiàn)出明顯的地域差異,西北部和東南部潛在蒸發(fā)量呈現(xiàn)下降趨勢,東北部及西南部呈增加趨勢,東北部增加趨勢尤為顯著,通過了0.01顯著性檢驗。在東北部7個站點中,有5個站點的潛在蒸發(fā)量表現(xiàn)為顯著增加趨勢,分別以2.89,3.56,3.37,3.04,5.73 mm/a的速率在增加。

    不同季節(jié)研究區(qū)潛在蒸發(fā)量變化也存在顯著差異(圖7),春季大部分地區(qū)潛在蒸發(fā)量呈增加趨勢,北部地區(qū)呈現(xiàn)顯著增加趨勢,有5個站點通過了顯著性檢驗(|Z|>2.58)。夏季潛在蒸發(fā)量變化大部分呈現(xiàn)減少趨勢,有9個站點通過了0.05顯著性檢驗,其中站點榆社通過了0.01顯著性檢驗。東北部地區(qū)出現(xiàn)增長趨勢,并且通過了顯著性檢驗(|Z|>1.96)。秋季除了站點河曲與大同呈現(xiàn)減少趨勢外,北部大部分地區(qū)呈現(xiàn)增加趨勢,并且有4個站點通過了0.05顯著性檢驗,南部大部分地區(qū)呈現(xiàn)減少趨勢,但是并未通過0.05顯著性檢驗。冬季整體呈現(xiàn)增加趨勢,有7個站點通過了顯著性檢驗(|Z|>2.58)。4個季節(jié)總體都表現(xiàn)為站點五臺地區(qū)呈現(xiàn)顯著增加趨勢,分別以5.86,5.13,5.69,6.89 mm/a的速率增加。

    圖5 山西省年潛在蒸發(fā)量空間分布

    圖6 潛在蒸發(fā)量Mann-KendallZ值的空間分布

    3.3 潛在蒸發(fā)量變化的成因分析

    影響潛在蒸發(fā)量的因素眾多,涉及到氣象因子、地理因子等[13]。本文在利用Penman-Monteith模型計算潛在蒸發(fā)量時用到的氣象因子有平均氣溫(Tave)、最高氣溫(Tmax)、最低氣溫(Tmin)、相對濕度(RH)、平均水氣壓(Ea)、風(fēng)速(W)、日照時數(shù)(SH),地理因子包括緯度(LAT)和海拔高度(H)。分析以上因子對山西省年潛在蒸發(fā)量的影響程度。主成分識別是以各因子作為原變量,通過計算變量方差和協(xié)方差矩陣的特征值,將多個變量通過降維轉(zhuǎn)化為少數(shù)幾個綜合變量,即對影響因子信息進行集中和提取,從中識別出起主導(dǎo)作用的成分[26]。表1為各影響因子解釋的總方差,從中可以看出前4個主成分的累計貢獻(xiàn)率達(dá)89.03%,說明前4個主成分已經(jīng)將9個因子所涵蓋的信息集中的體現(xiàn)了出來,選用前4個主成分代替原先的潛在蒸發(fā)量的影響因子。在前4個主成分中,第1主成分所涵蓋的信息最多,達(dá)到了33.56%。

    圖7 1964-2017年山西省不同季節(jié)潛在蒸發(fā)量Z值的空間變化

    表2為通過最大方差法得到旋轉(zhuǎn)后的因子載荷矩陣。因子載荷矩陣中絕對值的大小決定了其載荷量的多少,表明對該主成分的影響高低程度[27]。從表中可以看出,相對濕度和平均水氣壓在第1主成分上有較大載荷,載荷分別為0.28,0.29,并且這兩個因子與第1主成分呈負(fù)相關(guān)。第2主成分有較大載荷的因子分別是平均氣溫、最高氣溫和最低氣溫,這3個因子與第2主成分呈正相關(guān)。緯度在第3主成分上載荷較大,并且呈正相關(guān)。第4主成分上載荷較大的是風(fēng)速、日照時數(shù)和海拔高度,風(fēng)速和日照時數(shù)呈正相關(guān),海拔高度呈負(fù)相關(guān)。本文通過對潛在蒸發(fā)和影響因子做主成分分析,結(jié)果表明:各因子對潛在蒸發(fā)量的影響排序依次為平均水氣壓>相對濕度>最高氣溫>最低氣溫>平均溫度>緯度>海拔高度>風(fēng)速>日照時數(shù)。

    表1 山西省潛在蒸發(fā)量的各主成分及其特征值和貢獻(xiàn)率 %

    表2 山西省潛在蒸發(fā)量的因子載荷矩陣

    4 討 論

    利用日尺度氣象資料通過Penman-Monteith模型計算出潛在蒸發(fā)量可以對地區(qū)水資源合理利用、農(nóng)業(yè)生產(chǎn)發(fā)展提供科學(xué)依據(jù)。本文利用山西省17個氣象站點逐日氣象資料,分析了山西省潛在蒸發(fā)量的時空分布特征及其變化趨勢,結(jié)果表明:山西省潛在蒸發(fā)量呈現(xiàn)波動變化,總體表現(xiàn)為年潛在蒸發(fā)量呈上升趨勢,年均潛在蒸發(fā)量的最大最小值分別為(1 234,997 mm)。對比山西省與鄰近地區(qū)的潛在蒸發(fā)量分布特征可知(表3),研究區(qū)北部的內(nèi)蒙古地區(qū)[28]、研究區(qū)南部的河南地區(qū)[29]以及研究區(qū)東部的河北地區(qū)[30]潛在蒸發(fā)量在近50 a都呈下降趨勢,其年均潛在蒸發(fā)量的最大最小值分別為(834,790 mm),(1 247,418 mm),(1 488,1 190 mm)。研究者通過對黃土高原西部甘肅省[31]潛在蒸發(fā)量進行研究,得出黃土高原西部地區(qū)潛在蒸發(fā)量近50 a整體表現(xiàn)為上升趨勢。相同的潛在蒸發(fā)量變化趨勢(上升)在黃土高原中南部的陜西地區(qū)[32]也得到驗證,其年均潛在蒸發(fā)量的最大最小值分別為(1 153,895 mm)。黃土高原西部及中南部潛在蒸發(fā)量的變化趨勢與本文研究結(jié)論類似。

    一些研究者通過對于影響區(qū)域潛在蒸發(fā)量的驅(qū)動因素的分析結(jié)果顯示[12-13,32-34]:氣溫、風(fēng)速、日照時數(shù)、相對濕度、飽和水汽壓差和云量等諸多氣象要素是區(qū)域潛在蒸發(fā)量變化的控制因素,但不同區(qū)域不同的控制因素影響權(quán)重不盡相同[1,14,16,33]。山西省潛在蒸發(fā)量總體表現(xiàn)為在五臺地區(qū)較低,在太原盆地地區(qū)較高,這一現(xiàn)象主要是受了海拔高度的影響,五臺地區(qū)有“華北屋脊”之稱,海拔高度達(dá)2 208.3 m,而太原盆地地區(qū)地勢較低。潛在蒸發(fā)量的變化除了受氣象因子、地理因素影響外還受其他因素的間接影響,例如人類的活動[32]。因此,在今后還需深入研究山西省潛在蒸發(fā)量變化的原因,更好的為水資源管理、利用及農(nóng)業(yè)發(fā)展提供科學(xué)依據(jù)。

    表3 研究區(qū)鄰近區(qū)域潛在蒸發(fā)量變化趨勢及其最值

    5 結(jié) 論

    (1) 山西省多年潛在蒸發(fā)量在不斷上升,其平均值為1 148.6 mm。春季、秋季和冬季潛在蒸發(fā)量分別出現(xiàn)不同程度的上升,春季上升幅度最大,氣候傾向率為4.2 mm/10 a。夏季潛在蒸發(fā)量在不斷下降。山西省潛在蒸發(fā)量存在27 a左右的主周期變化以及9 a,45 a左右的次周期變化特征。

    (2) 山西省潛在蒸發(fā)量空間分布特征由中部向南北兩邊遞減,整體體現(xiàn)為地勢較低的地方潛在蒸發(fā)量高。

    (3) 山西省西北部和東南部潛在蒸發(fā)量呈現(xiàn)下降趨勢,東北部及西南部呈增加趨勢,東北部增加趨勢尤為顯著,通過了0.01顯著性檢驗。春季、秋季及冬季變化趨勢與年變化趨勢一致,夏季大部分地區(qū)呈現(xiàn)減少趨勢。

    (4) 平均水氣壓和相對濕度是影響山西省潛在蒸發(fā)量的主導(dǎo)因素。潛在蒸發(fā)量與平均水氣壓、相對濕度和海拔高度呈負(fù)相關(guān),與溫度、風(fēng)速、日照時數(shù)及緯度呈正相關(guān)。

    猜你喜歡
    蒸發(fā)量站點山西省
    山西省右玉縣的生態(tài)建設(shè)及其啟示
    山西省秋收已完成93.44%
    山西省打開農(nóng)產(chǎn)品供深通道
    基于Web站點的SQL注入分析與防范
    電子制作(2019年14期)2019-08-20 05:43:42
    2017~2018年冬季西北地區(qū)某站點流感流行特征分析
    1958—2013年沽源縣蒸發(fā)量變化特征分析
    1981—2010年菏澤市定陶區(qū)蒸發(fā)量變化特征分析
    新疆民豐縣地表水面蒸發(fā)量分析
    首屆歐洲自行車共享站點協(xié)商會召開
    中國自行車(2017年1期)2017-04-16 02:53:52
    山西省大水網(wǎng)工程建設(shè)如火如荼
    山西水利(2017年1期)2017-03-07 08:50:31
    男人和女人高潮做爰伦理| 免费av观看视频| 婷婷精品国产亚洲av| 日本五十路高清| 有码 亚洲区| 禁无遮挡网站| 精品人妻一区二区三区麻豆 | 久久人妻av系列| 国内精品一区二区在线观看| 国产一区二区在线观看日韩| 国产精品女同一区二区软件 | 五月伊人婷婷丁香| 国产精品自产拍在线观看55亚洲| 午夜激情福利司机影院| 亚洲午夜理论影院| 午夜福利18| 亚洲无线在线观看| 国内精品宾馆在线| 精品人妻一区二区三区麻豆 | 国产高清视频在线播放一区| 亚洲国产欧洲综合997久久,| 午夜精品久久久久久毛片777| 老女人水多毛片| 欧美日韩综合久久久久久 | 99热这里只有是精品在线观看| 成人一区二区视频在线观看| 亚洲专区中文字幕在线| 久久精品影院6| 国产精品一及| 国产伦精品一区二区三区四那| 亚洲第一电影网av| 91精品国产九色| 国产av一区在线观看免费| 国产高清三级在线| www.www免费av| 国产精品一及| 草草在线视频免费看| 黄色女人牲交| 亚洲精华国产精华精| 亚洲av电影不卡..在线观看| 亚洲欧美日韩东京热| 午夜福利18| 一级毛片久久久久久久久女| 色av中文字幕| 波野结衣二区三区在线| 88av欧美| 午夜激情欧美在线| 搞女人的毛片| 国产蜜桃级精品一区二区三区| 亚洲av熟女| 午夜福利成人在线免费观看| 黄色日韩在线| 亚洲av美国av| 久久午夜亚洲精品久久| 赤兔流量卡办理| 欧美性猛交╳xxx乱大交人| 国产成人aa在线观看| 可以在线观看毛片的网站| 精品福利观看| 韩国av在线不卡| 久9热在线精品视频| 黄色丝袜av网址大全| 不卡视频在线观看欧美| 国内少妇人妻偷人精品xxx网站| 久久久成人免费电影| or卡值多少钱| 精品久久久久久久末码| 真人一进一出gif抽搐免费| 亚洲成a人片在线一区二区| 一级av片app| 久久久久精品国产欧美久久久| 美女被艹到高潮喷水动态| 色综合色国产| 成人国产综合亚洲| 国产熟女欧美一区二区| 精品人妻偷拍中文字幕| 久久精品人妻少妇| 深夜a级毛片| 婷婷六月久久综合丁香| 久久中文看片网| 亚洲内射少妇av| 免费看美女性在线毛片视频| 亚洲精品影视一区二区三区av| 99riav亚洲国产免费| 国产淫片久久久久久久久| 欧美3d第一页| 大又大粗又爽又黄少妇毛片口| 日本黄大片高清| 长腿黑丝高跟| 久久99热6这里只有精品| 国产精品99久久久久久久久| 国产蜜桃级精品一区二区三区| 18禁裸乳无遮挡免费网站照片| 国产大屁股一区二区在线视频| 日本黄色视频三级网站网址| 精品久久国产蜜桃| 一级黄片播放器| 亚洲午夜理论影院| 亚洲成人中文字幕在线播放| 国产精品久久视频播放| 国产av麻豆久久久久久久| 嫩草影视91久久| 国产精品精品国产色婷婷| 无人区码免费观看不卡| 欧美不卡视频在线免费观看| 香蕉av资源在线| 国产亚洲精品久久久com| 国内精品久久久久精免费| 成人毛片a级毛片在线播放| 男女边吃奶边做爰视频| 欧美丝袜亚洲另类 | 国产一区二区在线av高清观看| 熟女人妻精品中文字幕| 亚洲va日本ⅴa欧美va伊人久久| eeuss影院久久| 国产精品久久视频播放| 中文字幕av在线有码专区| 国产男靠女视频免费网站| 日本爱情动作片www.在线观看 | 国产白丝娇喘喷水9色精品| 久久99热这里只有精品18| 一本久久中文字幕| 中文资源天堂在线| 国产精品av视频在线免费观看| 精品午夜福利在线看| 欧美日韩中文字幕国产精品一区二区三区| 国产69精品久久久久777片| 免费看光身美女| 国产主播在线观看一区二区| 久久午夜亚洲精品久久| 中文字幕高清在线视频| 国产精品av视频在线免费观看| 国产亚洲av嫩草精品影院| 国产欧美日韩精品一区二区| 国内揄拍国产精品人妻在线| 欧美xxxx性猛交bbbb| 尾随美女入室| 99久久中文字幕三级久久日本| 亚洲欧美日韩高清专用| 非洲黑人性xxxx精品又粗又长| 人人妻人人看人人澡| 精品久久久久久,| 欧美在线一区亚洲| 99久久精品一区二区三区| 亚洲成人精品中文字幕电影| 给我免费播放毛片高清在线观看| 国产精品亚洲美女久久久| 国内揄拍国产精品人妻在线| 少妇人妻精品综合一区二区 | 成熟少妇高潮喷水视频| 亚洲一区二区三区色噜噜| 九九热线精品视视频播放| 男人狂女人下面高潮的视频| 色哟哟哟哟哟哟| 99九九线精品视频在线观看视频| 中文在线观看免费www的网站| 两个人的视频大全免费| 国产私拍福利视频在线观看| 18禁黄网站禁片午夜丰满| 特级一级黄色大片| 五月玫瑰六月丁香| 女的被弄到高潮叫床怎么办 | 97超级碰碰碰精品色视频在线观看| av在线亚洲专区| 久久99热这里只有精品18| 午夜精品久久久久久毛片777| 欧美日韩国产亚洲二区| 国产成年人精品一区二区| 一级黄片播放器| 成人国产综合亚洲| 亚洲无线在线观看| 日本一本二区三区精品| 午夜福利欧美成人| bbb黄色大片| 人妻制服诱惑在线中文字幕| 精品久久久久久久久av| 亚洲欧美精品综合久久99| 动漫黄色视频在线观看| 男女之事视频高清在线观看| 99热精品在线国产| 国产主播在线观看一区二区| 99热这里只有是精品在线观看| 亚洲精品一卡2卡三卡4卡5卡| 精品久久久久久久久av| 久久精品国产鲁丝片午夜精品 | 老司机午夜福利在线观看视频| 九九在线视频观看精品| 国内毛片毛片毛片毛片毛片| 99热精品在线国产| 丰满的人妻完整版| 女人被狂操c到高潮| 日本成人三级电影网站| 听说在线观看完整版免费高清| 欧美高清性xxxxhd video| 亚洲成人久久性| 校园人妻丝袜中文字幕| 欧美潮喷喷水| 99热网站在线观看| 蜜桃亚洲精品一区二区三区| 大型黄色视频在线免费观看| 嫁个100分男人电影在线观看| 在线国产一区二区在线| 国产伦精品一区二区三区视频9| 国产淫片久久久久久久久| 亚洲最大成人中文| 露出奶头的视频| 蜜桃久久精品国产亚洲av| 久久久久久九九精品二区国产| 亚洲av第一区精品v没综合| 国产色爽女视频免费观看| 欧美又色又爽又黄视频| videossex国产| 精品久久久久久久人妻蜜臀av| 少妇丰满av| 久99久视频精品免费| 亚洲国产精品合色在线| 亚洲狠狠婷婷综合久久图片| 欧美一级a爱片免费观看看| 精品久久久久久,| 成年人黄色毛片网站| 国产伦一二天堂av在线观看| 国产精品,欧美在线| 久久久久九九精品影院| 级片在线观看| 看十八女毛片水多多多| 草草在线视频免费看| 欧美日韩乱码在线| 日本黄色片子视频| 国产毛片a区久久久久| 九色成人免费人妻av| 99热这里只有精品一区| 欧美日本亚洲视频在线播放| 成人av在线播放网站| 国产精品一区二区三区四区免费观看 | 日本精品一区二区三区蜜桃| 亚洲人成网站在线播| 此物有八面人人有两片| 欧美+日韩+精品| 日韩欧美一区二区三区在线观看| 在线看三级毛片| 丰满乱子伦码专区| 十八禁国产超污无遮挡网站| 男女下面进入的视频免费午夜| 久99久视频精品免费| 日本欧美国产在线视频| 网址你懂的国产日韩在线| 成人特级av手机在线观看| 亚洲天堂国产精品一区在线| avwww免费| 热99re8久久精品国产| av天堂中文字幕网| 岛国在线免费视频观看| 看黄色毛片网站| 在线免费观看的www视频| 国产单亲对白刺激| 亚洲真实伦在线观看| 久久国内精品自在自线图片| 欧美精品啪啪一区二区三区| 成人特级av手机在线观看| 国产又黄又爽又无遮挡在线| 在线观看一区二区三区| 色视频www国产| 悠悠久久av| 精品久久久久久久末码| 熟妇人妻久久中文字幕3abv| 国产一区二区在线av高清观看| 国产精品永久免费网站| 久久精品国产自在天天线| x7x7x7水蜜桃| 欧美一区二区精品小视频在线| 一个人观看的视频www高清免费观看| 看片在线看免费视频| 99热网站在线观看| 欧美国产日韩亚洲一区| 成人国产麻豆网| 最新在线观看一区二区三区| 精品人妻1区二区| 男人和女人高潮做爰伦理| 国产黄片美女视频| 亚洲黑人精品在线| 日本熟妇午夜| 久久久久久大精品| 国产欧美日韩精品一区二区| 变态另类丝袜制服| 久久久精品大字幕| 亚洲av免费在线观看| 国产精品久久电影中文字幕| 欧美日韩精品成人综合77777| x7x7x7水蜜桃| 久久欧美精品欧美久久欧美| 别揉我奶头~嗯~啊~动态视频| 亚洲人成网站在线播放欧美日韩| 性插视频无遮挡在线免费观看| 欧美人与善性xxx| 亚洲国产欧美人成| 亚洲三级黄色毛片| 亚洲国产欧洲综合997久久,| 亚洲在线自拍视频| 国产精品久久久久久精品电影| 国产v大片淫在线免费观看| 欧美又色又爽又黄视频| 很黄的视频免费| 国产女主播在线喷水免费视频网站 | 久久午夜亚洲精品久久| 亚洲第一区二区三区不卡| 精品久久久久久久人妻蜜臀av| 国产主播在线观看一区二区| 内地一区二区视频在线| 日本 av在线| 久久精品国产99精品国产亚洲性色| 男女啪啪激烈高潮av片| 国产美女午夜福利| 又粗又爽又猛毛片免费看| 精品欧美国产一区二区三| 亚洲人成网站在线播放欧美日韩| 18禁裸乳无遮挡免费网站照片| 国产白丝娇喘喷水9色精品| 欧美日韩乱码在线| 伊人久久精品亚洲午夜| 超碰av人人做人人爽久久| 国产av麻豆久久久久久久| 五月玫瑰六月丁香| 国产aⅴ精品一区二区三区波| 国产不卡一卡二| 蜜桃亚洲精品一区二区三区| 精品久久久久久久久亚洲 | 丰满的人妻完整版| 中文字幕人妻熟人妻熟丝袜美| 久久精品久久久久久噜噜老黄 | 免费不卡的大黄色大毛片视频在线观看 | АⅤ资源中文在线天堂| 国产探花在线观看一区二区| 精品久久久久久久久av| 亚洲av中文字字幕乱码综合| a级毛片a级免费在线| 国产激情偷乱视频一区二区| 日本精品一区二区三区蜜桃| 亚洲三级黄色毛片| 精品久久国产蜜桃| 午夜免费男女啪啪视频观看 | 一级黄片播放器| 国产爱豆传媒在线观看| 无人区码免费观看不卡| 精品人妻1区二区| 一级黄片播放器| 亚洲精品亚洲一区二区| 国产成人av教育| 成人av一区二区三区在线看| 日本欧美国产在线视频| 国产女主播在线喷水免费视频网站 | 久久精品国产亚洲网站| 人人妻,人人澡人人爽秒播| 三级毛片av免费| 午夜福利高清视频| 久久久国产成人免费| 不卡一级毛片| 99在线视频只有这里精品首页| 日本-黄色视频高清免费观看| 高清毛片免费观看视频网站| 精品一区二区三区视频在线观看免费| 国产精华一区二区三区| 悠悠久久av| 99在线视频只有这里精品首页| 欧美日韩黄片免| 国产亚洲av嫩草精品影院| 亚洲第一区二区三区不卡| 国产精华一区二区三区| 悠悠久久av| 国产真实乱freesex| 成人av一区二区三区在线看| 能在线免费观看的黄片| av天堂在线播放| 久久精品国产亚洲网站| 亚洲国产精品sss在线观看| 亚洲精品456在线播放app | 国产色爽女视频免费观看| 国产精品国产三级国产av玫瑰| 又爽又黄a免费视频| 亚洲av.av天堂| 久久久久精品国产欧美久久久| 丰满的人妻完整版| 欧美一级a爱片免费观看看| 国产欧美日韩精品亚洲av| 国产精品电影一区二区三区| 男女边吃奶边做爰视频| 韩国av一区二区三区四区| 舔av片在线| 亚洲无线观看免费| 久久久久久久精品吃奶| 国产精品亚洲一级av第二区| 国产亚洲精品久久久com| 999久久久精品免费观看国产| 国产亚洲av嫩草精品影院| 身体一侧抽搐| 无人区码免费观看不卡| 亚洲精品一区av在线观看| 成人鲁丝片一二三区免费| 国产 一区精品| 欧美极品一区二区三区四区| 最新中文字幕久久久久| 国产aⅴ精品一区二区三区波| 88av欧美| 小蜜桃在线观看免费完整版高清| 亚洲av一区综合| 不卡视频在线观看欧美| 久久午夜福利片| 国产亚洲精品久久久久久毛片| 高清毛片免费观看视频网站| 国产高清视频在线播放一区| 亚洲成a人片在线一区二区| 精品人妻视频免费看| 国产v大片淫在线免费观看| 午夜a级毛片| 级片在线观看| 床上黄色一级片| 中文字幕高清在线视频| 日本爱情动作片www.在线观看 | 久久婷婷人人爽人人干人人爱| 色尼玛亚洲综合影院| 日韩欧美三级三区| 波多野结衣巨乳人妻| av在线天堂中文字幕| 伦精品一区二区三区| 欧美zozozo另类| 国产老妇女一区| 精品欧美国产一区二区三| 久久精品人妻少妇| 极品教师在线视频| 99久久中文字幕三级久久日本| 久久久精品大字幕| 亚洲性夜色夜夜综合| 少妇熟女aⅴ在线视频| 毛片一级片免费看久久久久 | 欧美激情在线99| 三级男女做爰猛烈吃奶摸视频| 免费观看的影片在线观看| 国产精品综合久久久久久久免费| x7x7x7水蜜桃| 午夜a级毛片| 国产 一区 欧美 日韩| 国产精品一区二区三区四区久久| 中文字幕av在线有码专区| 制服丝袜大香蕉在线| 欧美区成人在线视频| 国产亚洲精品久久久com| 22中文网久久字幕| 夜夜夜夜夜久久久久| 欧美日韩精品成人综合77777| 久久精品国产鲁丝片午夜精品 | 国产成人一区二区在线| 亚洲精品影视一区二区三区av| 亚洲第一区二区三区不卡| 国产麻豆成人av免费视频| 村上凉子中文字幕在线| 亚洲美女搞黄在线观看 | 日韩亚洲欧美综合| 自拍偷自拍亚洲精品老妇| 中文字幕免费在线视频6| 国产精品国产三级国产av玫瑰| 日日摸夜夜添夜夜添小说| 俄罗斯特黄特色一大片| av天堂中文字幕网| 熟女人妻精品中文字幕| www日本黄色视频网| 在线播放国产精品三级| 九色成人免费人妻av| 91狼人影院| 少妇猛男粗大的猛烈进出视频 | 小说图片视频综合网站| 可以在线观看的亚洲视频| 久久人人精品亚洲av| 精品99又大又爽又粗少妇毛片 | 国产午夜精品论理片| 欧美成人a在线观看| 三级国产精品欧美在线观看| 两人在一起打扑克的视频| 国内揄拍国产精品人妻在线| 联通29元200g的流量卡| 免费大片18禁| 女同久久另类99精品国产91| 午夜爱爱视频在线播放| 亚洲精品日韩av片在线观看| 久久精品国产亚洲av天美| 美女高潮的动态| 99久久久亚洲精品蜜臀av| 真人一进一出gif抽搐免费| 国产男靠女视频免费网站| 久久精品国产99精品国产亚洲性色| 2021天堂中文幕一二区在线观| 国产中年淑女户外野战色| avwww免费| 日韩大尺度精品在线看网址| 免费观看人在逋| 小说图片视频综合网站| 国产中年淑女户外野战色| 午夜a级毛片| 日韩大尺度精品在线看网址| 搡老妇女老女人老熟妇| 一区二区三区激情视频| 国产日本99.免费观看| 国产精品av视频在线免费观看| 精品国内亚洲2022精品成人| 日本爱情动作片www.在线观看 | 亚洲av成人av| 国产精品99久久久久久久久| 春色校园在线视频观看| 国产男靠女视频免费网站| 久久精品国产鲁丝片午夜精品 | 国产毛片a区久久久久| 亚洲最大成人手机在线| 亚洲欧美日韩无卡精品| 精品久久久久久久久久免费视频| 我要搜黄色片| 国产熟女欧美一区二区| 麻豆av噜噜一区二区三区| 国产免费av片在线观看野外av| 真人做人爱边吃奶动态| 欧美zozozo另类| 12—13女人毛片做爰片一| 亚洲精品456在线播放app | 国产伦一二天堂av在线观看| 亚洲欧美清纯卡通| 欧美日本亚洲视频在线播放| 观看免费一级毛片| 乱人视频在线观看| 在线天堂最新版资源| 亚洲成a人片在线一区二区| 日日撸夜夜添| 色5月婷婷丁香| 男人的好看免费观看在线视频| 99国产极品粉嫩在线观看| 精品久久久久久久久久久久久| 亚洲精品国产成人久久av| 一个人看的www免费观看视频| 18禁黄网站禁片午夜丰满| 国产精品av视频在线免费观看| 看片在线看免费视频| 欧美黑人巨大hd| 久久热精品热| 成人二区视频| 国产亚洲91精品色在线| 国内精品一区二区在线观看| 成人性生交大片免费视频hd| 成人国产麻豆网| 午夜视频国产福利| 性插视频无遮挡在线免费观看| 91午夜精品亚洲一区二区三区 | 久久久久久久午夜电影| 亚洲专区国产一区二区| x7x7x7水蜜桃| 丰满人妻一区二区三区视频av| 国产伦一二天堂av在线观看| 久久午夜亚洲精品久久| 日韩中字成人| 日韩一区二区视频免费看| 国产精品久久久久久av不卡| 日韩欧美国产一区二区入口| 国产69精品久久久久777片| 午夜福利欧美成人| 日本撒尿小便嘘嘘汇集6| 观看免费一级毛片| 国产精品av视频在线免费观看| 黄色丝袜av网址大全| 久久亚洲精品不卡| 欧美高清性xxxxhd video| avwww免费| 老女人水多毛片| 国产欧美日韩精品亚洲av| av专区在线播放| 99九九线精品视频在线观看视频| 日韩,欧美,国产一区二区三区 | 99在线视频只有这里精品首页| 亚洲 国产 在线| 国产单亲对白刺激| 草草在线视频免费看| 亚洲av成人精品一区久久| 欧美日韩中文字幕国产精品一区二区三区| 国产午夜精品论理片| 一个人观看的视频www高清免费观看| 国产熟女欧美一区二区| 简卡轻食公司| 亚洲美女黄片视频| 亚洲美女搞黄在线观看 | 成人无遮挡网站| 成人国产一区最新在线观看| 中文在线观看免费www的网站| 国产伦人伦偷精品视频| 两个人的视频大全免费| 色尼玛亚洲综合影院| 超碰av人人做人人爽久久| 成人av一区二区三区在线看| 在线播放国产精品三级| 欧美日韩乱码在线| 亚洲专区中文字幕在线| 国产 一区 欧美 日韩| 欧美成人性av电影在线观看| 大型黄色视频在线免费观看| 九色国产91popny在线| 琪琪午夜伦伦电影理论片6080| 韩国av一区二区三区四区| 制服丝袜大香蕉在线| 国产亚洲精品久久久com| 两个人的视频大全免费| 欧美色欧美亚洲另类二区| 97人妻精品一区二区三区麻豆| 亚洲av不卡在线观看| 久久久久久伊人网av| 国产色婷婷99| 99九九线精品视频在线观看视频| 亚洲av五月六月丁香网| 好男人在线观看高清免费视频| 久久这里只有精品中国| 日日摸夜夜添夜夜添av毛片 | 国产极品精品免费视频能看的| 午夜福利成人在线免费观看|