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

    祁連山國家公園植被覆蓋變化地形分異效應(yīng)

    2021-08-09 13:25:46龔純偉
    水土保持通報(bào) 2021年3期
    關(guān)鍵詞:坡向祁連山覆蓋度

    李 娟, 龔純偉

    (1.甘肅省基礎(chǔ)地理信息中心, 甘肅 蘭州 730000; 2.昆明理工大學(xué) 國土資源工程學(xué)院, 云南 昆明 650032)

    植被在生態(tài)系統(tǒng)中扮演著重要角色,是連接陸地土壤、水體和大氣的重要介質(zhì)[1],起著能量交換和物質(zhì)循環(huán)的作用,對(duì)人類活動(dòng)的強(qiáng)弱程度反應(yīng)十分敏感[2],可以作為區(qū)域生態(tài)環(huán)境變化研究的“指示劑”[3]。植被覆蓋的變化導(dǎo)致了陸地表層屬性隨之發(fā)生改變,進(jìn)而影響水土保持和氣候調(diào)節(jié)等,從而生態(tài)系統(tǒng)的穩(wěn)定性也發(fā)生改變[4]。植被覆蓋度是衡量地表植被狀況的一個(gè)最重要的指標(biāo),也是影響土壤侵蝕與水土流失的主要因子,對(duì)于區(qū)域環(huán)境變化和監(jiān)測(cè)研究具有重要意義[5]。其中,植被覆蓋度以及植被長勢(shì)的最佳指示因子歸一化植被指數(shù)(NDVI),被廣泛應(yīng)用于植被遙感中,是區(qū)域植被變化和生態(tài)環(huán)境變化的最有效監(jiān)測(cè)指標(biāo)[6-8]。

    祁連山是青藏高原與西北內(nèi)陸干旱半干旱地區(qū)的分界線,是溫帶大陸性氣候向高原山地氣候的過渡地帶,也是河西走廊地區(qū)重要的生態(tài)安全屏障和主要的水資源涵養(yǎng)地,是中國生態(tài)環(huán)境最脆弱的區(qū)域之一[9]。植被作為祁連山生態(tài)系統(tǒng)的重要組成部分,其變化不僅影響著當(dāng)?shù)厣鐣?huì)經(jīng)濟(jì)和生態(tài)環(huán)境,同時(shí)對(duì)中國西部生態(tài)安全有著極其重要的影響[10]。對(duì)于山地來說,地形對(duì)其植被分布具有重要影響,由于山地的垂直高度存在較大差異,氣溫、降水等氣候條件受其影響也呈現(xiàn)出空間差異性,因此,植被垂直地域分布呈現(xiàn)出明顯差異;山地坡度影響著地表徑流和水分以及熱因子的再分配,進(jìn)而改變土壤的特性及分布,植被的分布態(tài)勢(shì)也會(huì)受到影響;地勢(shì)的陡峭起伏度及坡度的緩急程度不同等,不但會(huì)形成局部小氣候的變化,而且對(duì)水土的流失與積聚都有影響;另外,在不同的地形地勢(shì)條件下,人類活動(dòng)的劇烈程度不同,所有這些因素都對(duì)植物的生長發(fā)育及空間分布狀況產(chǎn)生很大影響,因此,進(jìn)行地形因子與植被覆蓋度變化關(guān)系研究可以很好地揭示山地植被的地域分布規(guī)律[11-12]。近年來,關(guān)于植被覆蓋變化與地形因子關(guān)系的相關(guān)研究較多,如陳學(xué)兄等[13]研究了太原市城區(qū)植被覆蓋變化地形分異效應(yīng);熊俊楠等[14]研究了汶川地震災(zāi)區(qū)植被覆蓋變化與地形因子的關(guān)系;趙婷等[11]研究了秦嶺山地植被覆蓋變化地形分異效應(yīng);但針對(duì)祁連山區(qū)域的相關(guān)研究非常有限,大多以祁連山NDVI變化的驅(qū)動(dòng)力及對(duì)氣候變化的響應(yīng)為主,如陳京華等[15]研究了1982—2006年祁連山植被覆蓋的時(shí)空變化情況,認(rèn)為氣溫和水是影響祁連山植被NDVI變化的主要因子;邱麗莎等[16]研究了2000—2017年祁連山植被的時(shí)空變化,認(rèn)為18 a間祁連山植被總體呈改善趨勢(shì),植被變化除受氣候變化影響外,人工修復(fù)和改善也是關(guān)鍵因素;付建新等[17]從海拔、坡度和坡向的角度探討了祁連山南坡NDVI變化及其與氣溫和降水的關(guān)系,認(rèn)為不同海拔、坡度、坡向的植被受到氣溫和降水的影響不同。

    過去,由于人為活動(dòng)劇烈,大規(guī)模礦產(chǎn)資源的開發(fā),水電項(xiàng)目的建設(shè),對(duì)祁連山生態(tài)環(huán)境造成了極大的破壞,2014年以來,祁連山生態(tài)環(huán)境變化引起了國家高度重視,開始采取一系列措施進(jìn)行祁連山生態(tài)修復(fù),在此過程中,地形對(duì)祁連山植被覆蓋變化是否有影響?本文基于祁連山2006,2014,2019年3期遙感影像,采用像元二分法模型進(jìn)行植被覆蓋度遙感估算,在充分考慮地形因子面積差異影響的基礎(chǔ)上,探討了該區(qū)域不同時(shí)間植被覆蓋的變化情況以及不同高程、坡度、坡向等地形因子條件下祁連山植被覆蓋變化空間分布規(guī)律,旨在全面地了解不同時(shí)段內(nèi)祁連山植被覆蓋變化趨勢(shì)及其地形分異效應(yīng),以期為祁連山國家公園的植被恢復(fù)和重建提供科學(xué)依據(jù)。

    1 研究區(qū)概況

    祁連山國家公園位于西北內(nèi)陸干旱半干旱區(qū),地處青藏高原、蒙古高原和黃土高原交匯地帶,跨青海省東北部和甘肅省西部區(qū)域,西起當(dāng)金山口,東至烏鞘嶺,南接柴達(dá)木盆地,北臨河西走廊(圖1),東西長約1 000 km,南北寬約300 km,總面積約5.02×104km2。祁連山地勢(shì)西北高,東南低,地形多為山地,海拔在1 760~5 800 m左右,從西北到東南走向的山脈包括大雪山、托來南山、疏勒南山、走廊南山、冷龍嶺、烏鞘嶺等。祁連山自然條件復(fù)雜,水熱條件差異大,東部濕度大,降水多,西部氣候干燥,降水少,屬典型的高原大陸性氣候,是內(nèi)陸河流石羊河、黑河和疏勒河的發(fā)源地和徑流形成區(qū),其地理位置對(duì)河西綠洲的形成以及該地區(qū)生態(tài)環(huán)境的維持具有重要意義[18]。在地形地勢(shì)及氣候的共同影響下,祁連山植被呈現(xiàn)出獨(dú)特的垂直地域分布特征,海拔由低到高,植被分布依次為荒漠帶、山地草原帶、山地森林草原帶、高山灌叢草甸帶、高寒草甸和高山亞冰雪稀疏植被帶。祁連山森林和草原分界明顯,森林主要分布在陰坡,覆蓋率約為14.2%,植被類型有溫帶落葉闊葉林、溫帶常綠針葉林、寒溫性針葉林及針闊混交林等;草原主要分布在陽坡,是祁連山生態(tài)的重要組成部分,主要以灌木及草本植被為主,隨著海拔的升高,在亞高山區(qū)廣泛分布著高寒落葉灌叢和高山蒿草草甸及高寒草甸等植被,在高山區(qū)則以高寒墊狀植被及菊科高山植物為主。

    圖1 祁連山國家公園區(qū)位及地形圖

    2 數(shù)據(jù)來源及研究方法

    2.1 數(shù)據(jù)來源及處理

    本研究使用的遙感影像為2006年7—9月獲取的中巴地球資源衛(wèi)星影像(20 m×20 m),2014年8—10月獲取的美國Landsat8衛(wèi)星影像(30 m×30 m),2019年6—10月獲取的資源三號(hào)及高分一號(hào)衛(wèi)星影像(2 m×2 m)共3期。3期遙感影像均攝于植被生長較好的月份,完全能夠代表研究區(qū)域的土地利用和植被生長狀況,可用于比較不同期植被變化情況。為保證研究區(qū)域數(shù)據(jù)的完整性以及研究結(jié)果的準(zhǔn)確性,分別對(duì)三期影像進(jìn)行大氣校正、去除云影及水體[6]、幾何精校正(借助DEM數(shù)據(jù)消除地形引起的幾何畸變)和鑲嵌處理等,為驗(yàn)證影像校正精度,在研究區(qū)盡量均勻選取8個(gè)明顯的地面控制點(diǎn),利用GPS接收機(jī)確定其精確地理位置,并與影像上相應(yīng)點(diǎn)位地理坐標(biāo)進(jìn)行比較,經(jīng)統(tǒng)計(jì)分析,3期影像定位誤差均不超過一個(gè)像元,且同名地物點(diǎn)位置相同,滿足研究要求;此外,為保證不同期植被覆蓋度的遙感估算結(jié)果具有比較意義,對(duì)鑲嵌后的3期影像進(jìn)行重采樣(分辨率為10 m×10 m),利用祁連山國家公園范圍對(duì)3期影像進(jìn)行掩膜和裁剪,得到研究區(qū)范圍完全相同的3期影像數(shù)據(jù)。

    DEM數(shù)據(jù)來源于甘肅省自然資源廳,格網(wǎng)間距為25 m×25 m,將該數(shù)據(jù)進(jìn)行重采樣(格網(wǎng)間距為10 m×10 m),方便與植被覆蓋度數(shù)據(jù)進(jìn)行疊加分析。

    利用像元二分模型反演3期影像的植被覆蓋度。像元二分模型的原理是將圖像中的遙感信息(植被信息或波段)分解為綠色植被部分與無植被覆蓋部分所貢獻(xiàn)的信息兩部分,則任一像元的信息值可表示為兩部分信息的線性加權(quán)之和,該模型保留了植被覆蓋的相關(guān)信息,而植被類型、土壤背景、大氣的影響被大大削弱[19-20]。根據(jù)該原理,假設(shè)一個(gè)像元的信息可以分為土壤與植被兩部分,則植被覆蓋度f的計(jì)算公式可以表示為:

    式中:NDVIsoil,NDVIveg分別為無植被覆蓋或裸土區(qū)域、純植被覆蓋區(qū)域的NDVI值。根據(jù)研究區(qū)域整幅影像上NDVI灰度值的分布情況,以置信度0.5%取NDVI的最小值與最大值近似代表NDVIsoil和NDVIveg的值[21]。利用實(shí)地調(diào)查與地面測(cè)量相結(jié)合的方式,大致均勻選取10個(gè)點(diǎn)位,利用林地冠層分析儀對(duì)點(diǎn)位及其鄰域進(jìn)行5~10次測(cè)量并取其平均值作為驗(yàn)證點(diǎn)的植被覆蓋度測(cè)量結(jié)果,以此對(duì)相應(yīng)點(diǎn)位遙感影像植被覆蓋度反演結(jié)果進(jìn)行精度定量評(píng)價(jià)(圖2)。

    圖2 祁連山驗(yàn)證點(diǎn)植被覆蓋度測(cè)量值與反演值結(jié)果對(duì)比

    2.2 研究方法

    2.2.1 植被覆蓋度變化趨勢(shì)分析 趨勢(shì)分析是采用最小二乘法逐像元對(duì)一組隨時(shí)間變化的序列數(shù)據(jù)采用趨勢(shì)線進(jìn)行回歸分析來擬合植被覆蓋度的年際變化趨勢(shì),其公式為:

    (1)

    式中:θslope為研究時(shí)段內(nèi)植被覆蓋度的變化趨勢(shì); NDVIi為第i年的植被覆蓋度;a為研究年份數(shù)目。

    對(duì)植被覆蓋變化趨勢(shì)率進(jìn)行計(jì)算,公式為:

    (2)

    按顯著性水平0.1對(duì)NDVI變化趨勢(shì)率進(jìn)行t檢驗(yàn),結(jié)合θslope值,將植被覆蓋度變化趨勢(shì)分為5類:顯著增加型(θslope>0,t<0.01);輕微增加型(θslope>0,0.010.05);顯著減少型(θslope<0,t<0.01);輕微減少型(θslope<0,0.01

    2.2.2 地形面積修正 當(dāng)?shù)匦巫兓瘯r(shí),會(huì)出現(xiàn)在某種地形條件下某種植被覆蓋變化類型面積很小,但其面積比例相對(duì)于整個(gè)研究區(qū)卻很大,從而影響植被覆蓋度評(píng)價(jià)結(jié)果的合理性[22]。因此,在分析植被覆蓋變化類型的空間分布狀況時(shí),采用地形面積修正系數(shù)v來消除地形的影響,其計(jì)算公式如下:

    (3)

    式中:Svt為在某地形條件t下某種植被變化類型v所占面積;Sv為某種植被變化類型v所占面積;St指某種地形t的面積,S總為整個(gè)研究區(qū)域的總面積。v>1,表示植被變化類型v在地形條件t下面積變化趨勢(shì)為增加;v=1,表示植被變化類型v在地形條件t下面積變化趨于平穩(wěn);v<1,表示植被變化類型v在地形條件t下面積變化趨勢(shì)為減少[23]。

    2.3 高程、坡度、坡向分級(jí)

    基于ArcGIS空間分析模塊,利用DEM數(shù)據(jù)生成高程、坡度、坡向分級(jí)圖。根據(jù)研究區(qū)域地形、地理狀況,并結(jié)合已有研究成果[17],將研究區(qū)高程分為<2 200,2 200~2 700,2 700~3 200,3 200~3 700,3 700~4 200,4 200~4 700 m和>4 700 m共7級(jí);坡度劃分依據(jù)自然資源部頒布的《第三次全國國土調(diào)查技術(shù)規(guī)程》,再結(jié)合祁連山地形狀況將坡度分為≤5°,5°~10°,10°~15°,15°~25°,25°~40°和>40°共6級(jí);將坡向劃分為平地、陽坡(135°~225°)、半陽坡(225°~315°)、陰坡(315°~45°)、半陰坡(45°~135°)共5類。坡度、坡向空間分布見圖3。

    圖3 祁連山高程、坡度、坡向空間分布特征

    2.4 植被覆蓋度分級(jí)

    根據(jù)水利部2008年頒布的《土壤侵蝕分類分級(jí)標(biāo)準(zhǔn)》(SL190-2007),將植被覆蓋度的不同等級(jí)與水土流失的不同等級(jí)相對(duì)應(yīng),同時(shí)結(jié)合祁連山地理狀況和地面調(diào)查情況,將研究區(qū)域植被覆蓋度分為Ⅵ級(jí):高植被覆蓋度(70%~100%);Ⅴ級(jí):較高植被覆蓋度(50%~70%);Ⅳ級(jí):中等植被覆蓋度(30%~50%);Ⅲ級(jí):較低植被覆蓋度(10%~30%);Ⅱ級(jí):低植被覆蓋(0%~10%);Ⅰ級(jí):無植被覆蓋或裸土(0)這6個(gè)等級(jí)[20-21]。祁連山3期植被覆蓋度分級(jí)見圖4。

    圖4 祁連山2006,2014,2019年植被覆蓋度及其分級(jí)

    3 結(jié)果與分析

    3.1 祁連山植被覆蓋度時(shí)空分布及變化特征

    3.1.1 祁連山植被覆蓋度時(shí)空分布特征 根據(jù)圖4可知,祁連山植被覆蓋度空間分布特征為“西北部低,東南部高”,Ⅵ級(jí)植被覆蓋區(qū)主要分布在烏鞘嶺的東南部區(qū)域;Ⅴ級(jí)植被覆蓋區(qū)主要分布在走廊南山東南部、冷龍嶺、烏鞘嶺及托來山東南部等區(qū)域,主要為森林草原覆蓋區(qū);Ⅳ級(jí)植被覆蓋區(qū)主要分布在托來山及走廊南山的西北部,主要為草原和高山植被混合覆蓋區(qū)域;荒漠與草原植被過渡區(qū)域?yàn)棰蠹?jí)植被覆蓋區(qū)域;祁連山西北部的大雪山、疏勒南山、黨河南山等區(qū)域主要為裸地、沙漠、碎石和冰川區(qū),植被覆蓋度為Ⅱ級(jí)和Ⅰ級(jí)。對(duì)祁連山各級(jí)植被覆蓋度進(jìn)行統(tǒng)計(jì)(表1)可知,2006年和2019年,祁連山植被覆蓋度均以Ⅲ級(jí)為主,分別占46.81%和42.7%,但2019年面積占比比2006年減少4.11%;Ⅱ級(jí)植被覆蓋度由2006年的14.95%減少至2 019年的12.09%,減少約2.86%;其他類型植被覆蓋度均有所增加,尤其以Ⅳ級(jí)植被覆蓋度區(qū)域面積增加最明顯,增幅約為3.18%,Ⅴ級(jí)和Ⅵ級(jí)植被覆蓋度增幅分別為1.07%和1.95%。但需要關(guān)注的是,祁連山無植被覆蓋和裸地(Ⅰ級(jí))面積略有增加,增加面積比例約為0.77%。

    表1 2006,2019年祁連山各級(jí)植被覆蓋度面積比例及其變化情況

    3.1.2 祁連山植被覆蓋度時(shí)空變化特征 對(duì)祁連山植被覆蓋度變化趨勢(shì)進(jìn)行分析(圖5),結(jié)果表明:2006—2019年祁連山植被覆蓋度整體的變化趨勢(shì)為增加,植被覆蓋度顯著增加和輕微增加區(qū)域面積分別占24.9%和21.8%,主要分布在黨河南山和疏勒南山西北部及烏鞘嶺東南部區(qū)域;植被覆蓋減少區(qū)域主要分布在疏勒南山東南部、走廊南山中部和托來山、大雪山等的部分區(qū)域,顯著減少和輕微減少區(qū)域面積分別占10.8%,22.5%。其中,2006—2014年,祁連山植被覆蓋度平均變化率為-0.016,總體呈現(xiàn)減少趨勢(shì),經(jīng)顯著性檢驗(yàn),植被覆蓋度顯著減少和輕微減少區(qū)域面積分別占總面積的18.2%和37.1%,主要分布在祁連山西北部大部分區(qū)域、走廊南山張掖段、托來山、托來南山及疏勒南山東南部等的部分區(qū)域;24.7%的區(qū)域未通過顯著性檢驗(yàn),為穩(wěn)定區(qū);輕微增加和顯著增加區(qū)域面積分別占14.6%和5.4%,主要分布在黨河南山、烏鞘嶺及冷龍嶺的部分區(qū)域。2014—2019年,祁連山植被覆蓋度平均變化率為0.028,總體呈現(xiàn)出增加趨勢(shì),植被處于恢復(fù)期,經(jīng)顯著性檢驗(yàn),植被覆蓋度顯著增加和輕微增加區(qū)域面積分別占總面積的27.4%和37.6%,主要分布在黨河南山、疏勒南山西北部及托來山北部、走廊南山、烏鞘嶺東南部、冷龍嶺的部分區(qū)域;未通過顯著性檢驗(yàn)的穩(wěn)定區(qū)域面積占23.8%;輕微減少和顯著減少區(qū)域面積分別占9.1%和2.1%,主要分布在走廊南山張掖段以及烏鞘嶺中部的個(gè)別區(qū)域。

    圖5 2006-2019年祁連山不同時(shí)段植被覆蓋度變化類型分布

    3.2 祁連山植被覆蓋度地形分異特征

    3.2.1 祁連山植被覆蓋度高程梯度分異 山地的氣候條件及土壤狀況等自然因素的變化由海拔高程決定,隨著海拔的升高,形成了明顯的植被垂直帶,不同植被垂直帶對(duì)人為活動(dòng)干擾以及氣候變化的響應(yīng)度有所不同,從而表現(xiàn)出不同的變化趨勢(shì)。圖6a為祁連山植被覆蓋度垂直分布狀況,由圖可知,隨著海拔的升高,祁連山植被覆蓋度整體呈現(xiàn)出先增大后減小的 “倒Ⅴ型”趨勢(shì)。在低海拔區(qū)域(2 200 m以下),植被生長所需的水熱條件充足,但受人為活動(dòng)干擾較大,植被覆蓋度較低,地表呈現(xiàn)出荒漠草原景觀;海拔升高,對(duì)人類活動(dòng)的限制逐步增大,植被生長干擾程度逐漸減小,加之水熱條件較適宜,植被覆蓋逐漸好轉(zhuǎn),在海拔2 200~3 200 m的中山區(qū),植被類型主要為森林草原;海拔3 200~3 700 m的亞高山區(qū),生境較濕潤,植被類型主要為灌叢草原;海拔3 700~4 200 m之間是植被覆蓋度最高的區(qū)域,年降水量較高,但氣溫較低,呈現(xiàn)出草甸草原景觀;隨著海拔的繼續(xù)升高,接收到的太陽輻射強(qiáng)度、空氣流動(dòng)等氣象條件也隨之發(fā)生改變,海拔升高至4 200 m以上,氣溫驟降,植被生長受熱量條件和陡峭地勢(shì)的限制,植被的適宜性變得較為脆弱,植被覆蓋度降低,植被類型為高寒草甸和高山亞冰雪稀疏植被。

    3.2.2 祁連山植被覆蓋度在坡度坡向上的空間分異 山地坡度影響著地表徑流和水分以及熱因子的再分配,進(jìn)而改變土壤的特性及分布,同時(shí)也會(huì)影響植被的分布態(tài)勢(shì)。圖6b為植被覆蓋度在坡度上的分布狀況,隨著坡度的增加,祁連山植被覆蓋度呈現(xiàn)出先降低后增高最后又降低的倒“N”趨勢(shì)。坡度在5°以下,地勢(shì)平緩,水土不易流失,植被覆蓋度相對(duì)較高,但祁連山區(qū)平地很少;在5°~15°相對(duì)比較平緩的區(qū)域,由于人為活動(dòng)的干擾,植被覆蓋度較低;隨著坡度的升高,人為活動(dòng)減少,植被覆蓋回歸正常趨勢(shì),植被覆蓋度增加;坡度>40°的區(qū)域,水土流失嚴(yán)重,并且得到的積溫較少,不利于植物生長,植被覆蓋度急劇降低。

    坡向即海拔改變量變化的方向。不同的坡向,接收的太陽輻射量及水分蒸散發(fā)量均存在差異,從而影響植被覆蓋的空間分布。由圖6c可知,祁連山陰坡由于日照的強(qiáng)度較低、獲得的太陽輻射總量相對(duì)較少,水分蒸散發(fā)量較低,相對(duì)于其他坡向更有利于植被生長,因此植被覆蓋度面積占比相對(duì)較高,達(dá)到29.1%;其他坡向植被覆蓋度面積占比差別不大,約為23%左右;平地由于人為活動(dòng)劇烈,植被覆蓋度很低。

    圖6 祁連山植被覆蓋度在不同地形條件下的空間分布特征

    3.3 祁連山植被覆蓋度變化地形分異特征分析

    3.3.1 高程梯度對(duì)祁連山植被覆蓋度變化的影響 將祁連山植被覆蓋度變化類型數(shù)據(jù)與高程分級(jí)數(shù)據(jù)進(jìn)行疊加,統(tǒng)計(jì)出不同高程范圍內(nèi)的植被變化類型,并計(jì)算與之對(duì)應(yīng)的地形面積改正系數(shù)v,詳細(xì)結(jié)果如表2所示。

    表2 祁連山2006-2019年不同植被覆蓋度變化類型v值垂直變化

    由表2可知,不同高程范圍內(nèi)各植被變化類型的空間分異較明顯。2006—2019年整個(gè)時(shí)段內(nèi),植被覆蓋顯著增加型隨高程的升高呈明顯減少趨勢(shì),在3 200 m以下區(qū)域?yàn)閮?yōu)勢(shì)分布(v>1.54),在2 200 m以下植被增加趨勢(shì)尤其明顯(v=2.99);輕微增加型先增加后減少,在3 200 m以下區(qū)域均為優(yōu)勢(shì)分布,在2 200~2 700 m之間區(qū)域優(yōu)勢(shì)分布較明顯(v=1.62);3 200~3 700 m之間分布相對(duì)穩(wěn)定(v略大于或略小于(1) ;顯著減少型隨高程升高先增加后逐漸減少,在3 700~4 200 m之間區(qū)域呈較弱優(yōu)勢(shì)分布(v=1.2);輕微減少型隨高程升高不斷增加,在4 200 m以上區(qū)域?yàn)橹鲗?dǎo)類型,在4 700 m以上分布優(yōu)勢(shì)明顯(v=2.64)。其中,2006—2014年,植被覆蓋顯著減少型隨高程升高逐漸減少,在4 200 m以下區(qū)域均為優(yōu)勢(shì)分布(v>1.0),尤其在2 200 m以下減少趨勢(shì)顯著(v=2.16);輕微減少型隨高程升高呈上下波動(dòng)趨勢(shì),在3 700 m以下及4 200 m以上區(qū)域均為優(yōu)勢(shì)分布,在2 200 m以下優(yōu)勢(shì)分布尤其明顯(v=3.38);顯著增加型及輕微增加型隨高程增加先增加后減少,在各高程范圍內(nèi)均為非優(yōu)勢(shì)分布(v<1.0)。2014—2019年,顯著增加型隨高程增加逐漸減少,在3 200 m以下區(qū)域呈優(yōu)勢(shì)分布,在2 200 m以下區(qū)域優(yōu)勢(shì)分布顯著(v=2.48);輕微增加型先增加后減少,在3 700 m以下區(qū)域?yàn)閮?yōu)勢(shì)分布,尤其在2 200~2 700 m范圍內(nèi)分布特別顯著(v=4.34),到3 700 m后v值迅速降低至1.0以下,呈非優(yōu)勢(shì)分布;顯著減少型呈先增加后減少趨勢(shì),輕微減少型呈逐漸增加趨勢(shì),兩者在各高程范圍內(nèi)均為非優(yōu)勢(shì)分布。

    3.3.2 坡度對(duì)祁連山植被覆蓋度變化的影響 祁連山地形復(fù)雜多樣,植被覆蓋變化受坡度影響較明顯。將祁連山植被覆蓋度變化數(shù)據(jù)與坡度分級(jí)數(shù)據(jù)疊加,統(tǒng)計(jì)出不同年份不同坡度范圍內(nèi)植被覆蓋變化情況,見表3。由表3可知,植被覆蓋度變化隨坡度的增加存在不同程度的差異,從2006—2019年整個(gè)時(shí)段來看,顯著減少型和輕微減少型在坡度25°以上區(qū)域均為優(yōu)勢(shì)分布,尤其在坡度40°以上區(qū)域優(yōu)勢(shì)分布明顯,v值分別為1.33,1.27;輕微增加型和顯著增加型在坡度15°以下區(qū)域占主導(dǎo)類型,在坡度5°以下區(qū)域分布優(yōu)勢(shì)尤其顯著,v值分別為1.34,1.37;15°~25°坡度范圍內(nèi),各植被覆蓋變化類型分布相對(duì)穩(wěn)定(v值略大于或略小于1) 。其中,2006—2014年,植被覆蓋顯著減少型在各坡度范圍內(nèi)均為優(yōu)勢(shì)分布,在坡度25°以上區(qū)域優(yōu)勢(shì)分布顯著(v>1.47);輕微減少型在坡度5°以上區(qū)域均為微弱優(yōu)勢(shì)分布,在坡度40°以上區(qū)域分布相對(duì)明顯(v=1.10);顯著增加類型在各坡度范圍內(nèi)均為非優(yōu)勢(shì)分布(v<1.0),輕微增加型僅在坡度5°以下區(qū)域呈微弱優(yōu)勢(shì)分布(v=1.07)。2014—2019年,輕微增加和顯著增加型隨坡度增加均逐漸減少,輕微增加型在坡度25°以下區(qū)域?yàn)閮?yōu)勢(shì)分布(1.03

    表3 祁連山2006-2019年不同植被覆蓋度變化類型v值隨坡度變化情況

    3.3.3 坡向?qū)ζ钸B山植被覆蓋度的影響 將祁連山植被覆蓋變化類型數(shù)據(jù)與坡向分類數(shù)據(jù)進(jìn)行疊加分析,統(tǒng)計(jì)出不同年份不同坡向植被覆蓋度變化情況見表4。

    由表4可知,除平地外,祁連山各植被覆蓋度變化類型在其他坡向分布差異不顯著。2006—2019年,在平地區(qū)域,植被覆蓋顯著減少型占主導(dǎo)類型(v=1.6),表明該區(qū)域受人為活動(dòng)干擾較大,植被覆蓋度較低;在陰坡區(qū)域,顯著增加型呈較弱優(yōu)勢(shì)分布(v=1.16);在陽坡區(qū)域,輕微增加型呈較弱優(yōu)勢(shì)分布(v=1.14);除此之外,各植被變化類型在其他坡向分布差異不明顯。其中,2006—2014年,在平地區(qū)域,植被覆蓋顯著減少型優(yōu)勢(shì)分布明顯(v=2.91),其他植被覆蓋變化類型呈非優(yōu)勢(shì)分布;在其他坡向,各植被覆蓋變化類型分布差異不明顯(v值略大于或略小于1)。2014—2019年,在平地區(qū)域,顯著增加型和輕微增加型均呈優(yōu)勢(shì)分布,v值分別為1.49,1.18,表明該區(qū)域植被覆蓋在逐漸恢復(fù);在其他坡向,各植被覆蓋變化類型分布差異不明顯。

    表4 祁連山2006-2019年不同植被覆蓋度變化類型v值隨坡向變化情況

    4 討論與結(jié)論

    4.1 結(jié) 論

    針對(duì)祁連山植被覆蓋度研究,已有成果多基于氣候?qū)ζ渥兓捻憫?yīng),從地形角度出發(fā)的研究較少。本文基于祁連山2006,2014,2019年3期影像,采用像元二分法模型進(jìn)行植被覆蓋度估算,在充分考慮地形面積差異影響的基礎(chǔ)上,探討了該區(qū)域不同時(shí)間植被覆蓋變化情況以及不同高程、坡度、坡向等地形條件下植被覆蓋的空間分異性,研究結(jié)果對(duì)祁連山生態(tài)修復(fù)具有參考意義。

    (1) 祁連山植被覆蓋度空間分布特征為“西北部低,東南部高”,總體以較低植被覆蓋度為主。2006—2019年祁連山植被覆蓋度整體的變化趨勢(shì)為增加,增加區(qū)域主要分布在黨河南山和疏勒南山西北部及烏鞘嶺東南部區(qū)域,尤其在祁連山西北部部分區(qū)域植被覆蓋增加顯著,增加面積約占46.7%,減少面積約占33.3%,植被恢復(fù)狀況較好。低、較低植被覆蓋度面積減少;中、較高和高植被覆蓋度面積均有增加。其中,2006—2014年,祁連山植被覆蓋度平均變化率為-0.016,總體呈現(xiàn)減少趨勢(shì);2014—2019年,祁連山植被覆蓋度平均變化率為0.028,總體呈現(xiàn)出增加趨勢(shì)。

    (2) 祁連山植被覆蓋度隨海拔的升高呈現(xiàn)出先增加后減少的趨勢(shì),其中,海拔3 700~4 200 m之間是植被覆蓋度最高的區(qū)域。從不同高程范圍來看,各植被覆蓋變化類型的空間分異較明顯,海拔越高,植被覆蓋度減少越明顯。2006—2019年,3 200 m以下的中低海拔區(qū)域植被覆蓋呈增加趨勢(shì),在2 200 m以下的低海拔區(qū)域的荒漠草原帶植被增加趨勢(shì)尤其明顯,根據(jù)相關(guān)資料,植被增加類型主要以草本植物及灌木為主;海拔3 700 m以上的中高海拔區(qū)域植被覆蓋則呈較弱的減少趨勢(shì),在4 700 m以上高海拔區(qū)域的高寒草甸帶植被減少趨勢(shì)更明顯,3 200~3 700 m高程范圍內(nèi)植被覆蓋變化相對(duì)穩(wěn)定。

    (3) 從坡度來看,植被覆蓋度變化情況與坡度的大小呈反比,坡度越大,植被覆蓋度減少趨勢(shì)越顯著。2006—2019年,在坡度15°以下區(qū)域植被覆蓋呈增加趨勢(shì),在坡度5°以下增加優(yōu)勢(shì)尤其顯著;在坡度25°以上區(qū)域植被覆蓋呈減少趨勢(shì),在坡度40°以上區(qū)域減少趨勢(shì)尤其明顯;在15°~25°坡度范圍內(nèi),各植被覆蓋變化類型為穩(wěn)定分布。

    (4) 從坡向來看,除平地外,各植被覆蓋變化類型在其他坡向分布差異不顯著。平地區(qū)域由于易受人為活動(dòng)的干擾,植被覆蓋度呈減少趨勢(shì);其他坡向?qū)χ脖桓采w變化影響不大。

    4.2 討 論

    祁連山作為西北地區(qū)生態(tài)環(huán)境最脆弱的區(qū)域之一,近年來,在全球氣候變暖的大背景下,其植被生長環(huán)境也呈現(xiàn)出較為明顯的變化,氣溫、降水等氣候條件受地形影響呈現(xiàn)出空間差異性,同時(shí),人類活動(dòng)的劇烈程度在不同的地形地勢(shì)條件下也不同,所有這些因素都對(duì)植物的生長發(fā)育及空間分布狀況產(chǎn)生很大影響,因此,在進(jìn)行祁連山生態(tài)保護(hù)及修復(fù)的同時(shí)關(guān)注地形對(duì)植被覆蓋變化的影響具有重要意義。

    本文研究結(jié)果表明,2006—2019年,祁連山植被覆蓋度總體上呈增加趨勢(shì),這是由于2014年以來,政府加強(qiáng)對(duì)祁連山地區(qū)的生態(tài)保護(hù)及生態(tài)環(huán)境治理修復(fù)工作,采取封山育林、禁牧、禁止濫砍濫伐、濫采濫挖等嚴(yán)格的生態(tài)保護(hù)和恢復(fù)措施,人類活動(dòng)大大減少,植被逐漸恢復(fù);另外,隨著全球氣溫升高,祁連山冰川積雪消融加速,也對(duì)植被生長發(fā)育有一定的促進(jìn)作用。但需要特別說明的是,祁連山中等、較高和高植被覆蓋度增加的同時(shí),疏勒南山東南部、走廊南山和大雪山等個(gè)別區(qū)域的無植被覆蓋或裸露土地也略有增加,分析其原因,可能是由于該區(qū)域過去人為活動(dòng)劇烈,礦產(chǎn)資源的開發(fā),水電項(xiàng)目的建設(shè),對(duì)自然生態(tài)環(huán)境造成了極大的破壞,植被的適應(yīng)性變得特別脆弱,導(dǎo)致植被覆蓋出現(xiàn)減少現(xiàn)象,具體的原因還需要進(jìn)一步的研究和探討才能確定。

    從祁連山生態(tài)植被變化的地形分布特征來看,其在高程、坡度上的差異明顯,植被覆蓋度隨高程和坡度增加呈減少趨勢(shì)。海拔3 700 m以上區(qū)域植被覆蓋減少,該區(qū)域雖然地勢(shì)高,人類活動(dòng)少,但植被類型主要以脆弱的灌草叢及高山草甸、高寒墊狀植被類型為主,脆弱的本底條件以及寒冷的氣候等自然因素可能是造成植被退化的主要原因,因而推斷氣候變化對(duì)祁連山高海拔區(qū)域的植被覆蓋影響較大,在日后的環(huán)境生態(tài)治理工作中需加強(qiáng)對(duì)高海拔區(qū)域的脆弱的草甸等植被類型的保護(hù)和管理;另外,坡度25°以上區(qū)域植被覆蓋呈減少趨勢(shì),究其原因,可能是由于該區(qū)域地勢(shì)陡峭,接收的太陽輻射強(qiáng)度較低,水土流失嚴(yán)重,植被生長適宜性脆弱,植被狀況容易發(fā)生變化,有必要引起足夠的重視,在進(jìn)行生態(tài)保護(hù)的同時(shí)加強(qiáng)該區(qū)域的水土涵養(yǎng)能力有利于植被覆蓋狀況的改善;需要特別說明的是,在分析祁連山植被變化在坡向上的差異時(shí)發(fā)現(xiàn),雖然植被覆蓋變化在不同坡向上存在一定差異,但變化幅度不明顯,可能是由于祁連山地勢(shì)走向?yàn)閺奈鞅钡綎|南,且整體海拔較高,地形狀況復(fù)雜,獨(dú)特的地形地勢(shì)條件形成了局地的小氣候,水熱條件組合情況多變,導(dǎo)致植被覆蓋變化受坡向的影響不大,而受高程和坡度的影響較大,不過這只是一種推斷,具體的原因和影響效應(yīng)還需在后續(xù)的工作中進(jìn)一步探討。

    此外,植被覆蓋度變化的影響因素較多,本文只從地形因子角度出發(fā)進(jìn)行了植被覆蓋度的變化分析,沒有考慮氣溫、降水等氣象因子以及政策等不確定因素所帶來的影響,在后續(xù)的研究中,綜合氣候、人類活動(dòng)等影響因素探討植被覆蓋度變化是非常必要的。

    猜你喜歡
    坡向祁連山覆蓋度
    圖志
    發(fā)展(2024年1期)2024-04-25 03:45:30
    呼和浩特市和林格爾縣植被覆蓋度變化遙感監(jiān)測(cè)
    基于NDVI的晉州市植被覆蓋信息提取
    祁連山下
    低覆蓋度CO分子在Ni(110)面的吸附研究
    DEM地表坡向變率的向量幾何計(jì)算法
    祁連山草原:如夢(mèng)如幻近高天
    祁連山
    黃河之聲(2016年20期)2016-02-21 11:55:33
    青藏高原東緣高寒草甸坡向梯度上植物光合生理特征研究
    土壤與作物(2015年3期)2015-12-08 00:46:55
    坡向和坡位對(duì)小流域梯田土壤有機(jī)碳、氮變化的影響
    久久性视频一级片| 一本一本久久a久久精品综合妖精| 精品人妻熟女毛片av久久网站| 老司机影院毛片| 视频区图区小说| 久久99热这里只频精品6学生| 亚洲国产精品成人久久小说| 国产精品久久久久久久久免| 天天躁夜夜躁狠狠久久av| 国产精品免费大片| 免费在线观看完整版高清| 麻豆乱淫一区二区| 哪个播放器可以免费观看大片| 国产乱来视频区| 精品一区在线观看国产| 亚洲精品日韩在线中文字幕| 欧美日韩亚洲高清精品| 午夜福利视频精品| 亚洲精品日韩在线中文字幕| 久久亚洲国产成人精品v| 黄色 视频免费看| 青青草视频在线视频观看| 青草久久国产| www.av在线官网国产| 男人舔女人的私密视频| 中文乱码字字幕精品一区二区三区| 欧美日韩视频高清一区二区三区二| 777米奇影视久久| 亚洲成人手机| 9色porny在线观看| 日日爽夜夜爽网站| 国产成人av激情在线播放| 久久 成人 亚洲| 久久国产亚洲av麻豆专区| 精品卡一卡二卡四卡免费| 老司机深夜福利视频在线观看 | 国产xxxxx性猛交| 丝瓜视频免费看黄片| 超碰成人久久| 99精国产麻豆久久婷婷| 精品人妻一区二区三区麻豆| 人人妻人人添人人爽欧美一区卜| 久久久久精品久久久久真实原创| 欧美日韩综合久久久久久| 免费黄频网站在线观看国产| 97人妻天天添夜夜摸| 国产视频首页在线观看| 国产精品 国内视频| 18禁国产床啪视频网站| 蜜桃国产av成人99| 精品国产一区二区久久| 天天添夜夜摸| 男人添女人高潮全过程视频| 777久久人妻少妇嫩草av网站| 精品一区二区三区av网在线观看 | 欧美成人精品欧美一级黄| 国产成人a∨麻豆精品| 一级毛片 在线播放| 一级片免费观看大全| 不卡视频在线观看欧美| 伊人久久国产一区二区| 男人爽女人下面视频在线观看| 青春草亚洲视频在线观看| 国产欧美日韩综合在线一区二区| 免费女性裸体啪啪无遮挡网站| 国产又色又爽无遮挡免| 亚洲精品国产av蜜桃| 国产亚洲一区二区精品| 在线观看免费高清a一片| 少妇的丰满在线观看| 只有这里有精品99| 亚洲人成77777在线视频| 免费少妇av软件| 久久久久国产精品人妻一区二区| 日本爱情动作片www.在线观看| 少妇人妻久久综合中文| 国产成人精品无人区| 久久精品人人爽人人爽视色| 久久狼人影院| 999久久久国产精品视频| 免费日韩欧美在线观看| 中文字幕av电影在线播放| 欧美少妇被猛烈插入视频| 性少妇av在线| 99香蕉大伊视频| 亚洲精品自拍成人| 老司机影院毛片| 深夜精品福利| 欧美在线黄色| 悠悠久久av| 国产亚洲av高清不卡| 不卡av一区二区三区| 9191精品国产免费久久| 欧美久久黑人一区二区| 男女边吃奶边做爰视频| 成人国产av品久久久| av有码第一页| 亚洲精品久久成人aⅴ小说| 亚洲美女搞黄在线观看| 在线观看www视频免费| 亚洲一区中文字幕在线| 久久精品久久精品一区二区三区| 欧美国产精品一级二级三级| 在线观看三级黄色| 国产精品秋霞免费鲁丝片| 欧美日韩成人在线一区二区| 大话2 男鬼变身卡| 亚洲国产日韩一区二区| 国产精品国产av在线观看| 国产精品香港三级国产av潘金莲 | 欧美日韩av久久| av国产久精品久网站免费入址| 黄片播放在线免费| 亚洲国产中文字幕在线视频| 亚洲欧美成人综合另类久久久| 汤姆久久久久久久影院中文字幕| 国产精品一区二区在线观看99| 久久国产精品男人的天堂亚洲| 亚洲精华国产精华液的使用体验| 啦啦啦 在线观看视频| 国产熟女欧美一区二区| 国产成人一区二区在线| 黄片小视频在线播放| 国产不卡av网站在线观看| 高清在线视频一区二区三区| 91成人精品电影| 久久精品国产a三级三级三级| 熟妇人妻不卡中文字幕| 天堂俺去俺来也www色官网| 欧美日韩av久久| 热99国产精品久久久久久7| 在线观看免费高清a一片| 黄色怎么调成土黄色| 国产99久久九九免费精品| 黑人欧美特级aaaaaa片| 久久久久精品国产欧美久久久 | 亚洲成av片中文字幕在线观看| av视频免费观看在线观看| 热99久久久久精品小说推荐| 国产精品国产av在线观看| 国产不卡av网站在线观看| 少妇被粗大的猛进出69影院| 日韩视频在线欧美| 精品亚洲乱码少妇综合久久| 欧美亚洲日本最大视频资源| 777久久人妻少妇嫩草av网站| 老司机深夜福利视频在线观看 | 热99久久久久精品小说推荐| 搡老岳熟女国产| 在线亚洲精品国产二区图片欧美| 免费黄色在线免费观看| 久久久久精品性色| 99久国产av精品国产电影| 男男h啪啪无遮挡| 欧美少妇被猛烈插入视频| xxx大片免费视频| 亚洲国产av影院在线观看| 黄片无遮挡物在线观看| 天天操日日干夜夜撸| 亚洲成人av在线免费| 日本黄色日本黄色录像| 国产成人午夜福利电影在线观看| 国产色婷婷99| 美女午夜性视频免费| 人人妻,人人澡人人爽秒播 | 亚洲国产中文字幕在线视频| 桃花免费在线播放| 亚洲av成人精品一二三区| 日本av免费视频播放| 亚洲av国产av综合av卡| 美女大奶头黄色视频| 国产成人精品无人区| 免费女性裸体啪啪无遮挡网站| 人人妻人人澡人人爽人人夜夜| 成年人免费黄色播放视频| 另类亚洲欧美激情| 午夜福利乱码中文字幕| 日韩 欧美 亚洲 中文字幕| 天天影视国产精品| 丁香六月天网| 国产在线一区二区三区精| 国产精品蜜桃在线观看| 一区二区三区激情视频| 韩国av在线不卡| 777久久人妻少妇嫩草av网站| 老司机深夜福利视频在线观看 | 青春草视频在线免费观看| 国产一区有黄有色的免费视频| av在线app专区| 韩国精品一区二区三区| 韩国精品一区二区三区| 伊人久久大香线蕉亚洲五| 伊人久久大香线蕉亚洲五| www.自偷自拍.com| 乱人伦中国视频| 精品亚洲乱码少妇综合久久| 乱人伦中国视频| 亚洲成人免费av在线播放| 男女午夜视频在线观看| 十八禁高潮呻吟视频| 久久久国产欧美日韩av| 纯流量卡能插随身wifi吗| 黄色毛片三级朝国网站| 日本猛色少妇xxxxx猛交久久| 哪个播放器可以免费观看大片| 欧美人与性动交α欧美精品济南到| 如日韩欧美国产精品一区二区三区| 午夜日本视频在线| 十八禁人妻一区二区| 极品人妻少妇av视频| 国产97色在线日韩免费| 日韩制服骚丝袜av| 久久这里只有精品19| 国产一区二区在线观看av| 操美女的视频在线观看| 中文字幕色久视频| 午夜激情久久久久久久| 多毛熟女@视频| 制服人妻中文乱码| 综合色丁香网| 美女福利国产在线| 精品卡一卡二卡四卡免费| 9191精品国产免费久久| 亚洲综合精品二区| 亚洲国产精品成人久久小说| 国产 精品1| 久久人人爽人人片av| 热re99久久精品国产66热6| 一边摸一边做爽爽视频免费| 亚洲欧美清纯卡通| 在现免费观看毛片| 最近2019中文字幕mv第一页| 九色亚洲精品在线播放| 国产日韩欧美视频二区| 热re99久久精品国产66热6| 黄色 视频免费看| 国产亚洲精品第一综合不卡| 精品视频人人做人人爽| 亚洲精品一二三| 啦啦啦 在线观看视频| 亚洲成人av在线免费| 精品免费久久久久久久清纯 | 搡老乐熟女国产| 免费观看性生交大片5| av国产久精品久网站免费入址| 老司机深夜福利视频在线观看 | 黑丝袜美女国产一区| 制服丝袜香蕉在线| 在线看a的网站| 婷婷色麻豆天堂久久| 街头女战士在线观看网站| 哪个播放器可以免费观看大片| 亚洲美女视频黄频| 国产一区亚洲一区在线观看| 亚洲欧洲日产国产| 午夜免费鲁丝| 欧美乱码精品一区二区三区| 欧美 亚洲 国产 日韩一| 国产精品一区二区精品视频观看| 亚洲熟女毛片儿| 中文字幕亚洲精品专区| 丝袜喷水一区| 国产又爽黄色视频| 免费在线观看完整版高清| 久久国产精品大桥未久av| 男女下面插进去视频免费观看| 只有这里有精品99| 日韩伦理黄色片| 精品久久蜜臀av无| 色视频在线一区二区三区| 99九九在线精品视频| 国产精品久久久久久精品电影小说| 国产亚洲一区二区精品| 女的被弄到高潮叫床怎么办| 男女午夜视频在线观看| 80岁老熟妇乱子伦牲交| 久久久久精品人妻al黑| 中文字幕人妻丝袜制服| 久久天躁狠狠躁夜夜2o2o | 国产有黄有色有爽视频| 日韩一本色道免费dvd| 观看av在线不卡| 少妇的丰满在线观看| 黑人欧美特级aaaaaa片| 久久韩国三级中文字幕| 老司机影院成人| 老司机靠b影院| 精品亚洲成国产av| 久久久精品免费免费高清| 19禁男女啪啪无遮挡网站| 你懂的网址亚洲精品在线观看| 日韩大码丰满熟妇| 99热网站在线观看| 男女高潮啪啪啪动态图| 亚洲婷婷狠狠爱综合网| 国产黄色视频一区二区在线观看| 老司机靠b影院| 免费少妇av软件| 精品少妇黑人巨大在线播放| 99国产综合亚洲精品| 亚洲精品日韩在线中文字幕| 久久久精品国产亚洲av高清涩受| 国产精品三级大全| 国产成人免费无遮挡视频| 美女国产高潮福利片在线看| 国产乱来视频区| 亚洲成人免费av在线播放| a级片在线免费高清观看视频| 人人妻人人澡人人爽人人夜夜| 国产精品久久久久久人妻精品电影 | 不卡av一区二区三区| 欧美日韩av久久| www.自偷自拍.com| 国产精品香港三级国产av潘金莲 | 亚洲视频免费观看视频| 久久久国产精品麻豆| 电影成人av| 最黄视频免费看| 在线观看三级黄色| 免费观看性生交大片5| 9热在线视频观看99| 三上悠亚av全集在线观看| 国产淫语在线视频| 国产成人精品在线电影| a级毛片在线看网站| 七月丁香在线播放| 色视频在线一区二区三区| 久久午夜综合久久蜜桃| 女人爽到高潮嗷嗷叫在线视频| 麻豆乱淫一区二区| 天天影视国产精品| 久久精品人人爽人人爽视色| 十八禁人妻一区二区| 狂野欧美激情性xxxx| 久久97久久精品| 欧美日韩av久久| 久久精品国产a三级三级三级| 国产精品国产av在线观看| 在线天堂最新版资源| 国产 一区精品| 国产亚洲欧美精品永久| av卡一久久| 另类精品久久| h视频一区二区三区| 免费高清在线观看日韩| 国产精品成人在线| 精品久久蜜臀av无| 99久久综合免费| 一边摸一边做爽爽视频免费| 国产一区有黄有色的免费视频| 亚洲国产日韩一区二区| 1024视频免费在线观看| 一本—道久久a久久精品蜜桃钙片| 色视频在线一区二区三区| 成年人午夜在线观看视频| 亚洲av在线观看美女高潮| 久久av网站| 青春草国产在线视频| 国产免费又黄又爽又色| 午夜福利在线免费观看网站| 制服诱惑二区| www.熟女人妻精品国产| 亚洲四区av| 精品亚洲成国产av| 大香蕉久久网| 在线观看免费日韩欧美大片| 免费av中文字幕在线| 美国免费a级毛片| 狠狠精品人妻久久久久久综合| 免费在线观看完整版高清| 美女大奶头黄色视频| 国产精品人妻久久久影院| 9191精品国产免费久久| 晚上一个人看的免费电影| 亚洲精华国产精华液的使用体验| 亚洲第一区二区三区不卡| 日本一区二区免费在线视频| 久久久久久久久免费视频了| 国产精品免费视频内射| 美女午夜性视频免费| 狠狠婷婷综合久久久久久88av| bbb黄色大片| 久久97久久精品| 啦啦啦啦在线视频资源| 国产又色又爽无遮挡免| 交换朋友夫妻互换小说| 亚洲成国产人片在线观看| 精品国产国语对白av| 狠狠精品人妻久久久久久综合| 热99久久久久精品小说推荐| 亚洲国产欧美一区二区综合| 狂野欧美激情性xxxx| 国产在视频线精品| 精品亚洲成国产av| 90打野战视频偷拍视频| 精品国产乱码久久久久久小说| 黄色怎么调成土黄色| 在线观看www视频免费| 亚洲欧美成人精品一区二区| a级毛片在线看网站| 高清av免费在线| 国产老妇伦熟女老妇高清| 少妇精品久久久久久久| 90打野战视频偷拍视频| 在线观看免费午夜福利视频| 国产成人av激情在线播放| 99热全是精品| 久久久精品免费免费高清| 电影成人av| 激情五月婷婷亚洲| 777久久人妻少妇嫩草av网站| 一区福利在线观看| 狠狠婷婷综合久久久久久88av| 国产成人91sexporn| 热re99久久国产66热| 亚洲美女搞黄在线观看| 国产精品.久久久| 午夜老司机福利片| 天堂8中文在线网| 免费看av在线观看网站| 亚洲久久久国产精品| 一区在线观看完整版| 亚洲色图综合在线观看| 国产精品久久久久久久久免| 大香蕉久久网| 性少妇av在线| 老司机影院成人| 国产又爽黄色视频| 久久人妻熟女aⅴ| 日本猛色少妇xxxxx猛交久久| 国产在视频线精品| 久久久久精品性色| 美女中出高潮动态图| 欧美老熟妇乱子伦牲交| 国产成人精品无人区| 一边亲一边摸免费视频| 狂野欧美激情性bbbbbb| 蜜桃在线观看..| 麻豆av在线久日| 久久综合国产亚洲精品| 国产精品久久久久久精品古装| 捣出白浆h1v1| 久久精品久久久久久噜噜老黄| 美女福利国产在线| 高清在线视频一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 久久久久视频综合| 在现免费观看毛片| 美女中出高潮动态图| 亚洲国产中文字幕在线视频| 女人爽到高潮嗷嗷叫在线视频| 无遮挡黄片免费观看| 国产成人精品无人区| 久久影院123| 婷婷成人精品国产| 国产一区二区激情短视频 | 午夜福利视频精品| 十八禁高潮呻吟视频| 男女免费视频国产| svipshipincom国产片| 青春草国产在线视频| 黄色视频不卡| 日韩一卡2卡3卡4卡2021年| 午夜激情久久久久久久| 国产福利在线免费观看视频| 午夜福利在线免费观看网站| 中国三级夫妇交换| 黄色视频在线播放观看不卡| 国产成人午夜福利电影在线观看| 久久精品亚洲av国产电影网| 丰满迷人的少妇在线观看| 黄色 视频免费看| 欧美最新免费一区二区三区| 国产精品人妻久久久影院| 亚洲欧美精品综合一区二区三区| 国产伦理片在线播放av一区| 精品国产乱码久久久久久小说| 王馨瑶露胸无遮挡在线观看| 国产日韩欧美在线精品| 亚洲精品自拍成人| 日本爱情动作片www.在线观看| 乱人伦中国视频| 精品福利永久在线观看| 国产人伦9x9x在线观看| 女的被弄到高潮叫床怎么办| 丝袜美腿诱惑在线| √禁漫天堂资源中文www| 国产无遮挡羞羞视频在线观看| 看十八女毛片水多多多| 国产精品人妻久久久影院| 狂野欧美激情性bbbbbb| 亚洲精品美女久久久久99蜜臀 | 亚洲av国产av综合av卡| xxxhd国产人妻xxx| 无遮挡黄片免费观看| 日韩av在线免费看完整版不卡| 天美传媒精品一区二区| 新久久久久国产一级毛片| 精品福利永久在线观看| 国产成人av激情在线播放| 涩涩av久久男人的天堂| 久久免费观看电影| 久久青草综合色| 欧美av亚洲av综合av国产av | 国产在视频线精品| 赤兔流量卡办理| 国产 一区精品| 高清欧美精品videossex| 精品一品国产午夜福利视频| 两性夫妻黄色片| av电影中文网址| 国产一卡二卡三卡精品 | 亚洲欧洲精品一区二区精品久久久 | 天天躁日日躁夜夜躁夜夜| 国产色婷婷99| 亚洲国产精品一区二区三区在线| 在线观看国产h片| 天天躁夜夜躁狠狠躁躁| 欧美精品一区二区大全| xxxhd国产人妻xxx| 久久久国产欧美日韩av| 久久久亚洲精品成人影院| 国精品久久久久久国模美| 一个人免费看片子| av在线观看视频网站免费| 精品少妇内射三级| 成人午夜精彩视频在线观看| 国产极品粉嫩免费观看在线| 老司机在亚洲福利影院| 午夜福利,免费看| 久久久久精品久久久久真实原创| av电影中文网址| 久久狼人影院| 看非洲黑人一级黄片| 五月天丁香电影| 国产欧美亚洲国产| 色综合欧美亚洲国产小说| 欧美成人精品欧美一级黄| 最近最新中文字幕大全免费视频 | 国产99久久九九免费精品| 欧美在线黄色| 18在线观看网站| 女人爽到高潮嗷嗷叫在线视频| 另类亚洲欧美激情| 午夜老司机福利片| 亚洲精品成人av观看孕妇| 波野结衣二区三区在线| 国产片内射在线| 精品人妻一区二区三区麻豆| 欧美日韩成人在线一区二区| 精品久久久精品久久久| 在线精品无人区一区二区三| 精品久久蜜臀av无| 亚洲精品国产av蜜桃| 男女边摸边吃奶| 欧美日韩国产mv在线观看视频| a级片在线免费高清观看视频| 中国三级夫妇交换| 晚上一个人看的免费电影| 交换朋友夫妻互换小说| 91精品伊人久久大香线蕉| 熟女av电影| 在线观看三级黄色| 女人被躁到高潮嗷嗷叫费观| 电影成人av| 日韩免费高清中文字幕av| 亚洲国产日韩一区二区| 黑人猛操日本美女一级片| 毛片一级片免费看久久久久| 久久天堂一区二区三区四区| 少妇的丰满在线观看| 国产欧美日韩一区二区三区在线| 久久久亚洲精品成人影院| 女人爽到高潮嗷嗷叫在线视频| 国产在线视频一区二区| 精品一区在线观看国产| 日本91视频免费播放| av免费观看日本| 老司机在亚洲福利影院| 国产精品av久久久久免费| 啦啦啦中文免费视频观看日本| 午夜福利乱码中文字幕| 亚洲一级一片aⅴ在线观看| 亚洲国产成人一精品久久久| 欧美精品人与动牲交sv欧美| 亚洲欧美一区二区三区久久| 不卡视频在线观看欧美| 丰满乱子伦码专区| 巨乳人妻的诱惑在线观看| av福利片在线| 国产亚洲av高清不卡| 久久99一区二区三区| 一区二区三区乱码不卡18| 伦理电影大哥的女人| 91老司机精品| 国产av一区二区精品久久| 国产一区亚洲一区在线观看| 熟妇人妻不卡中文字幕| 国产精品无大码| 亚洲国产精品成人久久小说| av.在线天堂| 啦啦啦在线免费观看视频4| 国产成人午夜福利电影在线观看| 久久久精品国产亚洲av高清涩受| 美女视频免费永久观看网站| 久久久久久久久久久久大奶| 亚洲四区av| 丝袜美腿诱惑在线| 国产免费视频播放在线视频| 精品少妇黑人巨大在线播放| 男女高潮啪啪啪动态图| 波多野结衣av一区二区av| 国产国语露脸激情在线看| 久久久久久人人人人人| 欧美日韩视频精品一区|