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

    遼寧省2001-2020年水源涵養(yǎng)量變化及其影響因素

    2022-04-02 06:14:18呂樂婷
    水土保持通報 2022年1期
    關鍵詞:產(chǎn)水量分異水源

    呂樂婷, 李 青, 楊 勇

    (遼寧師范大學 地理科學學院, 遼寧 大連 116029)

    生態(tài)系統(tǒng)服務是人類賴以生存和發(fā)展的基礎,水源涵養(yǎng)作為其提供的眾多服務之一,是水資源可持續(xù)利用的保障[1]。生態(tài)環(huán)境惡化及人類不合理開發(fā)利用,導致水資源短缺,已成為世界各國共同關注的問題。水源涵養(yǎng)是一個不斷發(fā)展的動態(tài)概念,它是指在一定時空范圍和條件下,將水分保持在系統(tǒng)內(nèi)的過程和能力,其功能表現(xiàn)為涵養(yǎng)水源、攔蓄降水、調(diào)節(jié)河川流量、影響降雨量、凈化水質(zhì)等[2],其內(nèi)涵隨著人們對生態(tài)系統(tǒng)的進一步研究而逐步豐富。國內(nèi)研究始于20世紀80年代,集中探討不同森林類型涵養(yǎng)水源的作用[3]。隨著社會發(fā)展及研究的深入,其研究范圍不斷擴大,不同土地利用下的水源涵養(yǎng)成為研究熱點[4]。水源涵養(yǎng)量計算方法主要有降水儲存法、綜合蓄水能力法、地下徑流增長法等,常用模型主要有基于水量平衡的InVEST模型、SWAT模型、MIKE SHE模型等。

    InVEST(integrated valuation of ecosystem services and tradeoffs)模型,是由斯坦福大學、大自然保護協(xié)會、世界自然基金會聯(lián)合開發(fā)的生態(tài)系統(tǒng)服務功能評估模型。該模型突破了傳統(tǒng)評估方法的局限性,可動態(tài)分析、定量評估生態(tài)系統(tǒng)服務功能并以地圖的形式表現(xiàn)出來,其陸地生態(tài)系統(tǒng)模塊包括產(chǎn)水量、生境質(zhì)量、碳儲存及生物多樣性等模塊[5]。InVEST模型在國外應用相對成熟,Babbar等[6]利用馬爾科夫鏈和InVEST模型對印度薩里卡老虎保護區(qū)固碳量進行評估,結(jié)果表明研究區(qū)2000—2018年由于人類過度砍伐森林導致固碳量逐年減少,經(jīng)預測2035年實施森林管理后,碳排放量減少。Ashok等[7]利用InVEST模型量化印度海得拉巴市綠地的防洪服務,結(jié)果表明城市綠地和開放空間保留了44%~50%的降雨,降雨強度稍微增加,就會產(chǎn)生巨大的徑流量,造成相應的經(jīng)濟損失。國內(nèi)研究集中在2010年以后,研究方向主要集中于生態(tài)系統(tǒng)服務功能的評估。吳哲等[8]、包玉斌等[9]采用InVEST模型等方法分別對海南島的產(chǎn)水量、陜北黃土高原的水源涵養(yǎng)進行研究模擬。結(jié)果表明土地利用類型中的林地、草地和氣候因子中的降水是影響產(chǎn)水量、水源涵養(yǎng)的主要因素。此外,模型在其他地理環(huán)境差異性較大的區(qū)域均取得了較好的效果[10-11]。

    遼寧省作為東北糧食基地和老工業(yè)基地之一,生產(chǎn)、生活用水需求量較大,城市化及產(chǎn)業(yè)集聚與水污染加重,生態(tài)系統(tǒng)水源涵養(yǎng)功能不斷退化,水資源嚴重匱乏。近年來,遼寧省水生態(tài)系統(tǒng)研究取得了新進展,薛超[12]、李紅振[13]對遼寧省水資源承載力、森林水源涵養(yǎng)能力進行研究。本研究基于遼寧省2001—2020年土地利用狀況,結(jié)合氣象,土壤,地形等數(shù)據(jù),使用InVEST模型評估遼寧省生態(tài)系統(tǒng)水源涵養(yǎng)量時空變化特征,并探索其與氣候因素、土地利用、生態(tài)、社會經(jīng)濟方面的相應關系,以期為遼寧省的生態(tài)環(huán)境保護、生態(tài)系統(tǒng)水源涵養(yǎng)功能恢復提供決策依據(jù)。

    1 研究區(qū)概況與數(shù)據(jù)來源

    1.1 研究區(qū)概況

    遼寧省位于118°53′—125°46′E,38°43′—43°26′N,地跨遼河、渾河、大凌河、太子河、繞陽河、鴨綠江六大水系,與吉林、內(nèi)蒙古及河北等地區(qū)接壤。屬于大陸性季風氣候,四季分明,夏季炎熱多雨,冬季寒冷漫長,2001—2020年多年均溫6—11 ℃,平均降水量為429.25~999.61 mm。氣溫由東部向西部遞增;降水量由沿海向內(nèi)陸遞減,東部的撫順、本溪和丹東市降水量較多,西部的朝陽、葫蘆島等市降水量較少。

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

    InVEST模型的產(chǎn)水量模塊(即水源供給模塊)所需參數(shù)有年降水量、年潛在蒸散發(fā)量、土壤最大根系埋藏深度、植物可利用水量和土地利用/覆被和生物物理參數(shù)表等。

    (1) 土地利用/覆被。本研究采用2005,2010,2015,2020年共4期土地利用/覆被數(shù)據(jù),分辨率為1 km×1 km,均來源于中國科學院資源環(huán)境科學數(shù)據(jù)中心(https:∥www.resdc.cn),數(shù)據(jù)生產(chǎn)制作是以各期Landsat TM/ETM遙感影像為數(shù)據(jù)源,通過人機交互目視解譯生成。根據(jù)研究需要將原土地利用/覆被類型數(shù)據(jù)歸類并調(diào)整為耕地、林地、草地、水域、建設用地、未利用地6類。遼寧省土地利用變化詳情見表1。耕地和林地為主要用地類型,占全省總面積的83%。2005—2020年,遼寧省各土地利用類型變化較大,其中耕地和草地面積分別減少4 097,4 653 km2,降幅分別為3.16%,3.22%;林地、建設用地增加,面積增長4 214,4 129 km2,增幅分別為2.53%,2.75%;水域和未利用地變化小,面積分別增加1 310,358 km2,增幅分別為0.86%,0.23%。

    表1 遼寧省土地利用類型結(jié)構

    (2) 氣象數(shù)據(jù)。本研究所需氣象數(shù)據(jù)均來源于中國氣象局數(shù)據(jù)中心(http:∥data.cma.cn)。2001—2020年平均降水量分布特征(圖1a)利用遼寧省27個氣象站通過克里金插值法對所有站點數(shù)據(jù)進行空間插值得到;年潛在蒸散量利用氣溫、風速、日照時數(shù)、氣壓日值數(shù)據(jù)基于聯(lián)合國糧農(nóng)組織FAO修正的Penman-Monteith公式計算[14],再進行空間插值(圖1b)。

    圖1 2001-2020年平均遼寧省降水量、潛在蒸散量空間分布特征

    降水量和潛在蒸散量的空間分布特征呈顯著差異,降水主要集中于遼寧省東部的丹東、撫順、本溪市,降水量在726.13~999.61 mm之間,西部的朝陽市、阜新市和錦州市是降水低值區(qū),降水范圍為429.25~550.67 mm。西部的朝陽市、阜新市和錦州市的潛在蒸散較強,蒸散量在977.71~1 061.91 mm之間,東部的丹東、本溪、撫順、鞍山市潛在蒸散較弱,蒸散范圍為773.54~887.04 mm。

    (3) 土壤數(shù)據(jù)。本研究所需土壤數(shù)據(jù)均來源于聯(lián)合國糧農(nóng)組織和維也納國際應用系統(tǒng)研究所構建的世界土壤數(shù)據(jù)庫,中國境內(nèi)數(shù)據(jù)源為第二次全國土地調(diào)查(2009年)南京土壤所提供的1∶100萬土壤數(shù)據(jù)。土壤厚度又稱根系最大埋藏深度,是指受物理或化學特性的影響而強烈阻礙根系穿透時的土層深度。土壤有效含水量計算是基于土壤屬性數(shù)據(jù),由美國華盛頓州立大學開發(fā)的土壤水特性軟件中的soil-water-characterstics(SWCT)模塊計算獲得。計算方法如下:根據(jù)黏土(caly)、砂(sand)、有機物(organic matter)、鹽度(salinity)、砂礫(gravel)等參數(shù)計算得出的凋萎系數(shù)和田間持水量,田間持水量和凋萎系數(shù)差值即土壤有效含水量。

    (4) 生物物理參數(shù)表。包括土地利用/覆被編碼、作物根系深度、蒸散系數(shù),其中作物根系深度來自其他學者對全球范圍內(nèi)植被最大根深的研究獲得[15],蒸散系數(shù)是根據(jù)聯(lián)合國糧農(nóng)組織作物蒸散系數(shù)指南、InVEST用戶指南獲得。

    2 研究方法

    2.1 水源涵養(yǎng)量的計算

    InVEST模型的產(chǎn)水模塊,是基于水量平衡原理,對產(chǎn)水量進行估算,其實質(zhì)為每一個柵格降水量和實際蒸散量的差值,計算方式如下:

    式中:Yxj為柵格單元x中土地覆被/利用類型j的年產(chǎn)水量,主要包括地表徑流量、土壤含水量、枯落物持水量、冠層截留量;Px為柵格單元x的年降水量; AETxj為柵格單元x中的土地覆被/利用類型j的年實際蒸散量,由(2) 式計算:

    式中:wx為自然氣候條件下表示土壤性質(zhì)的一個無量綱非物理參數(shù),由公式(3)計算:

    式中:Z為經(jīng)驗常數(shù),也叫季節(jié)性因子,代表區(qū)域降水分布及其他水文地質(zhì)特征。根據(jù)相關研究[16],對于冬季降水為主的地區(qū),Z值接近10,而對于降水均勻分布的濕潤地區(qū)和夏季降水為主的地區(qū),Z值接近1。AWCx為柵格單元x的土壤有效含水量,由土層深度和理化性質(zhì)決定。Rxj為Bydyko干燥指數(shù),由公式(4)計算:

    式中:ET0為柵格單元x中土地覆被類型j的植被蒸散系數(shù);Kxj為參考作物蒸散; ET0由公式(5)計算:

    ET0=

    式中:ET0為潛在蒸散量(mm/d);Rn為地表凈輻射〔MJ/(m·d)〕;G為土壤熱通量〔MJ/(m2·d)〕;γ為干濕表常數(shù)(kPa/℃);Tmean為日平均溫度(℃);u2為2 m高處風速(m/s);es為飽和水氣壓(kPa);ea為實際水氣壓(kPa)。

    水源涵養(yǎng)量的計算是InVEST模型在模擬出研究區(qū)年產(chǎn)水量的基礎上考慮地形因素、地表徑流(流速系數(shù)大小)、土壤性質(zhì)(土壤飽和導水率)等因素對其的影響,計算徑流在柵格上停留時間得出,具體計算如公式(6):

    式中:retention為單位面積水源涵養(yǎng)量(mm);V是流速系數(shù),采用模型參數(shù)表數(shù)據(jù);Ksat為土壤飽和導水率(mm/h); TI為地形指數(shù),由公式(7)計算:

    式中:Darea為集水區(qū)柵格數(shù)量; Soildep為土層深度(mm);Pslope為坡度比。

    2.2 基于地理探測器影響因素的分析

    地理探測器是探測地理現(xiàn)象空間分異以及其背后驅(qū)動因素的一種統(tǒng)計學方法,主要包括風險探測、因子探測、生態(tài)探測和交互探測4個模塊,其具有同時探測多個因子在不同空間單元的不同影響作用及其交互關系的優(yōu)勢[17]。本文主要通過因子探測和交互探測兩個模塊對遼寧省水源涵養(yǎng)空間分異影響因素進行定量分析,“因子探測”表明單個因子在一定程度上解釋生態(tài)系統(tǒng)服務功能的空間分異,“交互探測”識別不同自變量因子之間的交互作用。計算如公式(8):

    3 結(jié)果與分析

    3.1 模型校驗

    基于InVEST模型計算研究區(qū)產(chǎn)水量時,依據(jù)遼寧省水資源公報數(shù)據(jù)對模型參數(shù)進行率定,通過調(diào)整Z值(Zhang系數(shù))使模擬值更加接近實測值,模擬結(jié)果如圖2。2001—2020年,遼寧省多年平均徑流量為2.53×1010m3,與遼寧省水資源公報公布的多年平均徑流量2.61×1010m3較接近,平均誤差率為3.92%,說明模擬結(jié)果較好。

    圖2 2001-2020年遼寧省InVEST模型結(jié)果驗證

    3.2 遼寧省水源涵養(yǎng)量的年際變化特征

    2001—2020年遼寧省年均降水量、潛在蒸散量、產(chǎn)水量、水源涵養(yǎng)量變化見圖3。遼寧省多年平均降水量為644.35 mm,在研究時段內(nèi)呈顯著增加趨勢(p<0.05),增長速率為14.6 mm/(10 a);多年平均潛在蒸散量為931.88 mm,增長速率為49.92 mm/(10 a),未達顯著性水平;多年平均產(chǎn)水量和水源涵養(yǎng)量分別為177.22 mm和95.28 mm,最低值分別出現(xiàn)在2014年和2009年,分別為76.47 mm和40.52 mm,最高值均出現(xiàn)在2010年,分別為412.73 mm和186.64 mm。多年來,遼寧省產(chǎn)水量和水源涵養(yǎng)量均呈現(xiàn)不顯著增長趨勢,增長速率分別為19.36 mm/(10 a)和13.35 mm/(10 a)。

    圖3 2001-2020年遼寧省降水量、潛在蒸散量、產(chǎn)水量及水源涵養(yǎng)量的年際變化

    3.3 遼寧省水源涵養(yǎng)量的空間分布

    遼寧省多年平均產(chǎn)水量及水源涵養(yǎng)量空間變化見圖4。水源涵養(yǎng)總體呈從東—西遞減的空間分布特征,高值區(qū)包括撫順市、本溪市和丹東市,水源涵養(yǎng)量在180.73~245.24 mm之間。這是由于該區(qū)域降水充沛,林地分布廣,蒸散發(fā)為全省最低值,故產(chǎn)水量大,水源涵養(yǎng)能力強。中、西部的盤錦市、錦州市、阜新市、朝陽市、沈陽市和葫蘆島市是低值區(qū),水源涵養(yǎng)量在24.25~46.23 mm之間,該區(qū)域地形以平原、低山和丘陵為主,降水偏少,蒸散強,受地形與蒸散影響,產(chǎn)水量與水源涵養(yǎng)量均較少。處于東部向中西部過渡的大連市、營口市、鐵嶺市、遼陽市和鞍山市水源涵養(yǎng)量在71.45~111.04 mm之間。

    圖4 遼寧省產(chǎn)水量和水源涵養(yǎng)量空間分布特征

    3.4 水源涵養(yǎng)與不同因子的相關分析結(jié)果

    參考李怡穎等[18],王亞慧等[19]的研究及數(shù)據(jù)可獲取性,本研究分別從氣候、土地利用、生態(tài)、社會經(jīng)濟方面選擇降水量、蒸散量、氣溫、地塊破碎度、林地比例、地表水資源量、山區(qū)入滲量、城鎮(zhèn)化率等參數(shù),探索其與水源涵養(yǎng)量的關系。經(jīng)皮爾遜相關分析可得,水源涵養(yǎng)量與降水量、林地比例、地表水資源量和山區(qū)入滲量呈顯著正相關,與蒸散量、氣溫、建設用地比例、地塊破碎度呈顯著負相關(表2)。

    表2 遼寧省不同影響因素與水源涵養(yǎng)的相關關系

    水源涵養(yǎng)量與降水量、地表水資源量、林地比例和山區(qū)入滲量相關性較高,相關系數(shù)分別為0.95,0.67,0.61,0.56。遼寧省東部山區(qū)林地面積大,人類活動強度小,地塊較完整,植被覆蓋度高,降水、地表水資源量豐富,水源涵養(yǎng)量明顯高于其他地區(qū),水源涵養(yǎng)能力強。蒸散量、氣溫、地塊破碎度、建設用地比例與水源涵養(yǎng)量的相關系數(shù)分別為-0.93,-0.64,-0.72,-0.55,相關性較強,遼寧省中、西部地區(qū)地形以平原、低山和丘陵為主,氣溫高于東部山地,蒸發(fā)旺盛,建設用地擴張,地塊破碎度增加,水源涵養(yǎng)能力弱。此外,城鎮(zhèn)化率增長在一定程度上對水源涵養(yǎng)起促進作用,農(nóng)作物播種面積增加,人類過度使用耕地,灌溉率高,影響土壤持水功能,導致水源涵養(yǎng)功能退化。

    3.5 遼寧省水源涵養(yǎng)的空間分異驅(qū)動力分析

    遼寧省各單因子對水源涵養(yǎng)空間異質(zhì)性的解釋力見表3。

    表3 遼寧省單因子對水源涵養(yǎng)空間異質(zhì)性解釋力(q)

    降水、蒸散、建設用地比例、地表水資源量、氣溫和林地比例是影響遼寧省水源涵養(yǎng)空間分異的主導因素,解釋力分別為0.9,0.85,0.83,0.77,0.62,0.61。地塊破碎度、農(nóng)業(yè)灌溉率、山區(qū)入滲量、生態(tài)環(huán)境用水的解釋力在0.4~0.59之間,分別為0.59,0.49,0.44,0.43。其他因子如農(nóng)作物面積、城鎮(zhèn)化率的解釋力有限,分別為0.36,0.33。總體上,氣候因子對遼寧省水源涵養(yǎng)空間分異影響遠大于其他因子;土地利用方面,建設用地比例解釋力大于林地比例與地塊破碎度;生態(tài)因子與社會經(jīng)濟因子總體解釋力相當,但生態(tài)因子中的地表水資源量對水源涵養(yǎng)空間分異的影響更大,解釋力更強。

    遼寧省多因子的交互作用對水源涵養(yǎng)空間異質(zhì)性解釋力均高于單個因子的解釋力(表4)。氣候因子中的降水和蒸散與其他各因子的交互作用解釋力最強,均在0.9以上,氣溫與其他因子交互作用稍弱,與農(nóng)業(yè)灌溉率交互作用的解釋力最低,為0.72;土地利用中各因子的交互作用解釋力表現(xiàn)為:建設用地比例>地塊破碎度>林地比例;生態(tài)因子和社會經(jīng)濟因子交互作用解釋力相當。綜上可知,無論單因子作用還是多因子交互作用,氣候因子都是影響遼寧省水源涵養(yǎng)空間分異的主導因子。

    表4 遼寧省不同因子對水源涵養(yǎng)空間異質(zhì)性交互作用解釋力(q)

    4 討論與結(jié)論

    4.1 結(jié) 論

    (1) 2001—2020年遼寧省多年平均產(chǎn)水量和水源涵養(yǎng)量分別為177.22 mm和95.28 mm,增長速率分別為19.36 mm/(10 a)和13.35 mm/(10 a)。水源涵養(yǎng)量空間變化總體上從東—西遞減,高值區(qū)分布于丹東市、撫順市和本溪市,低值區(qū)位于盤錦市、錦州市、阜新市和朝陽市。

    (2) 水源涵養(yǎng)量與降水量、地表水資源量、林地比例和山區(qū)入滲量呈顯著正相關,與蒸散量、氣溫、地塊破碎度、建設用地比例呈顯著負相關。

    (3) 氣候因子與土地利用因子是遼寧省水源涵養(yǎng)空間分異的主導因子。降水、蒸散、建設用地比例、地表水資源量、氣溫和林地比例的解釋力0.6以上,地塊破碎度、農(nóng)業(yè)灌溉率、山區(qū)入滲量、生態(tài)環(huán)境用水的解釋力在0.4~0.59之間,農(nóng)作物面積、城鎮(zhèn)化率的解釋力有限。多因子交互解釋力均高于單個因子,其中降水、蒸散與各因子共同作用于水源涵養(yǎng)時部分解釋力可達0.92以上。

    4.2 討 論

    本研究結(jié)果表明遼寧省2001—2020年水源涵養(yǎng)總量呈波動上升變化。2001年水源涵養(yǎng)總量11.95×109m3,2020年上升至21.15×109m3,總體上呈上升趨勢。遼寧省水源涵養(yǎng)量多年均值自西部丘陵區(qū)向東部山區(qū)增加,造成這種分異規(guī)律的原因是東、西部地區(qū)在氣候、用地類型等方面的差異。東部地區(qū)地形以山地為主,林地約占全省林地面積50%,森林蓋度高且多為天然林,降水豐富,蒸散量低于其他地區(qū),產(chǎn)水能力強,水源涵養(yǎng)量大。中、西部地區(qū)與內(nèi)蒙古河北接壤,平原、低山和丘陵比例大,森林資源相對匱乏,植被覆蓋度低,氣候干燥,日照時數(shù)長,降水少且蒸散強,產(chǎn)水量低、水源涵養(yǎng)能力弱。

    氣候、土地利用類型與水源涵養(yǎng)相關性較強,是影響水源涵養(yǎng)變化的主要因素,其中降水與蒸散對水源涵養(yǎng)有直接影響。在6類土地利用類型中,林地具有根系深度大的特征,截留、下滲能力強,因而水源涵養(yǎng)量大,與水源涵養(yǎng)呈顯著正相關性,在城市化過程中應注重林地的開發(fā)與保護。無論單因子作用還是多因子交互作用,氣候因子都是遼寧省水源涵養(yǎng)空間分異的主要控制因素,降水量、蒸散量差異較大的地區(qū),水源涵養(yǎng)量不同,表現(xiàn)出明顯的空間差異性。此外,用地類型是水源涵養(yǎng)空間分異的驅(qū)動力之一,受人類活動影響較大,用地類型轉(zhuǎn)變影響水源涵養(yǎng)功能,因此優(yōu)化土地利用結(jié)構,對水源涵養(yǎng)功能恢復有重要作用。

    水源涵養(yǎng)量的計算是一個復雜的綜合性過程,參考可獲取的水資源公報及統(tǒng)計年鑒的相關數(shù)據(jù),對模型所需參數(shù)進行校驗,在今后的研究中應加強野外數(shù)據(jù)監(jiān)測,增加研究結(jié)果的可信度。相較于其他研究,本論文集中研究行政區(qū)尺度的水源涵養(yǎng)時空動態(tài)及歸因分析,在不同生態(tài)系統(tǒng)類型,例如森林生態(tài)系統(tǒng)方面研究較少,還需進一步完善。此外,Z值是表征多年平均降水特征的常數(shù),不同地區(qū)的年際降水特征存在差異,長時間序列情況下,Z值的確定需不斷調(diào)整,耗時較長。通過Z值調(diào)整提高模擬精度,使模擬結(jié)果更接近實測值。

    猜你喜歡
    產(chǎn)水量分異水源
    保護水源
    品牌研究(2022年18期)2022-06-29 05:33:24
    基于INVEST 模型的資水流域產(chǎn)水量及其對環(huán)境響應的評估分析
    贛江流域產(chǎn)水功能對土地利用變化的響應
    雅魯藏布江下游產(chǎn)水量時空演變及對氣候和土地利用變化的響應
    南水源頭清如許
    當代陜西(2019年24期)2020-01-18 09:14:18
    重慶市臭氧時空分異及其影響因素研究
    尋找水源
    平泉縣下營坊雜巖體分異演化及其成巖成礦
    煤層氣井長沖程、大泵徑排采設備的研究及應用
    北京市1989-2010年地表溫度時空分異特征分析
    岛国在线免费视频观看| 久久久久免费精品人妻一区二区| 人体艺术视频欧美日本| 亚洲av成人精品一二三区| 深爱激情五月婷婷| 乱系列少妇在线播放| 国产 一区 欧美 日韩| 人妻制服诱惑在线中文字幕| 我的老师免费观看完整版| 午夜a级毛片| 久久久久久久午夜电影| 国产亚洲精品久久久com| 直男gayav资源| 国产精品一区二区三区四区免费观看| 国产 一区精品| 舔av片在线| 三级毛片av免费| 青春草视频在线免费观看| 熟女电影av网| 免费不卡的大黄色大毛片视频在线观看 | 亚洲av成人av| 午夜福利视频1000在线观看| 亚洲av成人精品一区久久| 美女脱内裤让男人舔精品视频| 26uuu在线亚洲综合色| 欧美一区二区亚洲| 亚洲精品456在线播放app| 中文资源天堂在线| 成人三级黄色视频| 国产免费视频播放在线视频 | 两性午夜刺激爽爽歪歪视频在线观看| 欧美日本亚洲视频在线播放| 校园人妻丝袜中文字幕| 中文字幕制服av| 欧美+日韩+精品| 日本与韩国留学比较| 噜噜噜噜噜久久久久久91| 亚洲人与动物交配视频| 国产伦理片在线播放av一区| 国产 一区精品| 一个人免费在线观看电影| 青春草国产在线视频| 午夜视频国产福利| 汤姆久久久久久久影院中文字幕 | av视频在线观看入口| 婷婷色av中文字幕| 白带黄色成豆腐渣| 在线观看一区二区三区| 两个人视频免费观看高清| 如何舔出高潮| 最近中文字幕2019免费版| 亚洲怡红院男人天堂| 卡戴珊不雅视频在线播放| 国产一级毛片在线| 在线a可以看的网站| 久久婷婷人人爽人人干人人爱| 卡戴珊不雅视频在线播放| 精品久久久久久成人av| 国产老妇女一区| 免费av毛片视频| 最后的刺客免费高清国语| 在线播放国产精品三级| 你懂的网址亚洲精品在线观看 | 国产精品人妻久久久久久| 日韩欧美 国产精品| АⅤ资源中文在线天堂| 久久久久网色| 国产精品不卡视频一区二区| 精华霜和精华液先用哪个| 久久99热这里只频精品6学生 | 亚洲精品乱码久久久久久按摩| 欧美激情久久久久久爽电影| 亚洲av成人精品一二三区| 18禁在线无遮挡免费观看视频| 久久久久久久久久久丰满| 欧美最新免费一区二区三区| 七月丁香在线播放| 久99久视频精品免费| 丝袜喷水一区| ponron亚洲| 久久精品综合一区二区三区| a级毛色黄片| 麻豆成人av视频| 97人妻精品一区二区三区麻豆| 精品免费久久久久久久清纯| 亚洲天堂国产精品一区在线| 老司机影院毛片| 在线免费观看的www视频| 一区二区三区乱码不卡18| 国产综合懂色| 男女啪啪激烈高潮av片| 偷拍熟女少妇极品色| 国产单亲对白刺激| 美女脱内裤让男人舔精品视频| 国产欧美另类精品又又久久亚洲欧美| 国产精品一区www在线观看| 日韩一区二区视频免费看| 日韩精品青青久久久久久| 中文字幕熟女人妻在线| 成年av动漫网址| 色网站视频免费| 最近最新中文字幕免费大全7| 一区二区三区乱码不卡18| 国产免费又黄又爽又色| 久久99蜜桃精品久久| 校园人妻丝袜中文字幕| 男女那种视频在线观看| 2022亚洲国产成人精品| 国产成人a∨麻豆精品| 99久久人妻综合| 全区人妻精品视频| 久久久精品欧美日韩精品| 美女高潮的动态| 成人高潮视频无遮挡免费网站| 最新中文字幕久久久久| 国产老妇伦熟女老妇高清| 麻豆一二三区av精品| 国产精品国产三级专区第一集| 一级毛片久久久久久久久女| 高清毛片免费看| 熟女电影av网| 97在线视频观看| 建设人人有责人人尽责人人享有的 | 人人妻人人澡欧美一区二区| 精品久久久久久久末码| 成人午夜精彩视频在线观看| 国内揄拍国产精品人妻在线| 国产人妻一区二区三区在| 国产精品不卡视频一区二区| 男女下面进入的视频免费午夜| 亚洲国产最新在线播放| 赤兔流量卡办理| 久久99蜜桃精品久久| 欧美激情在线99| 99久久成人亚洲精品观看| 身体一侧抽搐| 91在线精品国自产拍蜜月| 亚洲国产精品久久男人天堂| 欧美三级亚洲精品| 国产高清国产精品国产三级 | 久久精品久久久久久噜噜老黄 | 日韩av不卡免费在线播放| 久久综合国产亚洲精品| 久久人妻av系列| 午夜福利高清视频| 久久6这里有精品| 老师上课跳d突然被开到最大视频| 日韩欧美 国产精品| 成人漫画全彩无遮挡| 亚洲av电影不卡..在线观看| 波野结衣二区三区在线| 国产精品一区二区三区四区久久| 男人狂女人下面高潮的视频| 亚洲国产高清在线一区二区三| 久久久久久国产a免费观看| 日本黄色片子视频| 国产成人精品久久久久久| 色综合站精品国产| 久久精品国产鲁丝片午夜精品| 日本黄大片高清| 成人漫画全彩无遮挡| 高清毛片免费看| 午夜福利网站1000一区二区三区| 99国产精品一区二区蜜桃av| 久热久热在线精品观看| 国产伦一二天堂av在线观看| 日本-黄色视频高清免费观看| 亚洲精品日韩av片在线观看| 久久久国产成人免费| 欧美日韩国产亚洲二区| 中文精品一卡2卡3卡4更新| 麻豆乱淫一区二区| 你懂的网址亚洲精品在线观看 | 亚洲经典国产精华液单| 久久久精品欧美日韩精品| 蜜桃久久精品国产亚洲av| 亚洲av日韩在线播放| 欧美性感艳星| 国产成人精品一,二区| av在线观看视频网站免费| 六月丁香七月| 亚洲精华国产精华液的使用体验| 日韩欧美三级三区| a级毛色黄片| 三级男女做爰猛烈吃奶摸视频| 有码 亚洲区| 国产精品蜜桃在线观看| 国产爱豆传媒在线观看| 国语对白做爰xxxⅹ性视频网站| 麻豆一二三区av精品| 午夜a级毛片| 水蜜桃什么品种好| 人体艺术视频欧美日本| 国产黄片视频在线免费观看| 少妇的逼水好多| 3wmmmm亚洲av在线观看| 热99re8久久精品国产| 看十八女毛片水多多多| 亚洲国产欧美人成| 日韩欧美国产在线观看| av免费观看日本| 亚洲四区av| 日韩成人伦理影院| 久久精品影院6| 一级av片app| 最近中文字幕2019免费版| 国产高清国产精品国产三级 | 亚洲精品国产成人久久av| 国产久久久一区二区三区| 国产免费一级a男人的天堂| 日韩成人av中文字幕在线观看| 最近最新中文字幕大全电影3| 久久欧美精品欧美久久欧美| 亚洲av免费高清在线观看| 村上凉子中文字幕在线| 中文乱码字字幕精品一区二区三区 | 亚洲成人av在线免费| 国产精品日韩av在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 我的老师免费观看完整版| 超碰av人人做人人爽久久| 老司机福利观看| 中国美白少妇内射xxxbb| 色尼玛亚洲综合影院| 永久网站在线| 国产私拍福利视频在线观看| 亚洲精品影视一区二区三区av| 久久久a久久爽久久v久久| 日韩av不卡免费在线播放| 亚洲国产精品成人综合色| 国产在视频线在精品| 亚洲av日韩在线播放| 国产高清国产精品国产三级 | 国产精品伦人一区二区| 免费看光身美女| 国产一区二区三区av在线| 超碰av人人做人人爽久久| 内地一区二区视频在线| 免费av不卡在线播放| 高清日韩中文字幕在线| 亚洲va在线va天堂va国产| 亚洲三级黄色毛片| 国产乱人视频| 免费观看人在逋| 特大巨黑吊av在线直播| 国产毛片a区久久久久| 国产成人91sexporn| 男人的好看免费观看在线视频| 国产成人精品一,二区| 人妻夜夜爽99麻豆av| 床上黄色一级片| 少妇人妻精品综合一区二区| 国产又色又爽无遮挡免| 久久国产乱子免费精品| 搡老妇女老女人老熟妇| 日本三级黄在线观看| 亚洲电影在线观看av| 久久久精品94久久精品| 少妇人妻精品综合一区二区| 国产精品99久久久久久久久| 在线观看一区二区三区| 男女视频在线观看网站免费| 内地一区二区视频在线| 在线观看66精品国产| 亚洲av福利一区| 最近的中文字幕免费完整| 欧美+日韩+精品| 少妇人妻一区二区三区视频| 免费观看性生交大片5| 超碰97精品在线观看| 黄色配什么色好看| 七月丁香在线播放| 熟女人妻精品中文字幕| 少妇的逼水好多| 亚洲国产精品sss在线观看| 一二三四中文在线观看免费高清| 亚洲不卡免费看| 国产精品1区2区在线观看.| av免费在线看不卡| 午夜精品在线福利| 国产成人a区在线观看| 边亲边吃奶的免费视频| 淫秽高清视频在线观看| 免费在线观看成人毛片| 精品少妇黑人巨大在线播放 | 特大巨黑吊av在线直播| 人体艺术视频欧美日本| 国产亚洲午夜精品一区二区久久 | 国产男人的电影天堂91| 国产美女午夜福利| 国产亚洲精品av在线| www日本黄色视频网| 在线观看美女被高潮喷水网站| 欧美高清性xxxxhd video| 国产爱豆传媒在线观看| 青春草亚洲视频在线观看| 青春草国产在线视频| 国产精品国产高清国产av| 日韩制服骚丝袜av| 最近手机中文字幕大全| 美女被艹到高潮喷水动态| 国产av在哪里看| 国产精品女同一区二区软件| 中文字幕免费在线视频6| 网址你懂的国产日韩在线| 久久婷婷人人爽人人干人人爱| 国产精华一区二区三区| 联通29元200g的流量卡| 国产三级在线视频| 99国产精品一区二区蜜桃av| av在线老鸭窝| 久久久亚洲精品成人影院| av在线播放精品| 直男gayav资源| 国产精品一区二区性色av| 我要搜黄色片| 久久久久久久久久久丰满| 久久久久久九九精品二区国产| 一级二级三级毛片免费看| 国产精品久久电影中文字幕| 最近中文字幕2019免费版| 国产精品日韩av在线免费观看| 晚上一个人看的免费电影| 中文亚洲av片在线观看爽| 国产一区二区在线av高清观看| 少妇人妻精品综合一区二区| 国产精品国产高清国产av| 久久久久久久久大av| 国产亚洲av嫩草精品影院| 国产精品野战在线观看| 深夜a级毛片| 国产免费视频播放在线视频 | 赤兔流量卡办理| 欧美日韩精品成人综合77777| 性色avwww在线观看| 赤兔流量卡办理| 久久这里只有精品中国| 亚洲国产精品成人久久小说| 久久这里有精品视频免费| 2021天堂中文幕一二区在线观| 免费搜索国产男女视频| 国产又黄又爽又无遮挡在线| 国产精品综合久久久久久久免费| 中国国产av一级| 亚洲精品乱码久久久v下载方式| 夜夜看夜夜爽夜夜摸| 激情 狠狠 欧美| 欧美极品一区二区三区四区| 永久网站在线| 嫩草影院精品99| 99久久成人亚洲精品观看| 九草在线视频观看| 亚洲一级一片aⅴ在线观看| 99热这里只有是精品在线观看| 国产精品一区二区三区四区免费观看| 日韩三级伦理在线观看| 韩国高清视频一区二区三区| 亚洲最大成人中文| 日韩制服骚丝袜av| 国产黄色视频一区二区在线观看 | 久久久久久久久大av| 亚洲欧美日韩高清专用| 久久久久免费精品人妻一区二区| 国产成人91sexporn| 国产视频内射| 欧美不卡视频在线免费观看| 五月玫瑰六月丁香| 韩国高清视频一区二区三区| 人体艺术视频欧美日本| 欧美一区二区亚洲| 亚洲久久久久久中文字幕| 亚洲精品aⅴ在线观看| 综合色丁香网| 国产综合懂色| 天天躁夜夜躁狠狠久久av| 精品免费久久久久久久清纯| 老女人水多毛片| 欧美日本亚洲视频在线播放| av国产久精品久网站免费入址| 啦啦啦观看免费观看视频高清| 国产精品无大码| 一区二区三区免费毛片| 可以在线观看毛片的网站| 97人妻精品一区二区三区麻豆| 最近手机中文字幕大全| 亚洲丝袜综合中文字幕| 免费播放大片免费观看视频在线观看 | 能在线免费看毛片的网站| 国产又色又爽无遮挡免| 五月伊人婷婷丁香| 亚洲精品国产av成人精品| 18+在线观看网站| 热99在线观看视频| 99热6这里只有精品| 黄色一级大片看看| 亚洲av福利一区| 国产亚洲一区二区精品| 国产极品精品免费视频能看的| 少妇裸体淫交视频免费看高清| 久久精品国产鲁丝片午夜精品| 在线观看av片永久免费下载| 卡戴珊不雅视频在线播放| 国产激情偷乱视频一区二区| 看非洲黑人一级黄片| 国内精品宾馆在线| 91精品一卡2卡3卡4卡| 少妇丰满av| 中文字幕精品亚洲无线码一区| 国产白丝娇喘喷水9色精品| 有码 亚洲区| 久久草成人影院| 中文天堂在线官网| 亚洲五月天丁香| 成人欧美大片| 亚洲欧洲国产日韩| 免费在线观看成人毛片| 日本黄色视频三级网站网址| 偷拍熟女少妇极品色| 日韩欧美国产在线观看| 自拍偷自拍亚洲精品老妇| 日韩高清综合在线| 日韩成人伦理影院| 亚洲国产欧洲综合997久久,| 九九爱精品视频在线观看| 精品人妻偷拍中文字幕| 岛国在线免费视频观看| 久久这里有精品视频免费| 久久婷婷人人爽人人干人人爱| 韩国av在线不卡| 精品国产露脸久久av麻豆 | 国产精品国产三级专区第一集| 男插女下体视频免费在线播放| 欧美日韩国产亚洲二区| 男女边吃奶边做爰视频| 亚洲av不卡在线观看| 久久精品人妻少妇| 男插女下体视频免费在线播放| 一个人看视频在线观看www免费| 亚洲一级一片aⅴ在线观看| 亚洲国产成人一精品久久久| 日韩av在线免费看完整版不卡| 久久久久九九精品影院| 伊人久久精品亚洲午夜| 亚洲熟妇中文字幕五十中出| 久久久久久久久久久丰满| 亚洲av日韩在线播放| 色哟哟·www| 国产又黄又爽又无遮挡在线| 国产不卡一卡二| 日韩 亚洲 欧美在线| 亚洲成人av在线免费| 国产黄色小视频在线观看| 夜夜爽夜夜爽视频| 色吧在线观看| 99热这里只有精品一区| 自拍偷自拍亚洲精品老妇| 国产精品一区www在线观看| 精品国产露脸久久av麻豆 | 久久久久精品久久久久真实原创| 欧美色视频一区免费| 免费av不卡在线播放| 69人妻影院| 国产伦理片在线播放av一区| 国产精品人妻久久久影院| 尤物成人国产欧美一区二区三区| 国产精品无大码| 免费看av在线观看网站| 22中文网久久字幕| 在现免费观看毛片| av福利片在线观看| 91在线精品国自产拍蜜月| 亚洲精华国产精华液的使用体验| 床上黄色一级片| 成人美女网站在线观看视频| 国产精品久久久久久精品电影| 99久国产av精品| 免费观看a级毛片全部| 成年免费大片在线观看| 久久久亚洲精品成人影院| 亚洲欧美日韩高清专用| 欧美三级亚洲精品| 一区二区三区免费毛片| 亚洲色图av天堂| 一夜夜www| 免费看光身美女| 免费观看人在逋| 亚洲精华国产精华液的使用体验| 国产精品一区二区三区四区久久| 精品人妻一区二区三区麻豆| 久久这里只有精品中国| 特大巨黑吊av在线直播| 国产 一区精品| 久久精品国产亚洲网站| 1000部很黄的大片| 精品无人区乱码1区二区| 全区人妻精品视频| 国语对白做爰xxxⅹ性视频网站| av福利片在线观看| 国产精品久久视频播放| 色综合色国产| 一二三四中文在线观看免费高清| 两性午夜刺激爽爽歪歪视频在线观看| 噜噜噜噜噜久久久久久91| 亚洲av中文字字幕乱码综合| 日韩欧美三级三区| 乱码一卡2卡4卡精品| 国产成人免费观看mmmm| 在线播放国产精品三级| 成人国产麻豆网| 亚洲不卡免费看| 国内揄拍国产精品人妻在线| 久久久国产成人精品二区| 2022亚洲国产成人精品| 免费黄色在线免费观看| 亚洲精品乱码久久久久久按摩| 亚洲国产精品成人久久小说| 天堂中文最新版在线下载 | 欧美激情国产日韩精品一区| 综合色丁香网| 国产免费男女视频| 激情 狠狠 欧美| 久久久久久久久中文| 狂野欧美激情性xxxx在线观看| 中文字幕av在线有码专区| 99在线人妻在线中文字幕| 精品国内亚洲2022精品成人| 国产精品一区二区三区四区久久| 十八禁国产超污无遮挡网站| 99久久中文字幕三级久久日本| 精品99又大又爽又粗少妇毛片| 亚洲精品成人久久久久久| 国产单亲对白刺激| 亚洲国产精品sss在线观看| 麻豆精品久久久久久蜜桃| 人人妻人人澡欧美一区二区| 青春草国产在线视频| 亚洲自拍偷在线| 干丝袜人妻中文字幕| 91av网一区二区| 中文字幕亚洲精品专区| 亚洲成人中文字幕在线播放| 久久精品国产99精品国产亚洲性色| 亚洲真实伦在线观看| 成人鲁丝片一二三区免费| 日本色播在线视频| 欧美成人免费av一区二区三区| 97人妻精品一区二区三区麻豆| 欧美成人精品欧美一级黄| 色噜噜av男人的天堂激情| 日本三级黄在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 最近中文字幕2019免费版| av国产久精品久网站免费入址| 在线天堂最新版资源| 在现免费观看毛片| 一区二区三区高清视频在线| 能在线免费观看的黄片| 男人舔奶头视频| 成人二区视频| 免费播放大片免费观看视频在线观看 | 午夜激情福利司机影院| 欧美日本亚洲视频在线播放| 免费看av在线观看网站| 亚洲成人精品中文字幕电影| 中文字幕制服av| 精品国产三级普通话版| 国产乱人偷精品视频| 免费电影在线观看免费观看| 丝袜喷水一区| 国产亚洲午夜精品一区二区久久 | 搡老妇女老女人老熟妇| 国产白丝娇喘喷水9色精品| 国产在视频线精品| 1000部很黄的大片| 国产中年淑女户外野战色| 人人妻人人澡欧美一区二区| 国产一级毛片在线| 亚洲欧美成人综合另类久久久 | 婷婷六月久久综合丁香| 在线播放国产精品三级| 国产成年人精品一区二区| 99久久精品一区二区三区| or卡值多少钱| 少妇熟女aⅴ在线视频| 六月丁香七月| 国产一级毛片七仙女欲春2| 能在线免费观看的黄片| 欧美日本亚洲视频在线播放| 成人二区视频| 久久99热这里只有精品18| 中文在线观看免费www的网站| 久久精品影院6| 在线播放无遮挡| 欧美激情在线99| 网址你懂的国产日韩在线| 国产亚洲精品久久久com| 天天躁夜夜躁狠狠久久av| 国产单亲对白刺激| 亚洲伊人久久精品综合 | 日韩欧美 国产精品| 麻豆精品久久久久久蜜桃| 黄片wwwwww| 中文在线观看免费www的网站| 国产高清视频在线观看网站| 可以在线观看毛片的网站| 国产精品精品国产色婷婷| 亚洲在久久综合| av卡一久久| 日本黄色视频三级网站网址| 午夜爱爱视频在线播放| 国产亚洲5aaaaa淫片| 91精品国产九色| 日日撸夜夜添| 国产探花在线观看一区二区| 狂野欧美激情性xxxx在线观看|