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

    基于Landsat OLI的綠洲灌區(qū)土壤鹽度最優(yōu)預(yù)測尺度分析

    2017-10-14 02:21:17魏陽丁建麗王飛
    中國農(nóng)業(yè)科學(xué) 2017年15期
    關(guān)鍵詞:環(huán)境變量鹽漬化鹽度

    魏陽,丁建麗,王飛

    ?

    基于Landsat OLI的綠洲灌區(qū)土壤鹽度最優(yōu)預(yù)測尺度分析

    魏陽,丁建麗,王飛

    (新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院/綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室,烏魯木齊830046)

    【目的】基于遙感提取的多尺度遙感指示因子和土壤實(shí)測電導(dǎo)率數(shù)據(jù),借助統(tǒng)計(jì)分析,試圖探尋適合干旱區(qū)典型綠洲灌區(qū)土壤鹽度變異的最佳觀測尺度和指征變量,為快速評估綠洲土壤鹽漬化提供備選方案?!痉椒ā恳孕陆几珊?庫車河綠洲為研究區(qū),以野外采集的土壤鹽度數(shù)據(jù)(采集0—10、10—20、20—40和40—60 cm土層土樣,制備土壤飽和溶液并測試電導(dǎo)率(ms·cm-1))并將其作為預(yù)測對象(n=87),借助Landsat OLI遙感影像數(shù)據(jù),利用柵格重采樣(30—1 000 m)和領(lǐng)域?yàn)V波(原始分辨率為30、60、90、120、150、180、210 m,濾波尺度為3×3至31×31)兩種方式,生成多個(gè)尺度若干種指示因子(主成分分析、纓帽變化、植被指數(shù)、濕度指數(shù)),共計(jì)獲得1 078個(gè)(其中,柵格重采樣生成352個(gè),領(lǐng)域?yàn)V波生成726個(gè))環(huán)境變量。在此基礎(chǔ)之上,利用線性和非線性曲線模型分別擬合上述兩種模式下土壤鹽度和環(huán)境變量之間的相關(guān)性,進(jìn)而找出最優(yōu)環(huán)境因子和預(yù)測尺度?!窘Y(jié)果】柵格重采樣模式下能夠較好響應(yīng)各層土壤變異性的皆為非線性模式。其次,該模式下,擬合精度隨著空間分辨率的降低而降低。此模式下最佳推理尺度為30 m,該尺度下最佳響應(yīng)變量除了40—60 cm處為三波段差分指數(shù)(Three-band Maximal Gradient Difference,TGDVI)外,其余深度皆為擴(kuò)展的歸一化指數(shù)(Extended Normalized Difference Vegetation Index,ENDVI)。領(lǐng)域?yàn)V波模式下的最佳推理尺度為180 m(濾波尺度3×3),同時(shí),各層最佳擬合變量皆為擴(kuò)展的增強(qiáng)型植被指數(shù)(Extented Enhanced Vegetation Index,EEVI)。相比較柵格重采樣模式,該模式下的擬合精度全面優(yōu)于前者,各層依次提高14.60%、34.40%、32.10%和21.70%?!窘Y(jié)論】基于領(lǐng)域?yàn)V波模式下,像元分辨率為180 m,窗口大小為3×3的ENDVI指數(shù)更適合預(yù)測本研究區(qū)土壤鹽度的空間變異性。

    土壤鹽漬化;Landsat OLI; 尺度分析;綠洲灌區(qū);非線性模型;新疆

    0 引言

    【研究意義】McCarthy指出,在地理學(xué)研究中,不能期望在某一個(gè)尺度上研究得出的結(jié)論能適用于其他尺度,尺度上的每一個(gè)變化都會引出新的問題,沒有理由假設(shè)在某一個(gè)尺度上的關(guān)聯(lián)在其他尺度上仍然存在[1]。不同的自然現(xiàn)象有不同的最佳觀測距離和尺度,需要適當(dāng)?shù)木嚯x和比例尺,才能得到有效完整的觀察[2]。遙感是多尺度下土壤數(shù)字制圖研究中被廣泛利用的預(yù)測工具[3]。然而,此數(shù)據(jù)展現(xiàn)出的多種分辨率影響著決策者利用此數(shù)據(jù)的方式[4]。遙感數(shù)據(jù)通常呈現(xiàn)出一定的尺度依賴性,并以獨(dú)立的像元作為分析的基本單元。因此,遙感數(shù)據(jù)的尺度分析直接涉及空間分辨率[5]。依靠尺度效應(yīng),學(xué)者嘗試測量或者模擬下墊面過程,即不同的分辨率可能在捕捉期望中的變異性研究中,擁有不同的效能[6]。環(huán)境變量空間尺度分析與柵格分辨率和領(lǐng)域分析(窗口大小的選擇)有關(guān)?!厩叭搜芯窟M(jìn)展】目前已有很多學(xué)者調(diào)查分析柵格大小對土壤屬性的影響[7-8],以及與既定尺度下生物物理過程的相互作用關(guān)系[9-10],多數(shù)研究中選擇窗口的大小為3×3,其他窗口尺度并未考察。此外,部分學(xué)者也證明領(lǐng)域大小的選擇的重要性以及其對土壤-景觀關(guān)系的影響。統(tǒng)計(jì)模型是預(yù)測土壤屬性空間變異性研究中的常用方法。特定土壤屬性和環(huán)境變量的定量描述可以豐富本地的土壤-景觀知識庫。多數(shù)統(tǒng)計(jì)模型建立的前提是假設(shè)土壤-環(huán)境變量之間呈線性相關(guān)。然而,土壤屬性空間變異性與潛在環(huán)境變量之間的關(guān)系異常復(fù)雜,以線性假設(shè)恐難精確擬合二者之間的關(guān)系。借助現(xiàn)有遙感數(shù)據(jù)派生出的環(huán)境變量,利用多元統(tǒng)計(jì)分析技術(shù),已對土壤有機(jī)質(zhì)、紋理分析、土壤深度等屬性進(jìn)行了初步的探索[7-9]。而在土壤鹽度這一領(lǐng)域并未有相關(guān)研究。即便是在分布最為廣泛的干旱區(qū),如新疆,也未曾涉及。新疆維吾爾自治區(qū)堪稱干旱區(qū)鹽堿土博物館,分布范圍廣,種類多樣,作為西北最大的農(nóng)業(yè)用地儲備區(qū),快速精確地評估本地土壤鹽度則有利于本地水資源管理、農(nóng)業(yè)用地規(guī)劃和防止土壤鹽漬化的擴(kuò)張和反復(fù)。新疆灌區(qū)鹽漬化土地面積占灌區(qū)總面積的比例高達(dá)32.07%[11],多數(shù)城市和縣級地區(qū)都有分布。整體而言,以天山為界,南疆地區(qū)比北疆地區(qū)嚴(yán)重。鹽漬化土地主要分布在該區(qū)洪積扇中下部,綠洲-荒漠交錯帶,新開墾地區(qū)或老灌區(qū)內(nèi)部的棄耕地[12]。放眼未來,隨著人口的增加,人類活動將不斷加劇,及時(shí)升級土壤環(huán)境數(shù)據(jù)(物理化學(xué)、生物、生態(tài)水文等)的需求相比從前更為強(qiáng)烈[13]。新疆天然的地理環(huán)境,加之水土資源的不合理利用,勢必會造成更多的土地面臨土壤鹽漬化問題。另外,干旱區(qū)生態(tài)環(huán)境脆弱,準(zhǔn)確高效地診斷本地土壤鹽度對于評估土壤健康和優(yōu)化管理尤為重要。相對于傳統(tǒng)技術(shù)而言,遙感數(shù)據(jù)因其覆蓋面積廣,尋訪周期快,且非侵入土壤等優(yōu)勢得以廣泛應(yīng)用?!颈狙芯壳腥朦c(diǎn)】目前,借助遙感光譜探測土壤鹽度的方式主要有兩種。首先,在植被覆蓋度較低(植被覆蓋度小于15%)或鹽漬化較為嚴(yán)重的裸土地區(qū),可利用遙感原始波段光譜直接測量其土壤光譜信息[14-15],并與土壤鹽分進(jìn)行關(guān)聯(lián),定量獲取其鹽度值?;谀壳耙延醒芯筷P(guān)于地面光譜所總結(jié)的理論基礎(chǔ),許多學(xué)者嘗試借助鹽分敏感波段,建立相關(guān)鹽度指數(shù)定量描述土壤鹽度信息。其次,鹽漬化土壤因其鹽分含量超過正常閾值進(jìn)而影響到植被的生理參數(shù),在不同像元尺度上表現(xiàn)為,植被類型豐富度、植被出現(xiàn)頻率、植被蓋度降低的趨勢。鑒于此發(fā)現(xiàn),諸多研究學(xué)者利用植被指數(shù)間接推理土壤中的鹽分含量[16-17]。而上述研究多在單一尺度上(如運(yùn)用最為廣泛的Landsat數(shù)據(jù),30 m分辨率)進(jìn)行,并未綜合考察不同尺度下土壤鹽度-環(huán)境景觀之間的響應(yīng)關(guān)系?!緮M解決的關(guān)鍵問題】本研究以渭干河-庫車河綠洲為例,利用線性和非線性曲線擬合模型,分析不同尺度下土壤鹽度-環(huán)境變量之間的響應(yīng)關(guān)系,以期找出最佳的敏感尺度,提高基于遙感推理土壤鹽度的精度。

    1 研究區(qū)概況

    渭干河-庫車河綠洲(41°06′—41°40′N,82°10′—83°50′E)位于塔里木盆地北麓中段。典型的溫帶大陸性干旱氣候,日平均溫度變化劇烈(夏天最高溫度超過40℃),年平均降水量67.5 mm。由于極端的氣候,地質(zhì)環(huán)境和不合理的農(nóng)業(yè)管理的影響,導(dǎo)致本地土壤鹽漬化現(xiàn)象普遍。自然條件下植被覆蓋度較低,優(yōu)勢植被群落包括蘆葦(),檉柳(),駱駝刺(),花花柴()和鹽爪爪()。嚴(yán)重的地區(qū)出現(xiàn)鹽殼風(fēng)化現(xiàn)象()。相比阜北地區(qū)和瑪納斯河綠洲而言,該地區(qū)鹽漬化現(xiàn)象更為典型和嚴(yán)重。

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

    2.1 數(shù)據(jù)

    2.1.1 采樣數(shù)據(jù) 本研究使用的樣本采樣時(shí)間為2015年7月10—17日,樣本量為87個(gè),每個(gè)樣點(diǎn)采集5個(gè)樣本進(jìn)行混合。研究團(tuán)隊(duì)自2005年至今對渭干河-庫車河綠洲進(jìn)行多次采樣(每年至少一次),本次樣點(diǎn)布局充分考慮到研究區(qū)的土地利用方式和不同鹽度土壤的分布格局。樣本覆蓋以下幾種土地利用類型:農(nóng)業(yè)用地、鹽土、荒漠、低覆蓋草地、中等覆蓋草地、灌叢。每個(gè)樣點(diǎn)采集5個(gè)樣本,覆蓋面積范圍為30—250 m2。查閱新疆各個(gè)典型綠洲關(guān)于土壤鹽分的空間變程后發(fā)現(xiàn),自相關(guān)距離在1.14—27.45 km。而本研究最小分辨率為1 km,小于上述區(qū)間。所以,可以認(rèn)定1 km以內(nèi)的土壤鹽度是相對同質(zhì)的,或者變異性較小。土樣采集深度為:0—10、10—20、20—40和40—60 cm。土樣去除雜質(zhì),經(jīng)過2 mm篩后,以20 g土樣中加入100 mL蒸餾水進(jìn)行充分混合,待過濾后,利用土壤溶液電導(dǎo)率儀(Merck Millipore, Billerica, MA, USA)測量土壤溶液電導(dǎo)率(ms·cm-1)。

    2.1.2 遙感數(shù)據(jù) 為了更好的擬合圖像數(shù)據(jù)與地面實(shí)測土壤鹽度數(shù)據(jù),且減少大氣對最終結(jié)果的影響,研究利用遙感軟件ENVI5.1(ITT VIS, 2014)中的FLAASH(Fast Line-of-Sight Atmospheric Analysis of Spectral Hypercubes)模塊對原始數(shù)據(jù)進(jìn)行大氣校正。模型中的主要參數(shù)包括大氣模式(設(shè)定為中緯度夏天(Mid-Latitude Summer)),氣溶膠模式(選擇鄉(xiāng)村(Rural))。最終將輻射率轉(zhuǎn)化為地表反射率,為計(jì)算指數(shù)提供相對精確的數(shù)據(jù)支持。本文使用的遙感數(shù)據(jù)為Landsat 8 OLI,影像獲取時(shí)間為2015年7月15日,行列號為145/31。

    2.2 方法

    2.2.1 環(huán)境變量 因地理環(huán)境的不同,借助遙感數(shù)據(jù)構(gòu)建的各種指數(shù)在其研究區(qū)表現(xiàn)出對土壤鹽度的高敏感性,在其他地區(qū)是否依然能夠保持這種特性還有待商榷。目前已有多個(gè)植被和鹽度指數(shù)被應(yīng)用到世界各地測試間接推理土壤鹽度的潛力,其中植被指數(shù)包括NDVI[18]、EVI[19]、GARI[20]、CRSI[20]、COSRI[21]、GDVI[22]、SAVI[23]、TGDVI[24]、ENDVI[17]、EEVI[17],鹽度指數(shù)包括SI-T[25]、S2[26]、SI2[27]、SI3[27]、NDSI[28]。然而,僅BRUNNER等[18]利用NDVI 在新疆的博斯騰湖地區(qū)評價(jià)過土壤鹽度,其他指數(shù)均較少應(yīng)用到新疆地區(qū)。而探尋指數(shù)與土壤鹽度之間的定量關(guān)系,對于識別潛在的土壤鹽漬化土地以及制定相關(guān)利用措施有重要意義。

    另外,研究還選用了圖像變換衍生因子(主成分分析和纓帽變換)和土壤濕度指數(shù)。其中后者,即歸一化紅外指數(shù)(即濕度指數(shù))(Normalized Difference Infrared Index, NDII),不僅與冠層濕度有關(guān),還與鹽度有關(guān),原因在于融合了Landsat TM的第五波段信息。其他非植被指數(shù)類環(huán)境協(xié)同變量,主成分分析(Principal Component Analysis, PCA)的前三個(gè)波段(PCA1、PCA2、PCA3)占整體信息的90%以上和纓帽變換(Tasseled Cap, TC)。TC變換將陸地表面的信息轉(zhuǎn)化成若干個(gè)主題指數(shù),如亮度(Tasseled Cap Brightness, BRI),植被綠度(Tasseled Cap Green, GRE)和土壤濕度(Tasseled Cap Wet, WET)信息。其中亮度指數(shù)可近似作為土壤反照率或者反射率的相關(guān)指示因子。

    2.2.2 尺度變換 本文主要基于升尺度方法研究不同像元分辨率對土壤鹽分-植被指數(shù)之間關(guān)系的影響。研究將采用兩種方法進(jìn)行升尺度變換。一是利用ENVI 5.1中的重采樣模塊(RESIZE DATA)將30 m的反射率數(shù)據(jù)進(jìn)行重采樣,重采樣方法選擇最鄰近(Nearest Neighbor)算法,此方法相對于雙線性內(nèi)插法(Bilinear)和三次卷積內(nèi)插法(Cubic Convolution)而言,更接近真實(shí)值。最終將30 m的原始反射率數(shù)據(jù)依次重采樣成60、90、150、 210、270、330、390、 450、510、630、690、750、810、870、930和1 000 m,將重采樣后的反射率數(shù)據(jù)分別計(jì)算2.2.1節(jié)中提及的指數(shù)和衍生因子,如此共計(jì)產(chǎn)生374個(gè)環(huán)境變量。二是借助GRASS 7.1 的r.param.scale尺度轉(zhuǎn)換模塊,按照n×n移動窗口分別計(jì)算3×3,5×5,7×7,…,31×31多個(gè)領(lǐng)域?yàn)V波,窗口大小為奇數(shù),數(shù)值間隔為2,具體計(jì)算方法詳見圖1。紋理計(jì)算中,初步選擇了5種原始分辨率,分別是30、60、90、150和210 m,紋理尺度窗口則選擇了(3×3)—(31×31)共計(jì)15個(gè)級別,如此則對應(yīng)33個(gè)尺度(分辨率)(圖1)。每種分辨率下都是優(yōu)先進(jìn)行原始反射率數(shù)據(jù)的紋理計(jì)算,之后所形成的新波段數(shù)據(jù)再計(jì)算文中選擇的22個(gè)指數(shù)和主題因子(2.2.1),由此共計(jì)獲取726個(gè)環(huán)境變量。

    圖1 原始分辨率經(jīng)過濾后對應(yīng)的空間分辨率

    2.2.3 統(tǒng)計(jì)分析 采用線性和非線性相關(guān)性方式檢驗(yàn)和分析土壤鹽分指數(shù)和植被指數(shù)對土壤鹽度的敏感性,具體包括線性模型、對數(shù)模型、逆模型、二次函數(shù)、立方函數(shù)、復(fù)合函數(shù)、冪函數(shù)、S函數(shù)、增長函數(shù)、指數(shù)函數(shù)、邏輯模型。

    為了評價(jià)上述函數(shù)的擬合精度,研究引用決定系數(shù)2和顯著性檢驗(yàn)(Sig.)(2-tailed),定量評估每個(gè)指數(shù)與不同深度土壤鹽度之間的敏感性。其中Sig.<0.05為顯著性,Sig.<0.01為極度顯著。數(shù)據(jù)分析和檢驗(yàn)均在SPSS 22.0中的曲線擬合函數(shù)完成。

    3 結(jié)果

    3.1 土壤鹽度統(tǒng)計(jì)分析

    根據(jù)實(shí)驗(yàn)統(tǒng)計(jì)結(jié)果,表1顯示了各層數(shù)據(jù)的基本分布信息。以表層土壤鹽度值為例,根據(jù)美國鹽度實(shí)驗(yàn)室的等級劃分(非鹽漬化:0—2 ms·cm-1;輕度鹽漬化:2—4 ms·cm-1;中度鹽漬化:4—8 ms·cm-1;重度鹽漬化:8—16 ms·cm-1;極度鹽漬化:>16 ms·cm-1),從非鹽漬化至極度鹽漬化排列,樣本所占比例依次為16.00%、6.90%、19.50%、17.20%和40.40%,樣本土壤已涵蓋各個(gè)等級。自上而下觀察土壤剖面的鹽度均值,25%分位數(shù)(將各層鹽度值由小到大依次排列,87個(gè)樣本中1/4處的數(shù)值即25%分位數(shù)),中值(50%)和75%分位數(shù)的分布發(fā)現(xiàn),25%分位數(shù)的鹽度值已是中度鹽漬化,剖面鹽度呈顯表聚態(tài)勢,以此判定該區(qū)土壤鹽漬化情況較為嚴(yán)重且普遍。此外,由實(shí)測數(shù)據(jù)得知(數(shù)據(jù)表未顯示),陰離子C1-的平均含量明顯高于SO42-和CO32-的含量。原因在于,易溶性鹽類上行過程中,氯化物表聚性最為強(qiáng)烈,硫酸鹽次之,碳酸鹽最弱。由此可得出,此時(shí)該綠洲土壤鹽類主要以氯化物為主,同時(shí)還伴有少量的硫酸鹽-氯化物類。

    表1 研究區(qū)土壤電導(dǎo)率統(tǒng)計(jì)分布特征

    3.2 尺度變換對土壤鹽度-環(huán)境變量的影響

    圖2—9顯示了本研究所選環(huán)境變量在多個(gè)尺度下(柵格重采樣和領(lǐng)域?yàn)V波)解析土壤鹽度(不同土壤深度)變異性。需要說明的是,圖中所有的橫坐標(biāo)為不同尺度環(huán)境變量,即不同尺度下的植被指數(shù)、鹽度指數(shù)及其他衍生因子。因數(shù)量較多,在此僅以數(shù)字表達(dá),排列順序則按照30—1 000 m依次排開,變量順序不變。另外,繪制的解釋力圖中,研究只選擇了顯著性Sig.<0.01的變量,其他指數(shù)則未參與顯示。由圖2—5可知,參與計(jì)算的11個(gè)模型中,有6個(gè)模型擬合精度相對較高,包括邏輯模型,成長模型,復(fù)合模型,指數(shù)模型,S模型和冪函數(shù)。同時(shí),前4個(gè)模型最大擬合精度要高于后2個(gè)模型,此規(guī)律適用于除20—40 cm深度之外的其他土層。圖6—9中,橫坐標(biāo)的變量排列順序?yàn)?0 m的3×3至31×31,隨后為60 m的3×3至17×17,如此依次排開,變量順序則不變。除了表層復(fù)合模型缺席之外,上述6個(gè)模型都能夠捕捉到各層土壤鹽度-環(huán)境變量之間的響應(yīng)關(guān)系,且響應(yīng)尺度較多??傮w而言,0—10 cm、10—20 cm兩個(gè)土層處,30、60、90、150和210 m等5個(gè)原始分辨率計(jì)算的指數(shù)擬合精度要高于其領(lǐng)域?yàn)V波計(jì)算的指數(shù),且隨濾波尺度增大而降低,但不絕對。20—40和40—60 cm兩個(gè)土層的擬合結(jié)果與上述兩個(gè)土層比較,規(guī)律不明顯。

    圖2 不同尺度環(huán)境變量(柵格重采樣模式)對0—10 cm土層鹽度的解釋力

    圖3 不同尺度環(huán)境變量(柵格重采樣模式)對10—20 cm土層鹽度的解釋力

    圖4 不同尺度環(huán)境變量(柵格重采樣模式)對20—40 cm土層鹽度的解釋力

    3.3 不同深度最佳響應(yīng)尺度和變量

    表2總結(jié)了兩種柵格模式下擬合精度數(shù)值和模型參數(shù),并列出各層鹽度擬合精度最優(yōu)尺度和模型。柵格柵格重采樣模式中,空間分辨率為30 m下的變量集相對更為敏感。同時(shí),0—40 cm深度的最佳敏感指數(shù)皆為擴(kuò)展性植被指數(shù)ENDVI,而40—60 cm處則為植被指數(shù)TGDVI。10—20 cm處擬合精度最高,2=0.46,最低值為2=0.35 (20—40 cm)。領(lǐng)域?yàn)V波模式下,擬合精度最高的尺度為180 m(領(lǐng)域?yàn)V波為3×3),響應(yīng)變量則全為擴(kuò)展性植被指數(shù)ENDVI。最大解釋力出現(xiàn)于10—20 cm處,與柵格重采樣模式結(jié)果相似,能夠解釋58.10%的鹽度變異性,最小值則出現(xiàn)在20—40 cm處,解釋力為42.5%。

    圖5 不同尺度環(huán)境變量(柵格重采樣模式)對土壤鹽度(40—60 cm)的解釋力

    圖6 不同尺度環(huán)境變量(領(lǐng)域?yàn)V波)對(0—10 cm)土壤鹽度的解釋力

    在分析上述兩種模式下最優(yōu)土壤鹽度-ENDVI散點(diǎn)圖擬合關(guān)系時(shí)發(fā)現(xiàn)(圖10和圖11),當(dāng)去除極端案例時(shí),擬合精度會進(jìn)一步提升(研究區(qū)極端鹽漬化的覆蓋面積相對較少,去除極端案例后的擬合精度更為貼合實(shí)際情況)。柵格模式下,10—20和40—60 cm處2分別提升至47.90%和43.20%。領(lǐng)域?yàn)V波模式下,10—20、20—40和40—60 cm處2分別提升至64.39%、46.21%和52.60%。

    圖7 不同尺度環(huán)境變量(領(lǐng)域?yàn)V波)對(10—20 cm)層鹽度的解釋力

    圖8 不同尺度環(huán)境變量(領(lǐng)域?yàn)V波)對(20—40 cm)土層鹽度的解釋力

    圖9 不同尺度環(huán)境變量(領(lǐng)域?yàn)V波)對(40—60 cm)土層鹽度的解釋力

    表2 不同深度土壤鹽度-植被指數(shù)的最佳擬合模型

    :代表土壤鹽度Soil salinity

    圖10 去除極端案例下不同深度處土壤鹽度-ENDVI/TGDVI的擬合關(guān)系(柵格重采樣模式)

    圖11 去除極端案例下不同深度處土壤鹽度-ENDVI的擬合關(guān)系(領(lǐng)域?yàn)V波模式)

    4 討論

    隨著空間分辨率的降低,兩種模式下的土壤鹽度-環(huán)境變量擬合精度隨之降低。MULLER和van NIEKERK在其研究中對比了多種分辨率下(Rapideye, IRS, SPOT6/7, CBERS4, SENTINEL-2, ASTER)植被指數(shù)與土壤鹽度之間的響應(yīng)關(guān)系發(fā)現(xiàn),隨著分辨率的降低2值依次降低(0.22—0.66)[29]。本文中的領(lǐng)域?yàn)V波模式下以表層鹽度擬合關(guān)系為例,30 m的濾波由3×3至31×31依次計(jì)算,發(fā)現(xiàn)擬合精度同樣依次降低。但也不絕對,如180 m的 3×3尺度,響應(yīng)關(guān)系明顯提升,這與研究區(qū)的地理特征和土壤屬性的空間變異性有關(guān)。整體而言,隨著土層深度的增加,該趨勢則不明顯。由此看來,空間分辨率對擬合精度的影響較為明顯。

    若只用30 m大小的像元去捕捉鹽分-地表景觀之間的關(guān)系,可能會忽視由于尺度影響帶來的不確定性。不同傳感器獲取的固定像元面積就是描述地表景觀變異的最佳尺度嗎?目前多數(shù)研究的地面覆蓋的采樣面積以像元大小為標(biāo)準(zhǔn)。但土壤屬性的空間變異性與地表像元代表的光譜變異性不是一一對應(yīng)的,這就需要進(jìn)行研究,到底多大的像元或者空間尺度(由此產(chǎn)生的地表反射率,植被指數(shù),植被覆蓋度等特性)可以描述鹽分空間的變異性,這就是本文的主題。尺度效應(yīng)問題同樣受到其他學(xué)者的關(guān)注,黃彥等[30]對國內(nèi)外基于遙感獲取植被理化參數(shù)因尺度效應(yīng)產(chǎn)生的原因和空間異質(zhì)性進(jìn)行了綜述。分析眾多文獻(xiàn)后發(fā)現(xiàn),當(dāng)多種地表覆被類型同時(shí)出現(xiàn)在單位像元內(nèi)時(shí),即混合像元現(xiàn)象,會發(fā)生尺度效應(yīng),由此帶來的異質(zhì)性是產(chǎn)生尺度效應(yīng)的主要因素。不同的像元(尺度)內(nèi)的地表景觀結(jié)構(gòu)異質(zhì)性,因分辨率的大小不同程度影響著土壤-植被之間的關(guān)系,同時(shí)伴隨著一定的不確定性。本研究結(jié)果顯示,像元分辨率180 m是觀測本地土壤鹽分-植被變化的最優(yōu)尺度。當(dāng)研究區(qū)發(fā)生變化,地表異質(zhì)性隨之改變時(shí),上述尺度還是否是適宜的觀測尺度還需要進(jìn)一步研究。但給予的啟發(fā)是,基于遙感描述上述關(guān)系時(shí),應(yīng)嘗試多種尺度。由郝金標(biāo)等[31]的研究可以看出,水鹽變化導(dǎo)致植被的類型、蓋度、頻率發(fā)生一定的規(guī)律性變化。在干旱區(qū),異型性下墊面非常普遍,單位像元下的地表覆被類型各有不同。當(dāng)?shù)叵滤患巴寥乐械乃},符合干旱區(qū)植被最佳生存環(huán)境時(shí),單位面積下的植被類型較為豐富,且類型多樣,蓋度也會相應(yīng)增加,此時(shí)的異質(zhì)性也隨之加大,而生態(tài)水位下降時(shí),生存環(huán)境變得苛刻,適應(yīng)的植被類型減少,即便存在,也多半是耐旱性植被。此時(shí)的空間異質(zhì)性出現(xiàn)下降趨勢,地表植被類型單一。裸露的土壤面積加大,異質(zhì)性的主導(dǎo)因子則轉(zhuǎn)為土壤類型,或者地表土壤的物化屬性,涉及地表粗糙度,土壤濕度,土壤顏色等特征。那么針對此問題,單一的尺度未必能很好地詮釋不同尺度下由于下墊面異質(zhì)性而改變的土壤-植被間的關(guān)系。

    參與計(jì)算的各種指數(shù)涉及下墊面的相關(guān)屬性,包括植被覆蓋度、土壤濕度、地上植被生物量、土壤亮度、土壤反照率、土壤溫度等,但以ENDVI為代表的植被指數(shù)響應(yīng)關(guān)系最佳,相對而言,并成明顯的負(fù)相關(guān)關(guān)系。其原因可能在于,ENDVI除了NIR和RED波段之外,第三方波段信息的加入,增加了其預(yù)測土壤鹽度的效能。首先,短波紅外的引入增加了植被的信息維度;其次,可見光波段之間和短波紅外波段之間呈高度相關(guān)性;最后,NIR與可見光和短波紅外之間相關(guān)性較低[31]。另外,2 100—2 300 nm區(qū)具有豐富的植被信息,可以更好地反映植被覆蓋度和生長狀況[15]。此外,目前該指數(shù)在黃河口區(qū)效果良好,2值為0.53(=113)。本研究結(jié)果,以表層為例,在考慮極端案例的情況下,2值為0.43(柵格重采樣)和0.49(領(lǐng)域?yàn)V波),樣本量為87個(gè),而在不考慮極端案例時(shí),=86,10—60 cm各土層擬合關(guān)系中,2值都有提升,但領(lǐng)域?yàn)V波模式更為明顯(2=0.64)。與前人在該區(qū)(渭干河流域)研究對比發(fā)現(xiàn),DING等[32]在構(gòu)建鹽度模型時(shí),綜合考慮變量SI(salinity index, SI)和NDVI,擬合精度2=0.43(=68),但需要注意的是,該研究中使用時(shí)的是EM38大地電導(dǎo)率儀測得的表觀電導(dǎo)率,而非飽和溶液電導(dǎo)率值,二者的相關(guān)性為2值為0.663(濕季,7—8月)[31]。對此,需要進(jìn)一步研究ENDVI與表觀電導(dǎo)率之間的關(guān)系。

    研究結(jié)果顯示,與線性關(guān)系相比,非線性關(guān)系更能夠表達(dá)土壤鹽度-植被之間的定量關(guān)系。這與MULLER 和 van NIEKERK[29]的研究結(jié)果相似,在其研究中顯示,土壤鹽度與EVI(6 m)之間最佳擬合模型為復(fù)合模型,屬于非線性模式。FERNáNDEZ- BUCES 等在其研究中利用構(gòu)建的聯(lián)合光譜響應(yīng)指數(shù)(combined spectral response index,COSRI)[21],即NDVI改進(jìn)型指數(shù)(NDVI(B1+B2)/(B3+B4)),與飽和溶液電導(dǎo)率ECe(ms·cm-1)擬合后發(fā)現(xiàn),二者呈現(xiàn)非線性關(guān)系(=348.10e-18.37x,2=0.83),同NDVI相比,擬合關(guān)系更好,但其與鹽度之間的非線性關(guān)系并沒有改變。該指數(shù)在本文中并未使用,但在一定程度上說明了上述非線性關(guān)系的存在。在查閱相關(guān)文獻(xiàn)資料后發(fā)現(xiàn),新疆地區(qū)的土壤鹽分變化與植被物種豐富程度、物種多樣性、植被分布頻率之間關(guān)系較為復(fù)雜。郝金標(biāo)等[31]在其研究中指出,土壤鹽分變化與植被物種豐富程度之間呈非線性關(guān)系,與物種多樣性之間呈線性關(guān)系,與植被分布頻率之間因植被類型差異即有線性又有非線性關(guān)系存在。此外,石瑞花等[33]以新疆焉耆盆地為研究區(qū),分析該地天然植被與地下水關(guān)系后發(fā)現(xiàn),喬木、灌木、草本蓋度因土壤鹽分以及地下水礦化度的變化表現(xiàn)出一定的規(guī)律性特征,且皆成指數(shù)負(fù)相關(guān)關(guān)系。魚騰飛等[34]研究中提到,胡楊的平均密度與水溶性鹽分含量的關(guān)系并非呈線性關(guān)系。通過分析上述結(jié)果,發(fā)現(xiàn)植被蓋度或者頻率與土壤鹽分之間多以負(fù)相關(guān)出現(xiàn),并呈非線性關(guān)系。上述關(guān)系可能與地下水分布關(guān)系密切,而地下水位和礦化度之間又呈指數(shù)負(fù)相關(guān)。在渭干河-庫車河綠洲地區(qū),當(dāng)?shù)叵滤恍∮? m時(shí),地下水的礦化度增加趨勢明顯。3 m有可能就是鹽分聚集速率的分界線,當(dāng)?shù)叵滤淮笥? m時(shí)的積鹽速率要小于地下水位少于3 m的地區(qū)。此外,因干旱區(qū),降水稀少,蒸發(fā)強(qiáng)烈,多數(shù)自然植被依靠地下水位生存。不同植被的出現(xiàn)都與地下水的臨界水位有關(guān)。張麗等[35]以生態(tài)適宜性理論為基礎(chǔ),根據(jù)塔里木河干流流域典型植物的隨機(jī)抽樣調(diào)查資料,建立了干旱區(qū)幾種典型植物生長與地下水位關(guān)系的對數(shù)正態(tài)分布模型。由建立的模型得出干旱區(qū)典型植物的最適地下水位、適宜地下水位區(qū)間及其對環(huán)境因子的忍耐度。結(jié)論表明,塔里木河干流流域中的幾種代表性植被的最適宜的地下水位為1—3 m,而整個(gè)干旱區(qū)的適宜性生態(tài)水位2—4 m。在上述水位范圍內(nèi)的植被長勢良好,而地下水位大于4 m時(shí),則植被長勢受水分虧缺,或受土壤鹽漬化的影響,植被的多樣性、蓋度和頻率都會降低。

    本研究以目前應(yīng)用最為廣泛的Landsat數(shù)據(jù)為依托研究尺度對土壤鹽度-環(huán)境變量擬合關(guān)系的影響,其空間分辨率為30 m,而隨著傳感器分辨率的不斷提升,如worldview-2數(shù)據(jù)(9個(gè)波段,多光譜空間分辨率為2 m,全色波段為0.5 m),為下一步的研究提供新的數(shù)據(jù)源。本文中柵格重采樣的步長為30 m(像元逐個(gè)加入),以此計(jì)算的多種重采樣數(shù)據(jù)可能存在漏選的尺度,進(jìn)而隨之計(jì)算的擬合關(guān)系不能精確刻畫土壤鹽度-環(huán)境變量之間動態(tài)耦合趨勢。高空間分辨率的加入縮短了步長,為改進(jìn)擬合精度提供了可能性。

    5 結(jié)論

    本文基于遙感分析干旱區(qū)綠洲灌區(qū)土壤鹽度最優(yōu)預(yù)測尺度,初步得出以下結(jié)論。即,領(lǐng)域?yàn)V波模式和空間重采樣模式相比,易于提高土壤鹽度-環(huán)境變量之間的響應(yīng)關(guān)系。在采用的眾多指示因子中,以擴(kuò)展的歸一化植被指數(shù)(ENDVI)表現(xiàn)最佳。此外,適于描述本地土壤鹽度-環(huán)境變量關(guān)系的最佳觀測尺度為180 m,領(lǐng)域窗口大小為3×3。需要注意的是,因地理環(huán)境的差異性,應(yīng)多嘗試幾種尺度,以便更好地了解本地的土壤-景觀之間的響應(yīng)關(guān)系。

    References

    [1] MCCARTHY ALBERT J P.The Irish national electrification scheme., 1957, 47(4): 539-554.

    [2] 李小文, 王祎婷. 定量遙感尺度效應(yīng)芻議. 地理學(xué)報(bào), 2013, 68( 9) : 1163-1169.

    LI X W, WANG Y T. Prospects on future developments of quantitative remote sensing., 2013, 68(9): 1163-1169. (in Chinese)

    [3] Mulder V L, de Bruin S, Schaepman M E, Mayr T R. The use of remote sensing in soil and terrain mapping — a review., 2011, 162: 1-19.

    [4] Woodcock C E, Strahler A H. The factor of scale in remote sensing., 1987, 21: 311-332.

    [5] CAMPBELL J B, WYNNE R H.. New York: Guilford Press, 2011: 235-239.

    [6] MOELLERING H, TOBLER W. Geographical variances., 1972, 4 (1): 34-50.

    [7] BRADLEY A, SYLVIA K, MARC W, MICHAEL S. Impact of multi-scale predictor selection for modeling soil properties., 2015 (239/240): 97-106.

    [8] PACHEPSKY Y, HILL R. Scale and scaling in soils., 2017(287): 4-30.

    [9] KAMAL M, PHINN S JOHANSEN K. Assessment of multi- resolution image data for mangrove leaf area index mapping.2016(176): 242-254.

    [10] 聶志文, 布仁倉, 李冰倫. 紅松存在/不存在在數(shù)據(jù)的多尺度空間分布格局研究. 生態(tài)學(xué)報(bào), 2017, 37(3): 1-9.

    NIE Z W, BU R C, LI B L. Multi-sclale spatial pattern analysis of Korean pine based on presence/absence data., 2017, 37(3): 1-9. (in Chinese)

    [11] 王芳芳, 吳世新, 喬木, 李和平, 楊涵, 李義玲. 基于3S技術(shù)的新疆耕地鹽漬化狀況調(diào)查與分析. 干旱區(qū)研究, 2009, 26(3): 366-371.

    WANG F F, WU S X, QIAO M, LI H P, YANG H, LI Y L. Investigation and analysis on the salinization degree of cultivated land in Xinjiang based on 3S technology.2009, 26(3): 366-371. (in Chinese)

    [12] WANG F, CHEN X, LUO G P, HAN Q F. Mapping of regional soil salinities in Xinjiang and strategies for amelioration and management., 2015, 25(3): 321-336.

    [13] GRUNWALD S, THOMPSON J A, BOETTINGER J L. Digital soil mapping and modeling at continental scales: finding solutions for global issues., 2011, 75(4): 1201-1213.

    [14] KHAN N M, RASTOSKUEV V V, SATO Y, SHIOZAWA S. Assessment of hydrosaline land degradation by using a simple approach of remote sensing indicators., 2005, 77(1/3): 96-109.

    [15] DOUAOUI A E K, NICOLAS H, WALTER C. Detecting salinity hazards within a semiarid context by means of combining soil and remote-sensing data., 2006, 134(1/2): 217-230.

    [16] ZHANG T T, QI J G, GAO Y, OUYANG Z T, ZENG S L, ZHAO B. Detecting soil salinity with MODIS time series VI data., 2015, 52: 480-489.

    [17] 陳紅艷, 趙庚星, 陳敬春, 王瑞燕, 高明秀. 基于改進(jìn)植被指數(shù)的黃河口區(qū)鹽漬土鹽分遙感反演. 農(nóng)業(yè)工程學(xué)報(bào), 2015, 31(5): 107-114.

    CHEN H Y, ZHAO G X, CHEN J C, WANG R Y, GAO M X. Remote sensing inversion of saline soil salinity based on modified vegetation index in estuary area of Yellow River., 2015, 31(5): 107-114. (in Chinese)

    [18] BRUNNER P, LI H T, KINZELBACH W, LI W P. Generating soil electrical conductivity maps at regional level by integrating measurements on the ground and remote sensing data., 2007, 28(15): 3341-3361.

    [19] SCUDIERO E, SKAGGS T H, CORWIN D L. Regional-scale soil salinity assessment using Landsat ETM+ canopy reflectance., 2015, 169: 335-343.

    [20] SCUDIERO E, SKAGGS T H, CORWIN D L. Regional scale soil salinity evaluation using Landsat 7, western San Joaquin Valley, California, USA., 2014(2/3): 82-90.

    [21] FERNáNDEZ-BUCES N, SIEBE C, CRAM S, PALACIO J L. Mapping soil salinity using a combined spectral response index for bare soil and vegetation: A case study in the former lake Texcoco, Mexico., 2006, 65(4): 644-667.

    [22] WU W. The Generalized Difference Vegetation Index (GDVI) for dryland characterization., 2014, 6(2): 1211-1233.

    [23] ALLBED A, KUMAR L, ALDAKHEEL Y Y. Assessing soil salinity using soil salinity and vegetation indices derived from IKONOS high-spatial resolution imageries: Applications in a date palm dominated region., 2014(230/231): 1-8.

    [24] JIAPAER G, CHEN X, BAO A. A comparison of methods for estimating fractional vegetation cover in arid regions., 2011, 151: 1698-1710.

    [25] TRIPATHI N K, BRIJESH K R D P. Spatial modelling of soil alkalinity in GIS environment using IRS data. Paper presented at the 18th Asian Conference in Remote Sensing, 1997.

    [26] ABBAS A, KHAN S. Using remote sensing techniques for appraisal of irrigated soil salinity. Paper presented at the Advances and Applications for Management and Decision Making Land, Water and Environmental Management: Integrated Systems for Sustainability MODSIM07, 2007.

    [27] DOUAOUI A E K, NICOLAS H, WALTER C. Detecting salinity hazards within a semiarid context by means of combining soil and remote-sensing data., 2006, 134: 217-230.

    [28] KHAN N M, RASTOSKUEV V V, SATO Y, SHIOZAWA S. Assessment of hydrosaline land degradation by using a simple approach of remote sensing indicators., 2005, 77 (1): 96-109.

    [29] MULLER S J, VAN NIEKERK A. Identification of WorldView-2 spectral and spatial factors in detecting salt accumulation in cultivated fields., 2016, 273: 1-11.

    [30] 黃彥, 田慶久, 耿軍, 王磊, 欒海軍. 遙感反演植被理化參數(shù)的光譜和空間尺度效應(yīng). 生態(tài)學(xué)報(bào), 2016, 36(3): 883-891.

    HUANG Y, TIAN Q J, GENG J, WANG L, LUAN H J. Review of spectral and spatial scale effects of remotely sensed biophysical and biochemical vegetation parameters., 2016, 36(3): 883-891. (in Chinese)

    [31] 郗金標(biāo), 張福鎖, 毛達(dá)如, 田長彥, 董振成, 王開芳. 新疆鹽生植物群落物種多樣性及其分布規(guī)律的初步研究. 林業(yè)科學(xué), 2006, 42(10): 6-12.

    XI J B, ZHANG F S, MAO D R, TIAN C Y, DONG Z C, WANG K F. Species diversity and distribution of halophytic vegetation in Xinjiang., 2006, 42(10): 6-12. (in Chinese)

    [31] 徐涵秋, 唐菲. 新一代 Landsat 系列衛(wèi)星: Landsat 8 遙感影像新增特征及其生態(tài)環(huán)境意義.生態(tài)學(xué)報(bào), 2013, 33(11) : 3249-3257.

    XU H Q, TANG F. Analysis of new characteristics of the first Landsat 8 image and their eco-environmental significance., 2013, 33(11): 3249-3257. (in Chinese)

    [32] DING J L, YU D L. Monitoring and evaluating spatial variability of soil salinity in dry and wet seasons in the Werigan-Kuqa Oasis, China, using remote sensing and electromagnetic induction instruments., 2014(235/236): 316-322.

    [33] 石瑞花, 李霞, 董新光, 楊鵬年, 劉磊. 焉耆盆地天然植被與地下水關(guān)系研究. 自然資源學(xué)報(bào), 2009, 24(12): 2097-2103.

    SHI R H, LI X, DONG X G, YANG P N, LIU L. Research on the relationship between natural vegetation growth and ground water in Yanqi Basin., 2009, 24(12): 2097-2013. (in chinese)

    [34] 魚騰飛, 馮起, 劉蔚, 司建華, 席海洋, 陳麗娟. 黑河下游土壤水鹽對生態(tài)輸水的響應(yīng) 及其與植被生長的關(guān)系. 生態(tài)學(xué)報(bào), 2012, 32(22): 7009-7017.

    YU T F, FENG Q, LIU W, SI J H, XI H Y, CHEN L J. Soil water and salinity in response to water deliveries and the relationship with plant growth at the lower reaches of Heihe River, Northwestern China., 2012, 32(22): 7009-7017. (in chinese)

    [35] 張麗, 董增川, 黃曉玲. 干旱區(qū)典型植被生長與地下水位關(guān)系的模型研究. 中國沙漠, 2004, 24(1): 110-113.

    ZHANG L, DONG Z C, HUANG X L. Modelling on relation between major plants growth and ground water depth in arid area., 2004, 24(1): 110-113. (in Chinese)

    (責(zé)任編輯 李云霞)

    Optimal Scale Analysis of Soil Salinity Prediction in Oasis Irrigated Area of arid Land Based on Landsat OLI

    WEI Yang,DING JianLi, WANG Fei

    (College of Resource and Environmental Science, Xinjiang University/ Laboratory of Oasis Ecosystems, Ministry of Education, Urumqi 830046)

    【Objective】 Based on the multi-scale remote sensing indicators and measured soil conductivity data, the objective of this study is to exploit the appropriate scales and indicators for inferring soil salinity in irrigation area of arid oasis, and to provide alternative schemes for rapid assessment of soil salinization in the study area. 【Method】Weigan-kuqa oasis located in southern Xinjiang was selected as the study area. Soil electrical conductivity of multiple soil layers (0-10 cm, 10-20 cm, 20-40 cm, 40-60 cm), as inference object, were analyzed. Raster resampling (30-1 000 m) and neighborhood extent (3×3 to 31×31) were employed as the method of scale transformation to generate multiple scales of environmental factors (with help of principal component analysis, tasseled cap, vegetation index, soil wet index) based on Landsat OLI image data, then, a total of 1 078 environmental variables (the former produced 352 variables, the latter produced 726 variables). On this basis, linear and nonlinear curve models were used to fit the correlation between soil salinity and environmental variables under the above mentioned method of scale transformation, and then the optimal prediction scale and environmental factor were explored. 【Result】Nonlinear curve models showed a significance between soil salinity and environmental factors compared to linear model under multiple scale (raster resampling model). Secondly, the fitting accuracy decreased when the spatial resolution become coarse. In this model, the best inference scale was 30 m, and the sensitive variable was three-band maximal gradient difference (TGDVI) at depth of 40-60 cm, and the remaining depths were implied by Extended Normalized Difference Vegetation index (ENDVI). The optimal inference scale in model of neighborhood extent was 180 m (filtering scale 3×3), and the best fit variables for each layer was EEVI. Compared with the grid resampling mode, the fitting accuracy of neighborhood extent was improved by 14.60%, 34.40%, 32.10% and 21.70% (from top layer to bottom one) compared to fomer model, respectively. 【Conclusion】Based on the model of neighborhood extent, the nonlinear model constructed by Extented Enhanced Vegetation Index(EEVI) (180 m, 3×3) is more suitable for the prediction of soil salinity variation in the study area.

    soil salinization; Landsat OLI; scale analysis; oasis irrigation area; nonlinear model; Xinjiang

    2016-12-24;接受日期:2017-03-31

    國家自然科學(xué)基金(41661046)、中國博士后面上基金(2016M602909)、自治區(qū)科技支疆項(xiàng)目(201591101)、新疆大學(xué)博士啟動基金(BS150248)、新疆維吾爾自治區(qū)重點(diǎn)實(shí)驗(yàn)室專項(xiàng)基金(2014KL005)、國家自然科學(xué)基金(新疆聯(lián)合基金本地優(yōu)秀青年培養(yǎng)專項(xiàng)(U1503302))

    魏陽,E-mail:wcyang0909@163.com。通信作者丁建麗,E-mail:Watarid@xju.edu.cn

    猜你喜歡
    環(huán)境變量鹽漬化鹽度
    蔬菜大棚土壤鹽漬化成因及防治措施
    土地質(zhì)量地球化學(xué)調(diào)查成果在判定土壤鹽漬化、沙化中的應(yīng)用
    甘肅蘇干湖濕地土壤鹽漬化、地下水位埋深及其對生態(tài)環(huán)境的影響
    從桌面右鍵菜單調(diào)用環(huán)境變量選項(xiàng)
    徹底弄懂Windows 10環(huán)境變量
    瑪納斯河流域土壤鹽漬化時(shí)空動態(tài)變化
    基于三階段DEA—Malmquist模型的中國省域城鎮(zhèn)化效率測度及其收斂分析
    鹽度和pH對細(xì)角螺耗氧率和排氨率的影響
    鹽度脅迫對入侵生物福壽螺的急性毒性效應(yīng)
    適用于高鹽度和致密巖層驅(qū)油的表面活性劑
    欧美激情 高清一区二区三区| 高清欧美精品videossex| www.999成人在线观看| 欧美亚洲日本最大视频资源| 不卡av一区二区三区| 正在播放国产对白刺激| 免费在线观看日本一区| 成人永久免费在线观看视频| 老熟妇乱子伦视频在线观看| 午夜视频精品福利| videos熟女内射| 久久久精品免费免费高清| 亚洲中文字幕日韩| 好男人电影高清在线观看| 国产免费男女视频| 精品国产超薄肉色丝袜足j| 大型黄色视频在线免费观看| 天堂动漫精品| a在线观看视频网站| 99热网站在线观看| 啦啦啦在线免费观看视频4| 亚洲国产精品合色在线| 黑人巨大精品欧美一区二区蜜桃| 在线免费观看的www视频| 丝袜美腿诱惑在线| 制服诱惑二区| 国精品久久久久久国模美| 一区二区三区激情视频| 国产精品免费视频内射| ponron亚洲| 国产真人三级小视频在线观看| av中文乱码字幕在线| 两个人免费观看高清视频| 国产一区二区三区综合在线观看| av超薄肉色丝袜交足视频| 男女床上黄色一级片免费看| 久久人妻福利社区极品人妻图片| 波多野结衣一区麻豆| 老司机福利观看| 精品少妇久久久久久888优播| 高清欧美精品videossex| 侵犯人妻中文字幕一二三四区| 12—13女人毛片做爰片一| 在线观看一区二区三区激情| 国产免费现黄频在线看| 91国产中文字幕| e午夜精品久久久久久久| 搡老乐熟女国产| 在线观看舔阴道视频| 一边摸一边做爽爽视频免费| 亚洲片人在线观看| 久久国产精品大桥未久av| a级毛片在线看网站| 午夜精品久久久久久毛片777| 亚洲免费av在线视频| √禁漫天堂资源中文www| 久久 成人 亚洲| 岛国在线观看网站| 久久精品国产a三级三级三级| 日韩精品免费视频一区二区三区| 老司机影院毛片| tocl精华| 制服人妻中文乱码| 97人妻天天添夜夜摸| 日韩欧美国产一区二区入口| 国产不卡一卡二| 国产免费现黄频在线看| 亚洲人成电影免费在线| 午夜成年电影在线免费观看| 国产免费男女视频| 999久久久精品免费观看国产| 国产成人系列免费观看| 国产男女超爽视频在线观看| 欧美日韩福利视频一区二区| 久久久久精品国产欧美久久久| 亚洲av成人一区二区三| 久久久久久人人人人人| 欧美激情 高清一区二区三区| 欧美激情高清一区二区三区| 免费在线观看黄色视频的| 欧美国产精品一级二级三级| 成年人午夜在线观看视频| 日本黄色视频三级网站网址 | 一级片'在线观看视频| 国产成人欧美在线观看 | 日韩熟女老妇一区二区性免费视频| 狠狠狠狠99中文字幕| 丝袜人妻中文字幕| 人妻久久中文字幕网| 精品午夜福利视频在线观看一区| 久久午夜综合久久蜜桃| 亚洲精品av麻豆狂野| 精品久久蜜臀av无| 精品一品国产午夜福利视频| 国产成人欧美| 麻豆成人av在线观看| 欧美激情久久久久久爽电影 | 天堂中文最新版在线下载| 国产成人精品在线电影| 久久香蕉国产精品| 99re在线观看精品视频| 他把我摸到了高潮在线观看| av视频免费观看在线观看| 成人av一区二区三区在线看| 老熟妇仑乱视频hdxx| 亚洲精品美女久久av网站| 久久精品国产亚洲av香蕉五月 | 大片电影免费在线观看免费| 岛国在线观看网站| 一区福利在线观看| 老司机影院毛片| 国产精品美女特级片免费视频播放器 | 国产成人精品久久二区二区91| 国产97色在线日韩免费| 国产99久久九九免费精品| 欧美日韩精品网址| 亚洲片人在线观看| 久久人妻熟女aⅴ| 超色免费av| 欧美成人免费av一区二区三区 | 精品一区二区三区视频在线观看免费 | 中国美女看黄片| 一级片'在线观看视频| 国产成+人综合+亚洲专区| 老熟女久久久| 国产av精品麻豆| 国产精品美女特级片免费视频播放器 | 一级a爱视频在线免费观看| 国产精品98久久久久久宅男小说| 亚洲精品一卡2卡三卡4卡5卡| 手机成人av网站| 一边摸一边做爽爽视频免费| 国产一区二区三区视频了| 国产精品1区2区在线观看. | 亚洲精品久久成人aⅴ小说| 久久久国产欧美日韩av| tube8黄色片| 欧美 日韩 精品 国产| 三上悠亚av全集在线观看| 亚洲精品国产区一区二| 亚洲欧美精品综合一区二区三区| 国产精品1区2区在线观看. | av天堂久久9| 亚洲人成伊人成综合网2020| 性少妇av在线| 中出人妻视频一区二区| 亚洲国产精品一区二区三区在线| 亚洲精品中文字幕一二三四区| 亚洲国产欧美一区二区综合| 亚洲男人天堂网一区| 亚洲久久久国产精品| 亚洲成人国产一区在线观看| 99国产精品免费福利视频| 国产免费现黄频在线看| 国产精品一区二区在线观看99| 免费观看a级毛片全部| 精品久久久久久久久久免费视频 | 久久国产亚洲av麻豆专区| 美女福利国产在线| 中文字幕制服av| 欧美日本中文国产一区发布| 亚洲人成电影免费在线| 成人精品一区二区免费| 精品人妻1区二区| 日韩熟女老妇一区二区性免费视频| 国产极品粉嫩免费观看在线| 两个人看的免费小视频| 成人特级黄色片久久久久久久| 国产高清激情床上av| 亚洲国产精品一区二区三区在线| 国产精品欧美亚洲77777| 午夜福利影视在线免费观看| videos熟女内射| 亚洲,欧美精品.| 亚洲国产欧美一区二区综合| 欧美久久黑人一区二区| 午夜两性在线视频| 亚洲av成人av| 亚洲欧美激情综合另类| 午夜91福利影院| 飞空精品影院首页| 国产一区在线观看成人免费| 在线观看午夜福利视频| 大片电影免费在线观看免费| 激情在线观看视频在线高清 | 热99re8久久精品国产| 两个人看的免费小视频| 最近最新中文字幕大全电影3 | 日韩免费高清中文字幕av| 不卡av一区二区三区| 久久久久久免费高清国产稀缺| 午夜免费成人在线视频| av视频免费观看在线观看| 女人久久www免费人成看片| a级毛片在线看网站| cao死你这个sao货| 欧美精品一区二区免费开放| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲一区高清亚洲精品| 可以免费在线观看a视频的电影网站| 免费在线观看黄色视频的| 俄罗斯特黄特色一大片| 51午夜福利影视在线观看| 亚洲精华国产精华精| 免费人成视频x8x8入口观看| 99riav亚洲国产免费| 成人黄色视频免费在线看| 无人区码免费观看不卡| 欧美国产精品一级二级三级| 欧美中文综合在线视频| 淫妇啪啪啪对白视频| 又黄又粗又硬又大视频| 久久久精品区二区三区| 999久久久国产精品视频| 日韩精品免费视频一区二区三区| 国产精品久久久久久人妻精品电影| 亚洲综合色网址| 9色porny在线观看| 国产99白浆流出| 亚洲精品美女久久av网站| 欧美人与性动交α欧美精品济南到| 人妻丰满熟妇av一区二区三区 | 九色亚洲精品在线播放| 久久精品国产亚洲av香蕉五月 | 搡老熟女国产l中国老女人| 老熟女久久久| 国产精品 国内视频| 两人在一起打扑克的视频| 色播在线永久视频| 在线视频色国产色| videosex国产| 老司机午夜十八禁免费视频| 天天躁日日躁夜夜躁夜夜| 亚洲成人免费av在线播放| 久久国产精品大桥未久av| 精品欧美一区二区三区在线| 亚洲精品成人av观看孕妇| 免费在线观看完整版高清| 桃红色精品国产亚洲av| 免费一级毛片在线播放高清视频 | 美女扒开内裤让男人捅视频| 真人做人爱边吃奶动态| 国产免费男女视频| 欧美日韩精品网址| 一边摸一边做爽爽视频免费| 久久婷婷成人综合色麻豆| 欧美+亚洲+日韩+国产| 欧美+亚洲+日韩+国产| 国产蜜桃级精品一区二区三区 | av视频免费观看在线观看| 黑丝袜美女国产一区| 丝袜美腿诱惑在线| 两性夫妻黄色片| 狂野欧美激情性xxxx| 大码成人一级视频| 欧美久久黑人一区二区| 最近最新中文字幕大全电影3 | 国产精品免费视频内射| 国产无遮挡羞羞视频在线观看| 999久久久精品免费观看国产| 精品一区二区三卡| 免费在线观看完整版高清| 成人亚洲精品一区在线观看| 午夜精品在线福利| 久久久久精品国产欧美久久久| 中文字幕精品免费在线观看视频| 在线观看一区二区三区激情| 香蕉丝袜av| 久久香蕉国产精品| 国产淫语在线视频| 国产精品一区二区在线观看99| 亚洲一码二码三码区别大吗| 免费不卡黄色视频| 成人18禁高潮啪啪吃奶动态图| 王馨瑶露胸无遮挡在线观看| 国产视频一区二区在线看| 999精品在线视频| 精品视频人人做人人爽| 女性生殖器流出的白浆| 麻豆av在线久日| 久久天躁狠狠躁夜夜2o2o| 一级作爱视频免费观看| av网站免费在线观看视频| 好男人电影高清在线观看| 亚洲 国产 在线| 欧美在线黄色| 人人妻人人澡人人爽人人夜夜| 亚洲精华国产精华精| e午夜精品久久久久久久| 免费一级毛片在线播放高清视频 | 国产人伦9x9x在线观看| 欧美日韩国产mv在线观看视频| 老司机靠b影院| 久久久久久人人人人人| 母亲3免费完整高清在线观看| 亚洲人成伊人成综合网2020| 亚洲欧美精品综合一区二区三区| 女性生殖器流出的白浆| 久久久久精品国产欧美久久久| 身体一侧抽搐| 欧美日韩乱码在线| 国产精品偷伦视频观看了| 亚洲综合色网址| 久久影院123| 成人亚洲精品一区在线观看| 热99re8久久精品国产| 大型av网站在线播放| 在线观看免费午夜福利视频| 一级毛片女人18水好多| 超碰97精品在线观看| 国产一区二区三区综合在线观看| 精品熟女少妇八av免费久了| xxxhd国产人妻xxx| 精品国产亚洲在线| 精品国产一区二区三区四区第35| 亚洲熟妇中文字幕五十中出 | 国产99白浆流出| 欧美激情 高清一区二区三区| 久久人妻熟女aⅴ| 日本撒尿小便嘘嘘汇集6| 狂野欧美激情性xxxx| 好看av亚洲va欧美ⅴa在| 国产男靠女视频免费网站| 午夜福利免费观看在线| 成年人黄色毛片网站| 一级毛片女人18水好多| 黑人操中国人逼视频| 黑人欧美特级aaaaaa片| 一级a爱视频在线免费观看| 久久久国产欧美日韩av| 捣出白浆h1v1| 亚洲国产欧美一区二区综合| 国产精品免费大片| 热99国产精品久久久久久7| 热re99久久精品国产66热6| av欧美777| 交换朋友夫妻互换小说| 老司机福利观看| 中文字幕精品免费在线观看视频| 亚洲美女黄片视频| 99精品久久久久人妻精品| 亚洲av第一区精品v没综合| 国产欧美日韩综合在线一区二区| 老司机在亚洲福利影院| 欧美日韩中文字幕国产精品一区二区三区 | 国产不卡一卡二| 精品无人区乱码1区二区| 亚洲专区国产一区二区| 91av网站免费观看| 亚洲精品av麻豆狂野| 日韩精品免费视频一区二区三区| 两人在一起打扑克的视频| 婷婷成人精品国产| 视频在线观看一区二区三区| www日本在线高清视频| 天堂√8在线中文| 久99久视频精品免费| 免费一级毛片在线播放高清视频 | 成人av一区二区三区在线看| 久久精品熟女亚洲av麻豆精品| 精品久久久久久久毛片微露脸| 如日韩欧美国产精品一区二区三区| av一本久久久久| 国产精品免费视频内射| 国产精品久久久av美女十八| 男女下面插进去视频免费观看| 精品久久久久久电影网| 亚洲精品av麻豆狂野| 999精品在线视频| 亚洲少妇的诱惑av| 午夜亚洲福利在线播放| 一级毛片女人18水好多| 亚洲第一青青草原| 国产成人av教育| 精品少妇一区二区三区视频日本电影| 精品乱码久久久久久99久播| 欧美成人免费av一区二区三区 | 妹子高潮喷水视频| 午夜亚洲福利在线播放| 成人影院久久| 妹子高潮喷水视频| 人人妻人人爽人人添夜夜欢视频| 午夜福利免费观看在线| 麻豆国产av国片精品| 国产免费男女视频| 亚洲精品久久成人aⅴ小说| 男人操女人黄网站| 咕卡用的链子| 另类亚洲欧美激情| 免费在线观看亚洲国产| 久久香蕉精品热| 国内毛片毛片毛片毛片毛片| 欧美大码av| 欧美av亚洲av综合av国产av| 欧美日韩视频精品一区| 天堂俺去俺来也www色官网| 久久久久国内视频| 精品少妇一区二区三区视频日本电影| 欧美成狂野欧美在线观看| 999久久久精品免费观看国产| 久久精品亚洲av国产电影网| 国产欧美日韩一区二区三区在线| 99热国产这里只有精品6| 18禁观看日本| 18禁裸乳无遮挡动漫免费视频| 精品久久久久久,| 最新在线观看一区二区三区| tube8黄色片| 欧美日韩亚洲国产一区二区在线观看 | 91麻豆精品激情在线观看国产 | 亚洲精品在线美女| 中文字幕人妻熟女乱码| 国产精品影院久久| 18禁裸乳无遮挡动漫免费视频| tube8黄色片| 欧美日韩亚洲国产一区二区在线观看 | 国产成人精品久久二区二区91| 国产精品国产高清国产av | 一进一出抽搐gif免费好疼 | 在线国产一区二区在线| 国产有黄有色有爽视频| 99国产精品一区二区三区| 国产成人精品无人区| 欧美激情久久久久久爽电影 | 人人妻人人澡人人爽人人夜夜| 99久久99久久久精品蜜桃| 久久中文字幕一级| 美女国产高潮福利片在线看| 丝瓜视频免费看黄片| 1024视频免费在线观看| 99久久99久久久精品蜜桃| 免费在线观看完整版高清| 狂野欧美激情性xxxx| av天堂在线播放| 男人操女人黄网站| 免费日韩欧美在线观看| 捣出白浆h1v1| 婷婷丁香在线五月| 国产成人欧美| 啪啪无遮挡十八禁网站| 黑人操中国人逼视频| 美女高潮喷水抽搐中文字幕| 可以免费在线观看a视频的电影网站| 亚洲国产欧美日韩在线播放| 老鸭窝网址在线观看| 精品卡一卡二卡四卡免费| a级片在线免费高清观看视频| 777米奇影视久久| 国产又色又爽无遮挡免费看| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品aⅴ一区二区三区四区| 亚洲欧洲精品一区二区精品久久久| 人妻一区二区av| 午夜免费鲁丝| 亚洲人成77777在线视频| 在线观看免费高清a一片| 极品教师在线免费播放| 色婷婷久久久亚洲欧美| 国产亚洲欧美在线一区二区| 欧美国产精品一级二级三级| 国产无遮挡羞羞视频在线观看| 999精品在线视频| 美女高潮喷水抽搐中文字幕| 亚洲成a人片在线一区二区| 亚洲欧美日韩另类电影网站| 欧美日韩亚洲综合一区二区三区_| 看片在线看免费视频| 一a级毛片在线观看| 精品无人区乱码1区二区| 女人被躁到高潮嗷嗷叫费观| 日韩人妻精品一区2区三区| 亚洲专区字幕在线| 叶爱在线成人免费视频播放| 一级毛片高清免费大全| 一二三四社区在线视频社区8| www.熟女人妻精品国产| 欧美精品高潮呻吟av久久| 久99久视频精品免费| 12—13女人毛片做爰片一| 中亚洲国语对白在线视频| 亚洲自偷自拍图片 自拍| 亚洲精品美女久久av网站| 黄片播放在线免费| 精品一品国产午夜福利视频| 亚洲欧美激情综合另类| 12—13女人毛片做爰片一| 一区在线观看完整版| 亚洲熟女毛片儿| 日韩人妻精品一区2区三区| 欧美乱色亚洲激情| 国产成人精品在线电影| 怎么达到女性高潮| 在线观看免费视频网站a站| 日本一区二区免费在线视频| 天堂俺去俺来也www色官网| 国产精品九九99| 久久影院123| 亚洲国产毛片av蜜桃av| 三级毛片av免费| 久久久久久免费高清国产稀缺| 国产成人精品无人区| 色播在线永久视频| 天堂√8在线中文| av一本久久久久| 国产精品久久久av美女十八| 多毛熟女@视频| 男女床上黄色一级片免费看| 国产99久久九九免费精品| 天堂动漫精品| 国产成人精品久久二区二区免费| 男女之事视频高清在线观看| 黄色成人免费大全| 国产欧美日韩一区二区精品| 涩涩av久久男人的天堂| 亚洲欧洲精品一区二区精品久久久| 国产午夜精品久久久久久| 9191精品国产免费久久| 最近最新中文字幕大全免费视频| 久久国产亚洲av麻豆专区| 啦啦啦 在线观看视频| 欧美精品av麻豆av| 日韩精品免费视频一区二区三区| 色婷婷av一区二区三区视频| 精品一区二区三区视频在线观看免费 | 精品久久久久久电影网| 亚洲成人免费电影在线观看| 久久中文字幕人妻熟女| 我的亚洲天堂| 精品国产一区二区三区四区第35| 免费观看a级毛片全部| 国产亚洲av高清不卡| 精品免费久久久久久久清纯 | 高清视频免费观看一区二区| 欧美成人免费av一区二区三区 | 国产三级黄色录像| 久久国产精品人妻蜜桃| 啦啦啦免费观看视频1| 十分钟在线观看高清视频www| 中文字幕制服av| 一边摸一边抽搐一进一小说 | 欧美日韩乱码在线| 夫妻午夜视频| 午夜福利视频在线观看免费| 国产极品粉嫩免费观看在线| 999久久久国产精品视频| bbb黄色大片| 国产精品综合久久久久久久免费 | 亚洲成人国产一区在线观看| 男男h啪啪无遮挡| 51午夜福利影视在线观看| 亚洲成人免费电影在线观看| 久久这里只有精品19| 国产又爽黄色视频| 最新美女视频免费是黄的| 国产成人精品在线电影| 国产精品98久久久久久宅男小说| 人人妻人人澡人人看| 国产又色又爽无遮挡免费看| 精品少妇久久久久久888优播| 久久中文字幕一级| 女同久久另类99精品国产91| 久99久视频精品免费| 美女午夜性视频免费| 黄色a级毛片大全视频| 亚洲人成77777在线视频| 手机成人av网站| 国产熟女午夜一区二区三区| 国产在线观看jvid| 9色porny在线观看| 天堂中文最新版在线下载| 妹子高潮喷水视频| 欧美另类亚洲清纯唯美| 亚洲熟妇熟女久久| 热99国产精品久久久久久7| 欧美乱妇无乱码| 天天影视国产精品| 黑人欧美特级aaaaaa片| 免费女性裸体啪啪无遮挡网站| 亚洲三区欧美一区| 狂野欧美激情性xxxx| 亚洲成av片中文字幕在线观看| 欧美乱色亚洲激情| 中文字幕人妻熟女乱码| 国产高清激情床上av| 在线国产一区二区在线| 大型av网站在线播放| 精品一区二区三区视频在线观看免费 | 交换朋友夫妻互换小说| 亚洲全国av大片| 国产片内射在线| 99国产综合亚洲精品| 亚洲一区中文字幕在线| 久久性视频一级片| 日韩免费av在线播放| 免费在线观看亚洲国产| 亚洲,欧美精品.| 国产精品久久久久成人av| 免费在线观看黄色视频的| videos熟女内射| 精品国产乱子伦一区二区三区| 老熟妇仑乱视频hdxx| 亚洲精品国产一区二区精华液| 国产男女超爽视频在线观看| 亚洲 国产 在线| 欧美日韩视频精品一区| 在线观看免费高清a一片| 18禁观看日本| 成人免费观看视频高清| 国产精品99久久99久久久不卡| 男人的好看免费观看在线视频 | 黄色视频,在线免费观看| 中文字幕精品免费在线观看视频| 母亲3免费完整高清在线观看| 国产视频一区二区在线看| 欧美日韩乱码在线|