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

    藏北高原礫石粒徑空間異質(zhì)性研究

    2023-03-16 01:55:02胡孟珂凌鵬飛
    干旱區(qū)研究 2023年2期
    關(guān)鍵詞:分異礫石土地利用

    徐 濤, 于 歡, 孔 博, 邱 霞,3, 胡孟珂, 凌鵬飛

    (1.成都理工大學(xué)地球科學(xué)學(xué)院,四川 成都 610059;2.中國(guó)科學(xué)院水利部成都山地與災(zāi)害研究所,四川 成都610041;3.四川省不動(dòng)產(chǎn)登記中心,四川 成都 610014)

    平均海拔4500 m以上的藏北高原,是我國(guó)的重要生態(tài)安全屏障[1]。但是,近年來(lái)由于氣候變暖以及人類活動(dòng)的加劇,藏北高原地區(qū)高寒草地退化嚴(yán)重,植被覆蓋度下降,生態(tài)環(huán)境被嚴(yán)重破壞[2],這直接影響著藏北地區(qū)農(nóng)牧民的生產(chǎn)生活,也對(duì)青藏高原乃至全國(guó)生態(tài)造成嚴(yán)重威脅[3]。因此,藏北高原生態(tài)環(huán)境恢復(fù)成為我國(guó)生態(tài)工作的重要內(nèi)容[4]。而在藏北高原生態(tài)領(lǐng)域研究中,多數(shù)學(xué)者以草原植被和土壤組分為主要研究對(duì)象[5-10],以高原礫石作為對(duì)象的研究相對(duì)較少。一方面,礫石是各種水文和侵蝕等過(guò)程綜合作用的產(chǎn)物,是草地和土壤退化,生態(tài)系統(tǒng)惡化的一個(gè)標(biāo)志;另一方面礫石反過(guò)來(lái)又影響水文和侵蝕過(guò)程,如入滲、蒸發(fā)、徑流和水蝕等[11]。因此,本文以地表礫石粒徑大小和空間位置為研究對(duì)象,深入分析礫石分異空間格局的效應(yīng)關(guān)系,以期為藏北高原生態(tài)修復(fù)提供參考。

    在對(duì)礫石粒徑空間分異的研究中,一方面在研究對(duì)象上多是針對(duì)沙丘、瀕海(河)的沙粒以及戈壁地區(qū)礫石的空間分異,如Okin等[12]在加利福利亞?wèn)|南部的馬尼克斯(Manix)研究了羽狀沙丘沙??臻g分異現(xiàn)象;董玉祥等[13]對(duì)河北黃金海岸沙粒的垂直分布模式進(jìn)行研究;陳西慶等[14]通過(guò)對(duì)長(zhǎng)江入河口底沙粒徑的變化進(jìn)行研究,結(jié)果表明研究區(qū)河口底沙粒徑存在繼續(xù)粗化的趨勢(shì)。另一方面在研究方法上多數(shù)利用經(jīng)典統(tǒng)計(jì)學(xué)分析海拔、粒度參數(shù)與粒徑的相關(guān)關(guān)系,或利用地統(tǒng)計(jì)學(xué)中的空間變異函數(shù)模型分析礫石粒徑的空間異質(zhì)性,如Quade[15]通過(guò)研究位于莫哈維戈壁表面的礫石分布隨海拔的變化,構(gòu)建了礫石覆蓋度隨海拔變化的函數(shù);曹曉陽(yáng)等[16]利用經(jīng)典統(tǒng)計(jì)學(xué)方法分析了噶順戈壁洪積扇區(qū)地表礫石的粒度特征;王利兵等[17]采用地統(tǒng)計(jì)學(xué)的理論和方法研究了渾善達(dá)克沙地沙粒粒徑組成的空間異質(zhì)性。

    上述研究從統(tǒng)計(jì)學(xué)及空間分析等方面對(duì)礫石粒徑空間分異展開相應(yīng)研究,但是,對(duì)高原礫石粒徑空間分異以及不同因子疊加貢獻(xiàn)率的研究報(bào)道較少。地理探測(cè)器模型是一種基于空間異質(zhì)性的研究方法,可以定量探測(cè)表達(dá)某一時(shí)空現(xiàn)象的主要驅(qū)動(dòng)因子以及不同驅(qū)動(dòng)因子之間交互作用,與其他空間異質(zhì)性探測(cè)工具相比地理探測(cè)器具有更高的解釋效率[18]。因此,本文在借鑒前人研究方法的基礎(chǔ)上,加入地理探測(cè)器模型開展藏北高原礫石粒徑空間分異的多因子定量歸因,并結(jié)合回歸分析建立礫石粒徑的空間預(yù)測(cè)模型,以期為藏北高原生態(tài)修復(fù)及環(huán)境保護(hù)提供科學(xué)依據(jù)。

    1 研究區(qū)與研究方法

    1.1 研究區(qū)概況

    由于藏北高原面積較大,且部分地區(qū)礫石分布不具代表性,為系統(tǒng)分析藏北高原礫石的空間異質(zhì)性,故研究區(qū)選擇礫石分布特征顯著的那曲地區(qū)南部地區(qū)(29°55′21″~33°53′22″N,87°22′21″~97°3′12″E),包括安多縣南部(除色務(wù)鄉(xiāng)外其余安多縣境內(nèi)區(qū)域)、班戈縣、申扎縣、雙湖縣南部(含多瑪鄉(xiāng)、措折羅瑪鎮(zhèn)、巴嶺鄉(xiāng)、協(xié)德鄉(xiāng))以及色尼區(qū)。研究區(qū)總面積約18.31 km2,平均海拔4100 m以上,地勢(shì)四周高、中間低,氣候寒冷干燥,大部分地區(qū)年均氣溫低于0 ℃。年均降水量從東部的548.8 mm 左右,逐漸減至西部的216.3 mm,且降水主要集中在6—9月。研究區(qū)土壤以寒鈣土和草氈土為主,植被類型較為單一,90%左右均為草甸。受到自然地理?xiàng)l件的影響,研究區(qū)經(jīng)濟(jì)發(fā)展水平較差,經(jīng)濟(jì)發(fā)展以放牧業(yè)與旅游業(yè)為主,人口密度較低,僅在東部地區(qū)有少量人口分布,地廣人稀。

    1.2 數(shù)據(jù)來(lái)源

    礫石粒徑空間分布數(shù)據(jù)來(lái)自野外采樣與室內(nèi)反演。野外采樣點(diǎn)位于藏北高原的那曲南部地區(qū)(圖1),海拔每升高100 m 設(shè)一條樣帶,每個(gè)樣帶設(shè)4~5個(gè)樣區(qū),每個(gè)樣區(qū)的4個(gè)邊角及中間位置設(shè)置1個(gè)0.5 m×0.5 m 的樣方(圖1d),由于區(qū)域情況復(fù)雜,部分地區(qū)受現(xiàn)實(shí)條件的限制,無(wú)法測(cè)得5個(gè)樣方,總計(jì)測(cè)得30個(gè)樣區(qū)(圖1b),145個(gè)樣方。每個(gè)樣方內(nèi)選取占主體粒級(jí)10個(gè)以內(nèi)的礫石,并用游標(biāo)卡尺對(duì)單個(gè)礫石的Feret 直徑進(jìn)行測(cè)量(圖1e),并記錄海拔、坡度、坡向、礫石數(shù)量等基本信息。室內(nèi)反演主要將145個(gè)樣方分為130個(gè)反演樣方和15個(gè)驗(yàn)證樣方,通過(guò)采用ImageJ 軟件計(jì)算130 個(gè)反演樣方中經(jīng)形貌輪廓提取后的野外采樣礫石的粒徑值,并對(duì)研究區(qū)Landsat 8 原始影像、光譜指數(shù)因子、DEM 數(shù)據(jù)等進(jìn)行了預(yù)處理并去除因子間多重共線性,分析粒徑值與上述因子間的相關(guān)性,篩選顯著變量因子,結(jié)合多元線性逐步回歸與基于PCA 的多元回歸分析構(gòu)建礫石特征參數(shù)反演模型,進(jìn)而得到礫石粒徑反演數(shù)據(jù)[19]。該數(shù)據(jù)經(jīng)15個(gè)驗(yàn)證樣方驗(yàn)證,均方根誤差為1.4,相對(duì)誤差在30%以內(nèi)的樣區(qū)數(shù)為11個(gè),約占總樣區(qū)數(shù)的73%,平均誤差為22.75%,反演精度良好。

    圖1 野外采樣與反演成果Fig.1 Field sampling and inversion results

    根據(jù)Wentworth[20]提出的礫石分級(jí)理論,將本次采集的礫石粒徑分為3 種等級(jí)。其中粒徑在2~4 mm之間的礫石分為極細(xì)礫,粒徑在4~8 mm之間的分為細(xì)礫,粒徑在8~16 mm 的礫石分為中礫。本研究反演得到的礫石粒徑分布如圖1c 所示,由圖可知,研究區(qū)粒徑主要為4~8 mm的細(xì)礫,分布于研究區(qū)的中西部,占總面積的76.99%,其次為8~16 mm的中礫,主要分布于研究區(qū)東部,占總面積的19.32%,而極細(xì)粒僅占總面積的3.69%。

    環(huán)境變量數(shù)據(jù)有自然因素、人文因素兩大類,共8個(gè)代理變量,數(shù)據(jù)信息如表1所示。

    表1 數(shù)據(jù)來(lái)源及介紹Tab.1 Data source and introduction

    1.3 研究方法

    1.3.1 數(shù)據(jù)預(yù)處理方法 由于地理探測(cè)器的正確運(yùn)行需要分類數(shù)據(jù),故對(duì)土地利用類型、土壤類型、植被類型采用原始分類,其余影響因子采用自然間斷點(diǎn)分類法進(jìn)行離散化,經(jīng)分異及因子探測(cè)器不斷調(diào)試,最終確定NDVI、年均降水分9類,夜間燈光強(qiáng)度分8類,高程分7類,人口密度分4類。

    利用ArcGIS 漁網(wǎng)工具將研究區(qū)劃分為10 km×10 km 的研究單元,然后采用雙線性插值法提取每個(gè)研究單元內(nèi)礫石粒徑值Y與分類后環(huán)境變量的類型量X,從而將兩種數(shù)據(jù)的屬性進(jìn)行匹配,最后將提取的數(shù)據(jù)帶入地理探測(cè)器模型進(jìn)行計(jì)算[21]。具體流程如圖2所示。

    圖2 數(shù)據(jù)處理流程Fig.2 Data processing flow

    1.3.2 空間異質(zhì)性方法

    (1)空間局域異質(zhì)性Moran’sI指數(shù)可以分析空間局域異質(zhì)性,即分析某區(qū)域內(nèi)礫石粒徑與周圍礫石粒徑的空間自相關(guān)程度。Moran’sI指數(shù)介于(-1,1),I值越接近于1,表示研究單元正相關(guān)性越強(qiáng);I值越接近于-1,表示研究單元負(fù)相關(guān)性越強(qiáng);I值趨近于0時(shí),表示研究單元呈隨機(jī)分布[22]。

    式中:I為Moran’sI指數(shù);n為研究礫石單元個(gè)數(shù);Yi和Yj分別為第i個(gè)礫石單元和第j個(gè)礫石單元的粒徑值;Wij為空間權(quán)重矩陣;S2為觀測(cè)值的方差;Yˉ為粒徑值的平均值[23]。

    (2)空間結(jié)構(gòu)異質(zhì)性 為進(jìn)一步分析空間自相關(guān)的范圍及空間結(jié)構(gòu)異質(zhì)性,引入空間變異函數(shù)探究其空間結(jié)構(gòu)異質(zhì)性,并計(jì)算各向同性及各向異性下礫石粒徑的最佳擬合模型及特征參數(shù)值,其公式為:

    式中:z(xi)和z(xi+h)分別為區(qū)域化隨機(jī)變量z在空間中xi和xi+h處的值;N(h)為抽樣間隔等于h時(shí)的點(diǎn)對(duì)數(shù);R(h)為變異函數(shù)值[24]。

    (3)空間分層異質(zhì)性 礫石的空間異質(zhì)性極易受環(huán)境條件的影響,為此本文引入地理探測(cè)器探測(cè)影響藏北高原礫石粒徑空間分異的主導(dǎo)因子。地理探測(cè)器可以將變量進(jìn)行分層,用來(lái)探究空間分層異質(zhì)性,并揭示變量間存在某種關(guān)聯(lián)的空間分析模型[25]。該模型被廣泛應(yīng)用于土地利用、區(qū)域經(jīng)濟(jì)、精準(zhǔn)扶貧等領(lǐng)域,是一種新型空間分異性分析工具[26]。在本研究中用以探測(cè)礫石粒徑的空間分異以及環(huán)境變量因子對(duì)礫石粒徑的解釋程度,用q值度量,表達(dá)式如下:

    式中:N為全區(qū)的單元數(shù);σ2為指標(biāo)的方差;h為變量的分區(qū),h=1,2,3,…,L;L表示分區(qū)數(shù)目。q的大小反映了空間分異的程度,q值越大,表示空間分層異質(zhì)性越強(qiáng),反之則空間分布的隨機(jī)性越強(qiáng)。

    2 結(jié)果與分析

    2.1 基于Moran’s I指數(shù)的空間局域異質(zhì)性

    Moran’sI指數(shù)可以用來(lái)探究研究區(qū)礫石粒徑的局域異質(zhì)性,研究區(qū)全局Moran’sI指數(shù)為0.481,Z值高于2.5 且通過(guò)1%的置信檢驗(yàn),表明研究區(qū)礫石粒徑呈顯著的正相關(guān),空間聚集性較強(qiáng),存在異質(zhì)性。局部Moran’sI指數(shù)結(jié)果顯示(圖3),研究區(qū)礫石粒徑聚集模式主要為高-高聚集模式、低-低聚集模式與隨機(jī)模式。在空間分布上,研究區(qū)東部礫石多呈現(xiàn)高-高聚集模式,研究區(qū)南部礫石則以隨機(jī)分布為主,在研究區(qū)西部和北部地區(qū),隨機(jī)分布和低-低聚集呈現(xiàn)相間分布。

    圖3 研究區(qū)礫石粒徑的空間聚集特征Fig.3 Spatial clustering characteristics of gravel size in the study area

    2.2 基于空間變異函數(shù)的空間結(jié)構(gòu)異質(zhì)性

    通過(guò)空間變異函數(shù)計(jì)算出研究區(qū)礫石粒徑的一系列變異參數(shù),并選擇最佳空間變異函數(shù)模型進(jìn)行擬合,從而得到研究區(qū)礫石粒徑的空間變異曲線(圖4a)??紤]到礫石粒徑在不同方向上的變化存在一定差異性,本文從各向異性的角度出發(fā),將研究區(qū)礫石粒徑劃分為0°、22.5°、45°、67.5°、90°、112.5°、135°及157.5°共8 個(gè)方向,各個(gè)方向上的點(diǎn)線圖如圖4b所示。

    圖4 礫石粒徑空間變異曲線Fig.4 Spatial variation curve of gravel size

    由圖4可知,從各向同性的角度出發(fā),總的礫石粒徑空間變異函數(shù)類型為線性,而從各向異性的角度來(lái)看,0°、22.5°、45°、157.5°方向上變異函數(shù)類型為指數(shù)模型,90°、112.5°、135°方向上變異函數(shù)類型為高斯模型,67.5°方向上為線性模型。

    根據(jù)變異函數(shù)曲線,得到研究區(qū)礫石粒徑空間結(jié)構(gòu)異質(zhì)性相關(guān)特征參數(shù)如表2所示。

    表2 礫石粒徑空間變異函數(shù)特征參數(shù)Tab.2 Characteristic parameters of spatial variogram of gravel size

    從各向同性角度分析,塊金值C0為0.391,具有一定的塊金效應(yīng),亦說(shuō)明隨機(jī)因素引起的空間異質(zhì)性不可忽略。偏基臺(tái)值代表結(jié)構(gòu)因素引起的空間變異,為0.708,大于塊金值,表明研究區(qū)礫石粒徑的空間異質(zhì)性主要由結(jié)構(gòu)因素引起。基臺(tái)值為1.099,代表了由隨機(jī)因素及結(jié)構(gòu)因素共同引起的空間變異,亦表明了研究區(qū)礫石粒徑的空間異質(zhì)性不可忽略。塊基比為0.356,表明研究區(qū)礫石粒徑空間異質(zhì)性的因素中,隨機(jī)因素占比0.356,結(jié)構(gòu)因素占比0.644,結(jié)構(gòu)因素占主導(dǎo)地位。變程為26.178 km,表明研究區(qū)礫石粒徑的空間自相關(guān)范圍為26.178 km。

    從各向異性的角度分析,各方向特征參數(shù)值存在差異性,但在任意方向上,偏基臺(tái)值均大于塊金值,說(shuō)明研究區(qū)礫石粒徑空間異質(zhì)性均由結(jié)構(gòu)因素主導(dǎo)??臻g變異函數(shù)變程的大小反映礫石粒徑沿某一方向變化速度的快慢,變程值越大,表明該參數(shù)在該方向上變化慢,即非勻質(zhì)性越弱,反之越強(qiáng)。由表可知,157.5°方向變程值最大,即該方向變異性最小,而90°方向變程值最小,即該方向變異性最大。故研究區(qū)礫石粒徑各向異性特征非常明顯,各向異性主軸為157.5°方向。

    2.3 基于地理探測(cè)器的空間分層異質(zhì)性

    因子探測(cè)器探測(cè)各因子對(duì)礫石粒徑的解釋度(q值)大小。探測(cè)結(jié)果如表3所示,通過(guò)顯著性檢驗(yàn)的因子中,按對(duì)礫石粒徑解釋度強(qiáng)弱排序?yàn)椋篨2(NDVI)>X3(土地利用類型)>X7(人口密度)>X5(植被類型)>X6(年均降水),對(duì)應(yīng)的q值分別為0.41、0.27、0.16、0.15、0.13。其中,NDVI 的因子解釋力最大(0.41),表明藏北高原礫石粒徑空間分布受NDVI的控制作用最為強(qiáng)烈,即NDVI 與礫石粒徑在空間分布上具有最強(qiáng)的一致性。次要因素為土地利用類型,因子解釋力為0.27,說(shuō)明土地利用類型也是影響礫石粒徑空間分布的重要因素。

    表3 礫石粒徑因子探測(cè)結(jié)果Tab.3 Gravel size factor detection results

    交互探測(cè)主要分析各因子對(duì)礫石粒徑空間分布是否存在交互作用。結(jié)果如表4所示,由表可知,除X3∩X6表現(xiàn)為非線性增強(qiáng)外,其余交互作用類型均表現(xiàn)為雙因子增強(qiáng)。其中X2(NDVI)∩X3(土地利用類型)交互作用值最強(qiáng),解釋力達(dá)到0.56,其次為X2(NDVI)∩X7(人口密度)與X3(土地利用類型)∩X6(年均降水),解釋力均為0.47,最小的為X6(年均降水)∩X7(人口密度),解釋力為0.27。各個(gè)因子間的交互作用共同促成了研究區(qū)礫石粒徑的空間分異格局。這說(shuō)明任意兩個(gè)環(huán)境因子交互后對(duì)礫石粒徑空間分布的因子解釋力均會(huì)顯著提升,也就是說(shuō)礫石粒徑空間分布是由各維度環(huán)境因子的共同作用形成。

    表4 礫石粒徑影響因子間交互作用Tab.4 Interaction between impact factors of gravel size

    風(fēng)險(xiǎn)探測(cè)器探測(cè)各因子不同分層對(duì)礫石粒徑空間分布是否存在顯著性差異及各分層中礫石粒徑均值大小。探測(cè)結(jié)果如圖5所示,NDVI 中,0.79~0.92 分區(qū)中粒徑均值最大,0.1~0.18 分區(qū)粒徑均值最小,總體粒徑均值與NDVI呈正相關(guān),除0.24~0.32與0.61~0.79 與0.79~0.92 分層外,其余區(qū)間均存在顯著性差異。在土地利用類型中,疏林地中的粒徑均值最大,鹽堿地中粒徑均值最小,絕大多數(shù)分層間均存在顯著性差異。植被類型中高山墊狀植被中的粒徑均值最大,除高山墊狀植被分層外,其余均不存在顯著性差異。年降水量中,年降水量在470.4~548.8 mm 區(qū)間內(nèi)粒徑均值最大,261.3~281.4 mm區(qū)間內(nèi)粒徑均值最小,總體兩者呈正相關(guān)。

    圖5 影響因子分區(qū)粒徑大小Fig.5 Influence factor partition particle size

    人口密度因子中,人口密度與礫石粒徑兩者呈現(xiàn)一定的正相關(guān)。但是,NDVI 與人口密度均呈現(xiàn)“東高西低”的分布態(tài)勢(shì),兩者具有空間分布上的相似性,且該區(qū)人口稀少,人類活動(dòng)對(duì)礫石粒徑的影響相對(duì)有限,結(jié)合因子探測(cè)器探測(cè)結(jié)果,NDVI 對(duì)礫石粒徑的解釋力遠(yuǎn)大于人口密度,NDVI 仍為影響礫石粒徑分布的主導(dǎo)因素。NDVI 一定程度上代表了植被覆蓋度,NDVI 大,則植被覆蓋度相對(duì)較高,水土保持能力強(qiáng),流水、風(fēng)力等外力侵蝕作用較弱,則礫石粒徑值較大;反之,礫石粒徑越大,則土壤水分的入滲能力越強(qiáng),則植被覆蓋度越高。因此,人口分布與礫石粒徑可能存在某種假相關(guān),即兩者空間分布的相似性可能是由于植被等外部環(huán)境引起的。

    2.4 多元線性回歸、隨機(jī)森林回歸與最優(yōu)尺度回歸驗(yàn)證

    由地理探測(cè)器的因子探測(cè)器及交互作用探測(cè)器可知,上述8 種影響因子均對(duì)礫石粒徑有不同程度的影響,但是DEM、土壤類型、夜間燈光強(qiáng)度P值過(guò)大,未通過(guò)0.05 的顯著性檢驗(yàn),故將NDVI、土地利用、植被類型、年均降水、人口密度利用多元線性回歸進(jìn)行擬合,擬合結(jié)果顯示各因子的VIF值均遠(yuǎn)<5,說(shuō)明各因子間不存在多重共線性,模型擬合調(diào)整后如圖6a所示,R2為0.426。將上述8個(gè)影響因子加入隨機(jī)森林回歸模型進(jìn)行擬合,結(jié)果如圖6b 所示,由圖可知,隨機(jī)森林回歸模型擬合R2為0.508,隨機(jī)森林的R2高于多元線性回歸。

    圖6 多元線性回歸與隨機(jī)森林回歸擬合Fig.6 Multiple linear regression and random forest regression fitting

    考慮到影響因素中有土地利用類型、植被類型等名義變量,利用多元線性回歸與隨機(jī)森林回歸模型進(jìn)行擬合,無(wú)法準(zhǔn)確地反映分類變量不同取值的距離,故本文引入最優(yōu)尺度回歸進(jìn)行擬合。將礫石粒徑、NDVI、年均降水、人口密度作為數(shù)字型變量,土地利用、植被類型作為名義變量輸入模型,結(jié)果顯示最優(yōu)尺度回歸調(diào)整后R2為0.491,擬合結(jié)果好于多元線性回歸,但次于隨機(jī)森林回歸,且人口密度未通過(guò)0.05的顯著性檢驗(yàn),原因可能為人口密度與礫石粒徑存在某種假相關(guān),考慮剔除人口密度再次進(jìn)行回歸,回歸結(jié)果如表5所示。

    表5 最優(yōu)尺度回歸模型Tab.5 Optimal scale regression model

    改進(jìn)后模型擬合精度由0.491 上升到0.518,通過(guò)0.01的顯著性檢驗(yàn),擬合結(jié)果好于多元線性回歸與隨機(jī)森林回歸,且各因子均通過(guò)0.05的顯著性檢驗(yàn),但該模型僅用于分析研究區(qū)的礫石粒徑,模型為:

    式中:Y為礫石粒徑的標(biāo)準(zhǔn)分;X2、X3、X5、X6分別為NDVI、土地利用類型、植被類型、年均降水的標(biāo)準(zhǔn)分。由地理探測(cè)器篩選出的因子除人口密度外均有較顯著的水平,通過(guò)公式(4)系數(shù)絕對(duì)值大小可知,在研究區(qū)影響礫石粒徑大小的因子中,NDVI、土地利用類型起主要作用,植被類型、年均降水起次要作用,結(jié)果與地理探測(cè)器吻合。

    3 討論

    由于自然環(huán)境與人類活動(dòng)的影響,礫石粒徑存在空間異質(zhì)性。本研究發(fā)現(xiàn),藏北高原礫石粒徑存在一定程度的空間分異,且這種空間分異是由結(jié)構(gòu)因素主導(dǎo)。NDVI、土地利用類型是影響研究區(qū)礫石粒徑空間分異的主要因素,年均降水、植被類型、人口密度是影響研究區(qū)礫石粒徑空間分異的次要因素。

    地理探測(cè)器探測(cè)結(jié)果顯示,NDVI 與土地利用類型是影響藏北高原礫石粒徑空間分異的主要因素,這主要是由當(dāng)?shù)氐纳鷳B(tài)環(huán)境決定的(圖7)。在研究區(qū)中,降水、NDVI 呈現(xiàn)“東高西低”的態(tài)勢(shì),植被類型、土地利用類型以高寒禾草、苔草草原,高寒嵩草、雜草類草甸為主,四者具有空間分布上的一致性。降水影響植被覆蓋度和生長(zhǎng)狀況,植被類型反映了植被的背景信息,二者均對(duì)NDVI 與土地利用類型產(chǎn)生影響。NDVI 反映了植被覆蓋度與植被生長(zhǎng)狀態(tài),其與土地利用類型又與水流、風(fēng)力等外力侵蝕息息相關(guān)(如高覆蓋度草地與低覆蓋度草地相比,水流、風(fēng)力等外力侵蝕作用較弱)。而礫石一方面是各種水文和侵蝕等過(guò)程綜合作用的產(chǎn)物,另一方面礫石又通過(guò)入滲、改變土壤儲(chǔ)水量等反過(guò)來(lái)影響植被和土壤肥力,進(jìn)而影響NDVI 與土地利用類型。

    圖7 環(huán)境變量對(duì)礫石粒徑的影響機(jī)制Fig.7 The influence mechanism of environmental variables on gravel size

    綜上所述,可以從礫石粒徑角度分析草地退化趨勢(shì),為草原牧草退化提供動(dòng)態(tài)監(jiān)測(cè)。人類活動(dòng)導(dǎo)致藏北高原天然草地植被急劇下降,中高覆蓋草地向低覆蓋草地轉(zhuǎn)變,居民和工業(yè)用地面積增加,土地利用類型發(fā)生變化,導(dǎo)致流域水土保持能力下降,土壤破壞更為嚴(yán)重。植被退化與土地荒漠化的相互作用,惡化了區(qū)域生態(tài)環(huán)境,嚴(yán)重浪費(fèi)了區(qū)域水資源。由于供水不足和缺水,灌溉草地的牲畜負(fù)重過(guò)重,草地退化、荒漠化和水土流失更加嚴(yán)重,形成了惡性循環(huán)。上述變化將會(huì)使水土保持能力下降,外力侵蝕日益嚴(yán)重,礫石數(shù)量增多,粒徑變小。但受研究資料限制,缺乏長(zhǎng)時(shí)間序列的礫石粒徑分布數(shù)據(jù),因此無(wú)法確定礫石粒徑變化狀態(tài),后期將在時(shí)間尺度上進(jìn)一步研究礫石粒徑的空間分異,以期為藏北高原生態(tài)環(huán)境修復(fù)提供科學(xué)的參考依據(jù)。

    4 結(jié)論

    (1)研究區(qū)礫石粒徑呈顯著的正相關(guān),空間聚集性較強(qiáng),存在異質(zhì)性;在研究區(qū)東部礫石多呈現(xiàn)高-高聚集,研究區(qū)南部礫石則為隨機(jī)分布,在研究區(qū)西部和北部地區(qū),隨機(jī)分布和低-低聚集呈現(xiàn)相間分布。

    (2)研究區(qū)礫石粒徑具有明顯的各向異性特征,且由結(jié)構(gòu)因素主導(dǎo)。從變異函數(shù)模型來(lái)看,總的變異函數(shù)模型為線性,其余方向的模型為指數(shù)、高斯、線性;從變程來(lái)看,157.5°方向變異性最小,90°方向變異性最大。

    (3)因子探測(cè)器顯示NDVI、土地利用類型為主要影響因素,植被類型、年均降水與人口密度為次要影響因素;交互探測(cè)器顯示各因子存在明顯的交互增強(qiáng)效應(yīng),NDVI∩土地利用類型的交互作用最高,年均降水∩人口密度的交互作用最??;風(fēng)險(xiǎn)探測(cè)器表明不同分級(jí)區(qū)間各因子對(duì)礫石粒徑空間分布具有一定的差異性,其中土地利用類型具有最強(qiáng)的區(qū)間影響差異。

    (4)回歸分析結(jié)果顯示,最優(yōu)尺度回歸精度最高,人口密度未通過(guò)顯著性檢驗(yàn),相關(guān)系數(shù)由大到小為:NDVI、土地利用類型、年均降水、植被類型,且均與礫石粒徑呈正相關(guān),除年均降水與植被類型與地理探測(cè)器略有不同外,其余結(jié)果均與地理探測(cè)器結(jié)果相吻合。

    猜你喜歡
    分異礫石土地利用
    考慮礫石顆粒形狀及含量影響的砂-礫石混合物離散元模擬直剪試驗(yàn)
    Task 3
    土地利用生態(tài)系統(tǒng)服務(wù)研究進(jìn)展及啟示
    重慶市臭氧時(shí)空分異及其影響因素研究
    平泉縣下?tīng)I(yíng)坊雜巖體分異演化及其成巖成礦
    濱??h土地利用挖潛方向在哪里
    熱采井礫石充填防砂篩管外擠受力分析及應(yīng)用
    斷塊油氣田(2014年6期)2014-03-11 15:34:04
    北京市1989-2010年地表溫度時(shí)空分異特征分析
    北京山區(qū)礫石覆蓋度和礫石含量的關(guān)系研究
    論低碳經(jīng)濟(jì)與轉(zhuǎn)變土地利用方式
    黄频高清免费视频| 国产一区二区三区在线臀色熟女| 久久天躁狠狠躁夜夜2o2o| 男人和女人高潮做爰伦理| 国产日本99.免费观看| 国产综合懂色| 国产精品国产高清国产av| 国产极品精品免费视频能看的| 无遮挡黄片免费观看| 久久精品国产清高在天天线| 99国产精品一区二区蜜桃av| 国产精品亚洲av一区麻豆| 日本一本二区三区精品| 国产精品99久久99久久久不卡| 在线观看免费视频日本深夜| 亚洲自拍偷在线| 亚洲精品美女久久av网站| 91久久精品国产一区二区成人 | 久久久久久久精品吃奶| 亚洲一区高清亚洲精品| 亚洲成av人片在线播放无| 亚洲熟女毛片儿| 日韩欧美精品v在线| 久久精品影院6| 中亚洲国语对白在线视频| 日本撒尿小便嘘嘘汇集6| 久久久久久久久免费视频了| av福利片在线观看| 啦啦啦免费观看视频1| 麻豆国产97在线/欧美| 国产精品 国内视频| 夜夜夜夜夜久久久久| 午夜精品一区二区三区免费看| 成人国产综合亚洲| 精华霜和精华液先用哪个| 久久中文看片网| 窝窝影院91人妻| 天堂动漫精品| 亚洲国产欧美网| 两人在一起打扑克的视频| 欧美日韩福利视频一区二区| av片东京热男人的天堂| 国产v大片淫在线免费观看| 久久久久国产一级毛片高清牌| 欧美中文综合在线视频| 久久国产精品人妻蜜桃| 在线观看日韩欧美| 久久久久久久精品吃奶| 性色avwww在线观看| 成人三级做爰电影| 国产精品一区二区三区四区免费观看 | 国产探花在线观看一区二区| 嫩草影视91久久| 国产成人福利小说| 在线看三级毛片| 久久久久久久久中文| 欧美+亚洲+日韩+国产| 一二三四社区在线视频社区8| 美女cb高潮喷水在线观看 | 床上黄色一级片| tocl精华| 日韩欧美在线乱码| 亚洲国产精品sss在线观看| 久久中文字幕人妻熟女| 成人三级做爰电影| 国产精品一区二区精品视频观看| 欧美午夜高清在线| 国产欧美日韩一区二区精品| 色综合婷婷激情| 日韩人妻高清精品专区| 免费大片18禁| 久久伊人香网站| 在线a可以看的网站| 九九久久精品国产亚洲av麻豆 | 国产伦在线观看视频一区| 2021天堂中文幕一二区在线观| 欧美性猛交╳xxx乱大交人| e午夜精品久久久久久久| 亚洲午夜理论影院| 午夜福利欧美成人| 禁无遮挡网站| 黑人欧美特级aaaaaa片| 日本免费一区二区三区高清不卡| 偷拍熟女少妇极品色| 欧美黄色淫秽网站| 日本成人三级电影网站| 久久久水蜜桃国产精品网| 最近在线观看免费完整版| 啦啦啦观看免费观看视频高清| 亚洲成人久久性| 国产三级黄色录像| 国产单亲对白刺激| 熟女电影av网| 国产av不卡久久| 在线观看日韩欧美| 欧美日韩瑟瑟在线播放| 免费电影在线观看免费观看| 国产美女午夜福利| 在线十欧美十亚洲十日本专区| 亚洲国产看品久久| www.www免费av| 99热这里只有是精品50| 国产亚洲av高清不卡| 久久国产乱子伦精品免费另类| 久久这里只有精品中国| 精品乱码久久久久久99久播| 日本免费a在线| 男女视频在线观看网站免费| 亚洲色图 男人天堂 中文字幕| 99久久国产精品久久久| www.精华液| 国产成人精品久久二区二区免费| 天天一区二区日本电影三级| 午夜福利在线观看免费完整高清在 | 欧美午夜高清在线| 日本a在线网址| 精品99又大又爽又粗少妇毛片 | 精品熟女少妇八av免费久了| 日本熟妇午夜| 亚洲av免费在线观看| 最近在线观看免费完整版| 中文字幕高清在线视频| 国产又色又爽无遮挡免费看| 欧美xxxx黑人xx丫x性爽| 最近最新中文字幕大全电影3| 夜夜看夜夜爽夜夜摸| 在线免费观看的www视频| a在线观看视频网站| 亚洲人成网站高清观看| 欧美又色又爽又黄视频| 欧美黑人巨大hd| 国产aⅴ精品一区二区三区波| 国产99白浆流出| 久久精品综合一区二区三区| 国产视频内射| 免费无遮挡裸体视频| 国产91精品成人一区二区三区| 久久久久九九精品影院| 中文字幕高清在线视频| 国产伦精品一区二区三区视频9 | 国产真实乱freesex| 成人一区二区视频在线观看| 午夜福利18| 舔av片在线| 成人鲁丝片一二三区免费| 搡老熟女国产l中国老女人| 色噜噜av男人的天堂激情| 精品久久久久久久末码| 国产伦精品一区二区三区四那| 人人妻,人人澡人人爽秒播| 偷拍熟女少妇极品色| 婷婷丁香在线五月| 日韩中文字幕欧美一区二区| 天天添夜夜摸| 国产成人啪精品午夜网站| 成人特级黄色片久久久久久久| 亚洲激情在线av| 国产精华一区二区三区| 久久久久免费精品人妻一区二区| 婷婷六月久久综合丁香| 狂野欧美白嫩少妇大欣赏| 久久精品国产综合久久久| 美女免费视频网站| 久久精品91无色码中文字幕| 天堂av国产一区二区熟女人妻| 成人三级黄色视频| 亚洲色图 男人天堂 中文字幕| 99久久99久久久精品蜜桃| 亚洲美女黄片视频| 后天国语完整版免费观看| 国产精品野战在线观看| 国产又色又爽无遮挡免费看| 国产精品一区二区精品视频观看| 日本精品一区二区三区蜜桃| 少妇的丰满在线观看| 国产精品1区2区在线观看.| 久久久成人免费电影| 亚洲五月婷婷丁香| 中文资源天堂在线| 亚洲av美国av| 精品日产1卡2卡| 在线观看66精品国产| 长腿黑丝高跟| 欧美成人免费av一区二区三区| 国产精品99久久99久久久不卡| 变态另类丝袜制服| 久久久久国产一级毛片高清牌| 夜夜看夜夜爽夜夜摸| 一个人免费在线观看的高清视频| 国产精品 欧美亚洲| 人妻夜夜爽99麻豆av| 国产精品精品国产色婷婷| www.熟女人妻精品国产| 久久人人精品亚洲av| 啪啪无遮挡十八禁网站| 91在线观看av| 国产成年人精品一区二区| 非洲黑人性xxxx精品又粗又长| 少妇熟女aⅴ在线视频| 国产精品久久久久久精品电影| 国产麻豆成人av免费视频| 成人高潮视频无遮挡免费网站| 99精品欧美一区二区三区四区| 精品一区二区三区四区五区乱码| 亚洲精品在线美女| 91av网站免费观看| 亚洲 欧美一区二区三区| 美女免费视频网站| 男人和女人高潮做爰伦理| 国产免费男女视频| 国产精品九九99| 亚洲国产中文字幕在线视频| 老熟妇仑乱视频hdxx| 久久久久久久精品吃奶| 国产亚洲精品综合一区在线观看| 免费av毛片视频| 免费在线观看亚洲国产| 日本黄大片高清| 男女之事视频高清在线观看| 中文字幕av在线有码专区| 亚洲国产看品久久| 亚洲成a人片在线一区二区| 亚洲成人免费电影在线观看| 久久久精品欧美日韩精品| 伦理电影免费视频| 精品人妻1区二区| 18禁裸乳无遮挡免费网站照片| 丁香六月欧美| 18禁国产床啪视频网站| 99久久精品国产亚洲精品| 在线a可以看的网站| 最近视频中文字幕2019在线8| 最新美女视频免费是黄的| 午夜福利18| 男人和女人高潮做爰伦理| 国产精品乱码一区二三区的特点| 久久久久久大精品| 高潮久久久久久久久久久不卡| 亚洲国产精品sss在线观看| 亚洲欧美日韩高清专用| 在线观看免费视频日本深夜| 制服丝袜大香蕉在线| 男女那种视频在线观看| 国产主播在线观看一区二区| 热99re8久久精品国产| 亚洲色图av天堂| 欧美性猛交黑人性爽| 亚洲天堂国产精品一区在线| 巨乳人妻的诱惑在线观看| 1024手机看黄色片| 三级男女做爰猛烈吃奶摸视频| 九九热线精品视视频播放| 国产人伦9x9x在线观看| 巨乳人妻的诱惑在线观看| 可以在线观看的亚洲视频| 免费av不卡在线播放| 久久久精品欧美日韩精品| 国产在线精品亚洲第一网站| 国内揄拍国产精品人妻在线| 成人鲁丝片一二三区免费| 亚洲成人久久性| 国产精品久久视频播放| 激情在线观看视频在线高清| 国产精品女同一区二区软件 | 亚洲av日韩精品久久久久久密| 成人18禁在线播放| 男女下面进入的视频免费午夜| 狂野欧美白嫩少妇大欣赏| 久久香蕉精品热| 51午夜福利影视在线观看| 国产激情欧美一区二区| 亚洲美女视频黄频| 国产精品亚洲av一区麻豆| 国产高潮美女av| 嫩草影院入口| 国产毛片a区久久久久| 久久久久久国产a免费观看| 成人18禁在线播放| 啦啦啦免费观看视频1| 午夜激情福利司机影院| 麻豆久久精品国产亚洲av| 欧美在线一区亚洲| 国产成+人综合+亚洲专区| www.999成人在线观看| 三级毛片av免费| 亚洲专区字幕在线| 九九久久精品国产亚洲av麻豆 | 国产av麻豆久久久久久久| www.自偷自拍.com| 日韩欧美免费精品| 久久伊人香网站| 999久久久精品免费观看国产| 亚洲精品国产精品久久久不卡| 日本黄色片子视频| 国产黄片美女视频| 国产av不卡久久| 欧美中文日本在线观看视频| 三级男女做爰猛烈吃奶摸视频| 老汉色av国产亚洲站长工具| 非洲黑人性xxxx精品又粗又长| 精品久久久久久久人妻蜜臀av| 国产野战对白在线观看| 国产高清有码在线观看视频| 精品一区二区三区av网在线观看| 亚洲av成人不卡在线观看播放网| 制服丝袜大香蕉在线| 国产高清videossex| 成人av一区二区三区在线看| 亚洲九九香蕉| 久久伊人香网站| 两个人的视频大全免费| 听说在线观看完整版免费高清| 精品久久蜜臀av无| 老熟妇乱子伦视频在线观看| 国产免费av片在线观看野外av| 午夜亚洲福利在线播放| 校园春色视频在线观看| 不卡一级毛片| 午夜精品在线福利| 我的老师免费观看完整版| 这个男人来自地球电影免费观看| 99久久99久久久精品蜜桃| 久久这里只有精品中国| 真人一进一出gif抽搐免费| 亚洲中文日韩欧美视频| 欧美日本亚洲视频在线播放| 国产高清有码在线观看视频| 日韩欧美国产在线观看| 午夜两性在线视频| 人妻丰满熟妇av一区二区三区| 日本与韩国留学比较| 色播亚洲综合网| 久久久精品大字幕| 亚洲真实伦在线观看| 欧美黑人巨大hd| 国产欧美日韩精品一区二区| 村上凉子中文字幕在线| 亚洲精华国产精华精| 国产精品电影一区二区三区| 欧美av亚洲av综合av国产av| 18禁裸乳无遮挡免费网站照片| 国产成人精品久久二区二区免费| 村上凉子中文字幕在线| 国产日本99.免费观看| 嫩草影视91久久| 亚洲成人久久爱视频| 12—13女人毛片做爰片一| 亚洲五月婷婷丁香| 色哟哟哟哟哟哟| 国产亚洲av嫩草精品影院| 在线观看免费视频日本深夜| 禁无遮挡网站| 免费在线观看影片大全网站| 男插女下体视频免费在线播放| 国产亚洲av高清不卡| 在线观看午夜福利视频| 国产亚洲精品久久久久久毛片| 欧美成人一区二区免费高清观看 | 这个男人来自地球电影免费观看| 一区二区三区激情视频| 男插女下体视频免费在线播放| 麻豆成人av在线观看| 国产精品久久视频播放| 女同久久另类99精品国产91| 毛片女人毛片| 我的老师免费观看完整版| 男人和女人高潮做爰伦理| 国产极品精品免费视频能看的| 午夜福利视频1000在线观看| 国产极品精品免费视频能看的| 亚洲五月婷婷丁香| 人妻夜夜爽99麻豆av| 亚洲欧美日韩高清专用| 国产精品综合久久久久久久免费| 久久这里只有精品19| 国产精品久久久久久亚洲av鲁大| 成年免费大片在线观看| 亚洲成人中文字幕在线播放| 欧美极品一区二区三区四区| 一区福利在线观看| 亚洲第一电影网av| 精品乱码久久久久久99久播| 99热精品在线国产| 欧美性猛交黑人性爽| xxx96com| 18禁黄网站禁片免费观看直播| 国产成人aa在线观看| 久久香蕉精品热| 日本免费a在线| 一级毛片女人18水好多| 日韩 欧美 亚洲 中文字幕| 国产欧美日韩一区二区精品| 久久久久国产精品人妻aⅴ院| 99热只有精品国产| 老司机福利观看| 日韩免费av在线播放| 男女视频在线观看网站免费| 久久久久久久久中文| 欧美黄色片欧美黄色片| 99久国产av精品| 两个人的视频大全免费| 亚洲国产欧美一区二区综合| 一区二区三区国产精品乱码| 亚洲激情在线av| 国产v大片淫在线免费观看| 精品电影一区二区在线| 中亚洲国语对白在线视频| 国产亚洲精品av在线| 五月玫瑰六月丁香| 成人三级做爰电影| av中文乱码字幕在线| 国产真实乱freesex| 免费一级毛片在线播放高清视频| 国产成年人精品一区二区| svipshipincom国产片| 亚洲成人中文字幕在线播放| 亚洲人成电影免费在线| 夜夜看夜夜爽夜夜摸| 精品国产乱码久久久久久男人| 99久久综合精品五月天人人| 不卡一级毛片| 免费人成视频x8x8入口观看| 中文字幕av在线有码专区| 亚洲精品色激情综合| 亚洲avbb在线观看| 男女午夜视频在线观看| 亚洲国产色片| 黑人巨大精品欧美一区二区mp4| 国产高清有码在线观看视频| 99久久99久久久精品蜜桃| 亚洲熟妇中文字幕五十中出| 精品久久久久久成人av| av天堂在线播放| 成人性生交大片免费视频hd| 又紧又爽又黄一区二区| 观看美女的网站| 色综合婷婷激情| 麻豆国产97在线/欧美| netflix在线观看网站| 免费观看人在逋| 国产黄片美女视频| 1024手机看黄色片| 日本与韩国留学比较| 成熟少妇高潮喷水视频| 欧美黄色片欧美黄色片| 在线国产一区二区在线| 国产成人精品无人区| 国产精品亚洲一级av第二区| 欧美日韩一级在线毛片| 色在线成人网| svipshipincom国产片| 欧美日韩综合久久久久久 | 久久久久免费精品人妻一区二区| 99久久无色码亚洲精品果冻| 日本黄大片高清| 久久这里只有精品19| 18禁黄网站禁片免费观看直播| 丰满人妻熟妇乱又伦精品不卡| 久久久久久久久久黄片| 色综合婷婷激情| 岛国在线免费视频观看| 亚洲 欧美 日韩 在线 免费| 精品一区二区三区四区五区乱码| 男女下面进入的视频免费午夜| 国产aⅴ精品一区二区三区波| 国产一区二区在线av高清观看| 亚洲成av人片在线播放无| 欧美zozozo另类| 亚洲 欧美 日韩 在线 免费| aaaaa片日本免费| avwww免费| 91老司机精品| 亚洲av电影在线进入| 两性午夜刺激爽爽歪歪视频在线观看| 精品久久久久久久末码| 999久久久国产精品视频| 欧美色欧美亚洲另类二区| 国内揄拍国产精品人妻在线| 99热精品在线国产| 亚洲av第一区精品v没综合| 国语自产精品视频在线第100页| 男女视频在线观看网站免费| 老熟妇乱子伦视频在线观看| 一夜夜www| 国产精品乱码一区二三区的特点| 极品教师在线免费播放| 国产午夜精品论理片| 视频区欧美日本亚洲| 国产在线精品亚洲第一网站| 1024香蕉在线观看| 国产人伦9x9x在线观看| 久久精品国产亚洲av香蕉五月| 99久国产av精品| 国产av在哪里看| 亚洲成人久久爱视频| 1000部很黄的大片| www.www免费av| 日韩精品中文字幕看吧| 制服丝袜大香蕉在线| 一本综合久久免费| 嫩草影院精品99| 国产成人精品久久二区二区免费| 国产午夜精品久久久久久| 一级毛片女人18水好多| 女同久久另类99精品国产91| 夜夜爽天天搞| 免费无遮挡裸体视频| 天天躁日日操中文字幕| 亚洲电影在线观看av| 欧美色欧美亚洲另类二区| 国产精品久久久人人做人人爽| 国产精品野战在线观看| 亚洲黑人精品在线| 每晚都被弄得嗷嗷叫到高潮| 老司机福利观看| 十八禁网站免费在线| 老司机深夜福利视频在线观看| 亚洲五月婷婷丁香| 亚洲自拍偷在线| 九色国产91popny在线| 国产精品精品国产色婷婷| 一级黄色大片毛片| 很黄的视频免费| 男女做爰动态图高潮gif福利片| 国产av麻豆久久久久久久| 国产成人精品无人区| 成人av一区二区三区在线看| 人妻久久中文字幕网| 国产av不卡久久| 欧美精品啪啪一区二区三区| 久久久国产成人免费| 啦啦啦免费观看视频1| 怎么达到女性高潮| 亚洲七黄色美女视频| 三级国产精品欧美在线观看 | 久久精品国产清高在天天线| 韩国av一区二区三区四区| 俺也久久电影网| 男女午夜视频在线观看| 国产精品一区二区免费欧美| 黄频高清免费视频| 免费在线观看影片大全网站| 欧美极品一区二区三区四区| 国产精品九九99| 国语自产精品视频在线第100页| 国产美女午夜福利| 久久久成人免费电影| 亚洲av片天天在线观看| 18禁裸乳无遮挡免费网站照片| 最好的美女福利视频网| 日本一本二区三区精品| 18禁美女被吸乳视频| 三级国产精品欧美在线观看 | 免费高清视频大片| 少妇熟女aⅴ在线视频| 一级作爱视频免费观看| 午夜福利欧美成人| 天堂av国产一区二区熟女人妻| 亚洲国产欧美一区二区综合| 波多野结衣高清无吗| 一本精品99久久精品77| 少妇裸体淫交视频免费看高清| 日本黄大片高清| 色av中文字幕| 成人特级av手机在线观看| 麻豆成人av在线观看| 国产亚洲欧美98| 五月玫瑰六月丁香| 色老头精品视频在线观看| 国产美女午夜福利| 精品国产美女av久久久久小说| 国产午夜精品论理片| tocl精华| 午夜a级毛片| 亚洲中文日韩欧美视频| 嫁个100分男人电影在线观看| 欧美日本亚洲视频在线播放| 天堂网av新在线| 91久久精品国产一区二区成人 | 一区福利在线观看| 亚洲国产欧美网| 美女cb高潮喷水在线观看 | 淫妇啪啪啪对白视频| 国产成人精品久久二区二区91| 国产精品综合久久久久久久免费| 国产精品永久免费网站| 极品教师在线免费播放| av在线蜜桃| 精品久久久久久成人av| 嫩草影院精品99| 国产99白浆流出| 99热这里只有精品一区 | 国产成人av激情在线播放| 69av精品久久久久久| 欧美又色又爽又黄视频| 精品久久久久久久久久免费视频| 国产成人精品无人区| 国产精品影院久久| 精品一区二区三区av网在线观看| 亚洲国产精品成人综合色| 婷婷精品国产亚洲av在线| 亚洲人成网站高清观看| 91久久精品国产一区二区成人 | 在线永久观看黄色视频| 又黄又粗又硬又大视频| 亚洲成人久久爱视频| 久久中文看片网| www日本在线高清视频| 国产成人一区二区三区免费视频网站| 国产精品,欧美在线| 91九色精品人成在线观看| 99在线视频只有这里精品首页| 国产蜜桃级精品一区二区三区| 91在线观看av| 99在线视频只有这里精品首页|