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

    基于多源數(shù)據(jù)的干旱區(qū)鹽漬化農(nóng)田精準(zhǔn)管理分區(qū)研究

    2022-04-15 01:26:14白建鐸王玉珍柳維揚(yáng)李洪義
    干旱區(qū)研究 2022年2期
    關(guān)鍵詞:鹽漬化表觀電導(dǎo)率

    白建鐸, 彭 杰, 史 舟, 王玉珍, 柳維揚(yáng), 李洪義

    (1.塔里木大學(xué)農(nóng)學(xué)院,新疆 阿拉爾 843300;2.浙江大學(xué)農(nóng)業(yè)遙感與信息技術(shù)應(yīng)用研究所,浙江 杭州310058;3.江西財經(jīng)大學(xué)旅游與城市管理學(xué)院,江西 南昌 330000)

    土壤鹽漬化不僅是全球面臨的生態(tài)問題[1],同時也是影響全球灌溉農(nóng)業(yè)最重要的因素之一[2],而且不同程度鹽漬化土壤對不同植物和農(nóng)作物生長發(fā)育影響也不盡相同。因此,掌握土壤鹽漬化信息可在一定程度上指導(dǎo)農(nóng)業(yè)合理灌溉和合理調(diào)整作物布局。近年來,隨氣候的變化,特別是全球氣候變暖將顯著影響全球水循環(huán),導(dǎo)致農(nóng)業(yè)用水量急劇減少[3]。在干旱區(qū)、半干旱區(qū),土壤鹽漬化問題和農(nóng)業(yè)灌溉缺水問題尤為突出。因此在干旱區(qū)、半干旱區(qū)急需一種根據(jù)鹽漬化程度進(jìn)行變量灌溉和調(diào)整作物布局的技術(shù),而精準(zhǔn)管理分區(qū)是解決此類問題的根本途徑,同時也是當(dāng)前精準(zhǔn)農(nóng)業(yè)研究的熱點(diǎn)之一[4]。

    精準(zhǔn)管理分區(qū)是根據(jù)土壤中各種制約產(chǎn)量的影響因子,劃分出具有相同均質(zhì)的子區(qū)域,來調(diào)整農(nóng)業(yè)投入品和作物布局的措施[5-6]。分區(qū)之間實(shí)施差異化管理措施有利于提高農(nóng)業(yè)生產(chǎn)投入效率、節(jié)約生產(chǎn)成本、保護(hù)生態(tài)環(huán)境,從而實(shí)現(xiàn)農(nóng)業(yè)的可持續(xù)發(fā)展。目前,在精準(zhǔn)管理分區(qū)方面,國內(nèi)外學(xué)者大多研究采用土壤養(yǎng)分?jǐn)?shù)據(jù)和產(chǎn)量相關(guān)數(shù)據(jù)進(jìn)行插值,再利用模糊聚類算法或自然斷點(diǎn)法進(jìn)行分區(qū)管理[7-9],但該方法在數(shù)據(jù)獲取方面需要耗費(fèi)大量的人力和物力,且不易大面積推廣。近年來,隨著遙感技術(shù)的發(fā)展,許多學(xué)者開始借助遙感影像數(shù)據(jù),再利用面向?qū)ο蠖喑叨确指钏惴ɑ蚰:垲惖人惴ㄟM(jìn)行分區(qū)管理[10-12],但此方法只能間接的通過光譜指數(shù)反映出土壤表層屬性數(shù)據(jù),而不能直接獲取到深層土壤屬性數(shù)據(jù)。同時也有部分學(xué)者使用土壤屬性數(shù)據(jù)與遙感影像數(shù)據(jù)相結(jié)合的多源數(shù)據(jù),利用面向?qū)ο蠖喑叨确指钏惴ㄟM(jìn)行分區(qū)管理[13],但該方法在土壤屬性數(shù)據(jù)的獲取方面也存在費(fèi)時又費(fèi)力的缺點(diǎn),只能在較小區(qū)域展開研究。上述研究均是針對土壤養(yǎng)分差異進(jìn)行管理分區(qū),基于土壤鹽漬化差異進(jìn)行管理分區(qū)的研究甚少。干旱區(qū)或半干旱區(qū),特別是新疆的南部地區(qū),土壤鹽漬化是影響作物產(chǎn)量和品質(zhì)的關(guān)鍵因子,需根據(jù)土壤鹽漬化程度采取不同的管理措施,如作物種植結(jié)構(gòu)的合理布局、冬灌或春灌的灌溉量定額。因此在這種土壤鹽漬化問題突出且鹽漬化變異強(qiáng)的南疆區(qū)域,應(yīng)以土壤鹽漬化差異作為分區(qū)關(guān)鍵因子進(jìn)行精準(zhǔn)管理分區(qū),才能提高農(nóng)業(yè)效益。傳統(tǒng)的土壤鹽漬化信息獲取方法主要依靠田間調(diào)查結(jié)合室內(nèi)分析,存在費(fèi)時、費(fèi)力、成本高等方面的不足,再考慮到鹽漬化強(qiáng)烈的時空變異特征,該方法很難滿足區(qū)域尺度土壤鹽漬化的調(diào)查,而以電磁感應(yīng)為原理的EM38-MK2大地電導(dǎo)率儀,可快速、高效地直接獲取不同深度土壤的表觀電導(dǎo)率數(shù)據(jù),土壤表觀電導(dǎo)率數(shù)據(jù)與土壤鹽分含量具有高度正相關(guān)性[14],可有效反映土壤鹽漬化程度,使用EM38-MK2大地電導(dǎo)率儀可顯著提高土壤鹽漬化信息的獲取效率,特別是針對區(qū)域尺度的調(diào)查而言。自20世紀(jì)末以來,使用表觀電導(dǎo)率來表征土壤鹽漬化的方法在土壤鹽漬化監(jiān)測和評價中已得到了廣泛應(yīng)用[15-18],但直接應(yīng)用表觀電導(dǎo)率數(shù)據(jù)進(jìn)行農(nóng)田管理分區(qū)的報道鮮見。此外,目前關(guān)于農(nóng)田管理分區(qū)的研究普遍只采用了土壤的基本信息,缺乏綜合利用土壤信息和作物信息進(jìn)行管理分區(qū)的報道,而衛(wèi)星遙感技術(shù)可以快速獲取大面積作物的動態(tài)信息,為這種研究思路的實(shí)施提供了有力條件。因此可使用土壤表觀電導(dǎo)率和遙感影像數(shù)據(jù)進(jìn)行土壤鹽漬化的分區(qū)管理。

    綜上所述,本研究以南疆干旱區(qū)典型綠洲阿拉爾墾區(qū)農(nóng)田為研究對象,以研究區(qū)不同深度土壤的表觀電導(dǎo)率數(shù)據(jù)、地形數(shù)據(jù)和衛(wèi)星遙感影像數(shù)據(jù)為數(shù)據(jù)源,以土壤表層的表觀電導(dǎo)率數(shù)據(jù)為主變量,在分析土壤表觀電導(dǎo)率的空間異質(zhì)性和空間自相關(guān)性的基礎(chǔ)上,通過相關(guān)分析篩選出與表觀電導(dǎo)達(dá)到極顯著相關(guān)的衛(wèi)星遙感光譜指數(shù),并將其作為輔助變量,再利用面向?qū)ο蠖喑叨确指钏惴ㄟM(jìn)行不同輸入量下土壤鹽漬化的分區(qū),并對不同輸入量下的管理分區(qū)結(jié)果進(jìn)行評價對比,找到最優(yōu)管理分區(qū),并提取出研究區(qū)各田塊的邊界,從而根據(jù)田塊所在的鹽漬化分區(qū)指導(dǎo)農(nóng)田的精準(zhǔn)灌溉和作物的合理布局,以實(shí)現(xiàn)作物優(yōu)質(zhì)高產(chǎn)的目的。

    1 研究區(qū)概況

    研究區(qū)位于新疆維吾爾族自治區(qū)一師阿拉爾市十二團(tuán)境內(nèi)(圖1),地理坐標(biāo)為81°19′31″E,40°29′52″N,面積為1144 hm2,地處塔里木盆地西北部、塔克拉瑪干沙漠北緣,天山南麓。境內(nèi)主要有葉爾羌河、和田河、阿克蘇河和塔里木河四大水系。冬季嚴(yán)寒,夏季炎熱,日照時間長,降水稀少,蒸發(fā)強(qiáng)烈。該區(qū)多年平均降水量為48.5 mm,多年平均蒸發(fā)量為1988 mm,蒸降比高達(dá)40.99,系典型的暖溫帶極端大陸性干旱荒漠氣候。研究區(qū)地勢平坦,土壤質(zhì)地以砂壤土為主,區(qū)域內(nèi)農(nóng)田主要種植棉花,棉花生育期內(nèi)均采用膜下滴灌方式進(jìn)行灌溉,而春灌和冬灌則采用大水漫灌的方式。大水漫灌旨在將整個棉花生育期內(nèi)運(yùn)移到土壤表層的鹽分淋洗到耕層以下,以此降低下茬作物種植的鹽害脅迫,其灌溉用水主要來自塔里木河。研究區(qū)土壤鹽漬化程度差異較大,種植結(jié)構(gòu)單一,冬春灌采用統(tǒng)一灌溉量進(jìn)行灌溉,不僅易造成局部田塊壓鹽不徹底導(dǎo)致鹽害現(xiàn)象產(chǎn)生或過度淋洗浪費(fèi)水資源,同時也造成研究區(qū)作物布局不合理導(dǎo)致作物經(jīng)濟(jì)效益下降,因此,研究區(qū)迫切需要根據(jù)不同鹽漬化程度來進(jìn)行變量灌溉和種植結(jié)構(gòu)的調(diào)整。

    圖1 研究區(qū)位置及樣點(diǎn)分布示意圖Fig.1 Location of the study area and distribution of sample points

    2 數(shù)據(jù)的獲取與研究方法

    2.1 多源數(shù)據(jù)的獲取

    2.1.1 電磁感應(yīng)數(shù)據(jù)的獲取 土壤表觀電導(dǎo)率數(shù)據(jù)采集儀器為大地電導(dǎo)率儀(EM38-MK2),在棉花收獲后冬灌前的11 月上旬進(jìn)行土壤表觀電導(dǎo)率數(shù)據(jù)的采集。每次采集表觀電導(dǎo)率時,首先將大地電導(dǎo)率(EM38-MK2)設(shè)置為手動測定模式,并對儀器進(jìn)行預(yù)熱和歸零校驗,待儀器校驗完畢后,進(jìn)行表觀電導(dǎo)率的獲取;在每個樣點(diǎn)采集時,儀器將自動記錄每個樣點(diǎn)的地理坐標(biāo),每個采樣點(diǎn)同時采集垂直模式和水平模式兩套數(shù)據(jù);水平模式采集的為0~0.375 m 和0~0.75 m 的2 個深度土層的表觀電導(dǎo)率數(shù)據(jù),即ECh0.375和ECh0.75,垂直模式則為0~0.75 m 和0~1.5 m 的兩個深度土層表觀電導(dǎo)率數(shù)據(jù),即ECv0.75和ECv1.5,相鄰采樣點(diǎn)間距小于300 m。研究區(qū)共隨機(jī)均勻采集約1200 組表觀電導(dǎo)率數(shù)據(jù)(圖1),采集數(shù)據(jù)描述性統(tǒng)計見表1。

    表1 不同深度表觀電導(dǎo)率描述Tab.1 Descriptive statistics of apparent electrical conductivity at different depths

    2.2.2 遙感影像數(shù)據(jù)和地形數(shù)據(jù)的獲取 研究采用來源于美國地質(zhì)勘查局(USGS, http://glovis.usgs.gov)的Landsat8遙感影像,衛(wèi)星重訪周期為16 d,空間分辨率為30 m。本文選取2020 年棉花生長旺盛期(9 月5 日)的影像和棉花收獲冬灌前(11 月4 日)的裸土影像。影像經(jīng)輻射定標(biāo)、大氣校正、圖像裁剪等預(yù)處理后,使用ENVI 5.3 軟件計算得出9 月5日影像的9個植被指數(shù)(RVI、DVI、GDVI、GRVI、NDVI、EVI、NNIR、NR、NG)和11 月4 日影像的13 個土壤鹽分指數(shù)(NDSI、S1、S2、S3、S5、S6、SI-T、SI、SI1、SI2、SI3、Int1、Int2)。同時通過地理空間數(shù)據(jù)云(http://www.gscloud.cn)下載空間分辨率為30 m 的DEM 數(shù)據(jù),使用ArcGIS 10.7 軟件完成洼地填充、影像裁剪、幾何精度校正和高程數(shù)據(jù)提取。

    2.2 分割方法

    2.2.1 基于面向?qū)ο蟮亩喑叨确指?面向?qū)ο蠖喑叨确指钍怯跋駥ο蟮囊粋€兩兩合并的過程,合并從均勻分布在影像上的影像對象開始[19]。對于每個影像對象,找到一個使它們之間的異質(zhì)性變化最小相鄰影像對象,在異質(zhì)性變化最小處,2個對象合并成1 個更大的對象,每個對象在每個循環(huán)中處理1次,直到分割對象的異質(zhì)性超過用戶給定的閾值,合并過程停止[20]。本研究嘗試基于多源數(shù)據(jù)的土壤鹽漬化空間變異進(jìn)行分區(qū),在eCognition 軟件中首先確定合適的形狀因子(Shape)和緊致度因子(Compactness),然后通過調(diào)節(jié)不同的分割尺度參數(shù)(Scale parameter),使分區(qū)之間表觀電導(dǎo)率異質(zhì)性最大,分區(qū)內(nèi)部電導(dǎo)率同質(zhì)性最大。

    2.2.2 確定最優(yōu)管理分區(qū) 由于面向?qū)ο蠓指畹脑瓌t為盡量使分割對象之間的異質(zhì)性高,同時使分割對象內(nèi)部的同質(zhì)性高[21],所以為提高分區(qū)內(nèi)部的同質(zhì)性和分區(qū)之間的異質(zhì)性需要尋找最優(yōu)分區(qū)尺度;根據(jù)學(xué)者研究,可使用平均分割評價指數(shù)的變化來確定最優(yōu)分割尺度,將平均分割評價指數(shù)和尺度參數(shù)繪制成曲線圖,在平均分割評價指數(shù)開始產(chǎn)生明顯變化時的前一個尺度參數(shù)為最優(yōu)分割尺度參數(shù)[11]。通過在多尺度分割算法中調(diào)節(jié)分割尺度參數(shù)將研究區(qū)分為m個區(qū)域,并計算出不同分割尺度下各區(qū)域的標(biāo)準(zhǔn)差和各區(qū)域與鄰域的平均差分絕對值,其同質(zhì)性指數(shù)用標(biāo)準(zhǔn)差(δ)來表示,異質(zhì)性指數(shù)用平均差分絕對值(ΔCL)來表示,計算如公式(1)、(2)所示;同時構(gòu)建對象的分割評價指數(shù)(SEI)和平均分割評價指標(biāo)(ASEI),計算如公式(3)和(4)所示[22]。

    式中:n為某區(qū)域內(nèi)像元的個數(shù);CLi為該區(qū)域內(nèi)第i個像元的屬性值;CL為該區(qū)域的屬性平均值;δ表示該區(qū)域的標(biāo)準(zhǔn)差;L為該區(qū)域的周長;N表示與該區(qū)域相鄰區(qū)域的個數(shù);Li表示與第i個相鄰區(qū)域公共邊的長度;CLp為第p個相鄰區(qū)域的屬性平均值;ΔCL表示該區(qū)域與鄰域的平均差分絕對值;SEI表示分割評價指數(shù);A表示研究區(qū)的總面積,Ai為i個區(qū)域的面積;m為被劃分的區(qū)域個數(shù);SEIi為第i個區(qū)域的分割評價指數(shù);ASEI表示平均分割評價指數(shù)。

    2.3 精準(zhǔn)管理分區(qū)的評價方法

    2.3.1 分區(qū)內(nèi)部的同質(zhì)性 根據(jù)分區(qū)內(nèi)像元屬性值的標(biāo)準(zhǔn)差和平均值,分別計算出各個分區(qū)內(nèi)的變異系數(shù),并考慮到各分區(qū)面積大小不同而賦予相應(yīng)的權(quán)重,從而用分區(qū)內(nèi)部變異系數(shù)(CV)的平均值表示同質(zhì)性,計算如公式(5)所示:

    式中:CVi為第i個分區(qū)的變異系數(shù);ai為第i個分區(qū)的面積;m為被劃分的區(qū)域個數(shù)。CV越小,分區(qū)內(nèi)變異越小,代表分區(qū)內(nèi)異質(zhì)性較小,同質(zhì)性較大。

    2.3.2 分區(qū)之間的異質(zhì)性 采用全局莫蘭指數(shù)(I)判斷空間相關(guān)性,根據(jù)莫蘭指數(shù)大小判斷分區(qū)間異質(zhì)性高低,其計算如公式(6)所示:

    式中:m為被劃分的區(qū)域個數(shù);yi為分區(qū)i的屬性平均值;yˉ為整個研究區(qū)的屬性平均值;Wij為i和j之間的空間權(quán)重,若i與j相鄰,則Wij=1,否則Wij=0。I表示空間相關(guān)程度,I越接近于0,表示分區(qū)之間相關(guān)程度越低,即分區(qū)之間異質(zhì)性越高,而根據(jù)分區(qū)評價原則可知,分區(qū)之間異質(zhì)性越高,分區(qū)結(jié)果越好。

    3 結(jié)果與分析

    3.1 土壤鹽漬化的空間異質(zhì)性

    根據(jù)土壤鹽漬化來劃分管理分區(qū),需證明該研究區(qū)土壤表觀電導(dǎo)率具有空間異質(zhì)性,因此需使用地統(tǒng)計學(xué)中半變異函數(shù)進(jìn)行分析。在計算半變異函數(shù)時,數(shù)據(jù)需符合正態(tài)分布,否則會出現(xiàn)比例效應(yīng)。經(jīng)分析,原始數(shù)據(jù)呈偏態(tài)分布,經(jīng)平方根轉(zhuǎn)化后基本符合正態(tài)分布。通過GS+軟件獲得變異函數(shù)的最優(yōu)模型,其參數(shù)如表2所示。

    表2 不同深度表觀電導(dǎo)率空間異質(zhì)性Tab.2 Spatial heterogeneity analysis of apparent electrical conductivity at different depths

    由表2 可知,各個土層表觀電導(dǎo)率的最優(yōu)模型均為指數(shù)模型,決定系數(shù)均大于0.93,且殘差均小于1.27,說明該模型可較好地反映研究區(qū)土壤鹽漬化的空間分布情況。塊金值表示人為隨機(jī)因素造成的空間異質(zhì)性,各土層表觀電導(dǎo)率的塊金值在3.05~5.33 之間,表明研究區(qū)表觀電導(dǎo)率存在由人為隨機(jī)因素所造成的空間異質(zhì)性?;_值表示人為隨機(jī)因素和自然結(jié)構(gòu)因素造成的最大程度空間異質(zhì)性,各土層表觀電導(dǎo)率的基臺值介于6.27~10.79 之間,表明各土層的表觀電導(dǎo)率均存在著顯著的空間變異。塊金值與基臺值之比表示隨機(jī)因素引起的空間異質(zhì)性占總空間異質(zhì)性的大小,各土層表觀電導(dǎo)率的塊金值與基臺值之比均接近于50%,表明隨機(jī)性因素和結(jié)構(gòu)性因素共同影響研究區(qū)土壤表觀電導(dǎo)率的空間異質(zhì)性。變程表示空間自相關(guān)的最大距離,各土層表觀電導(dǎo)率的變程介于1383~2579 m之間,并且采集數(shù)據(jù)時,各樣點(diǎn)距離控制在300 m之內(nèi),均小于變程,表明合理的采樣距離不僅可以提高插值的精度,還可以準(zhǔn)確地揭示研究區(qū)土壤鹽漬化的空間異質(zhì)性。因此可根據(jù)土壤鹽漬化空間異質(zhì)性劃分管理分區(qū)。

    3.2 表觀電導(dǎo)率與光譜指數(shù)和地形數(shù)據(jù)的相關(guān)性分析

    為探究表層土壤表觀電導(dǎo)率(ECh0.375)與地形數(shù)據(jù)(DEM)及經(jīng)不同時期遙感影像提取所得光譜指數(shù)之間的相關(guān)性,將DEM 數(shù)據(jù)及所提取的22 個光譜數(shù)據(jù)與表層土壤表觀電導(dǎo)率數(shù)據(jù)在Excel 中進(jìn)行相關(guān)系數(shù)計算,結(jié)果如表3 所示。通過分析棉花生長旺盛期9 月的9 個植被指數(shù)可知,DVI 和GDVI 與土壤表觀電導(dǎo)率未達(dá)到顯著相關(guān),而RVI、GRVI 和NDVI 等7 個植被指數(shù)與表觀電導(dǎo)率達(dá)到極顯著相關(guān)。對棉花收獲后冬灌前11月裸土的13個土壤鹽分指數(shù)分析可知,NDSI、S3 和S5 等8 個土壤鹽分指數(shù)與土壤表觀電導(dǎo)率達(dá)到顯著相關(guān),其中NDSI、S5和SI-T 與土壤表觀電導(dǎo)率達(dá)到極顯著正相關(guān)。此外,由表3 可知,植被指數(shù)中除NG 外其余均與表觀電導(dǎo)率呈負(fù)相關(guān),而土壤鹽分指數(shù)卻呈相反結(jié)果,除S6外其余均與表觀電導(dǎo)率呈正相關(guān);這與彭杰[23]對多光譜波段和光譜指數(shù)與電導(dǎo)率的相關(guān)性關(guān)系研究結(jié)論相一致。從地形數(shù)據(jù)分析可知,DEM數(shù)據(jù)與表觀電導(dǎo)率未達(dá)到顯著性水平,這主要原因包括兩方面,一是研究區(qū)地勢平坦,削弱了地形條件對鹽分空間分布的支配作用,二是灌溉、土壤改良等人類管理活動改變了地形條件控制下鹽分空間分布的原始格局。本研究在相關(guān)性分析基礎(chǔ)上篩選出相關(guān)性達(dá)到極顯著,且相關(guān)系數(shù)r≥0.10的植被指數(shù)和土壤鹽分指數(shù),將其作為精準(zhǔn)管理分區(qū)的輔助變量,可有效提高分區(qū)內(nèi)同質(zhì)性和分區(qū)間異質(zhì)性。

    表3 表觀電導(dǎo)率與光譜指數(shù)和地形數(shù)據(jù)的相關(guān)性分析Tab.3 Correlation analysis table of apparent electrical conductivity with spectral index and DEM data

    3.3 基于表層表觀電導(dǎo)率與不同變量下的精準(zhǔn)管理分區(qū)

    基于各土層表觀電導(dǎo)率空間結(jié)構(gòu)的分析,為了更加直觀地揭示研究區(qū)土壤鹽分的空間異質(zhì)性特征,使用普通克里格插值法繪制出各土層的表觀電導(dǎo)率分布圖。因各變量間單位不同,需將篩選出的植被指數(shù)(RVI、GRVI、EVI)、土壤鹽分指數(shù)(NDSI、S5、SI-T)及各土層土壤表觀電導(dǎo)率插值圖進(jìn)行標(biāo)準(zhǔn)化處理,將變量控制在0~1 之間。然后將各土層表觀電導(dǎo)率插值圖重采樣為30 m×30 m 的像元,保持與其他變量像元大小一致。在棉花收獲后,土壤中大部分的鹽分在生育期內(nèi)會通過蒸發(fā)作用將深層鹽分隨水運(yùn)移到表層土壤,因此本文選取表層表觀電導(dǎo)率(ECh0.375)為主變量,其他變量為輔助變量進(jìn)行管理分區(qū)。同時,通過eCognition 軟件提取出農(nóng)田邊界,進(jìn)而進(jìn)行對農(nóng)田精準(zhǔn)管理分區(qū)。

    通過多次多尺度分割試驗確定最為合適的形狀因子為0.3,緊致度因子為0.5。在控制形狀因子與緊致度因子不變的條件下,通過最優(yōu)分區(qū)評價指數(shù)曲線(圖2)可知,基于ECh0.375的最優(yōu)分割尺度為19,分區(qū)數(shù)為19;基于ECh0.375和深層表觀電導(dǎo)率的最優(yōu)分割尺度為19,分區(qū)數(shù)為16;基于ECh0.375和植被指數(shù)的最優(yōu)分割尺度為14,分區(qū)數(shù)為17;基于ECh0.375和土壤鹽分指數(shù)的最優(yōu)分割尺度為14,分區(qū)數(shù)為13。不同變量下的管理分區(qū)如圖3所示。由圖3可知,不同變量下的分區(qū)結(jié)果,均可將不同程度電導(dǎo)率區(qū)域進(jìn)行分區(qū),表明基于多源數(shù)據(jù)的面向?qū)ο蠖喑叨确指畹墓芾矸謪^(qū)具有可行性。

    圖2 精準(zhǔn)管理分區(qū)評價指數(shù)Fig.2 Segmentation evaluation index

    3.4 不同輸入量下管理分區(qū)對比與評價

    由圖3對比基于不同數(shù)據(jù)下的4種最優(yōu)管理分區(qū)發(fā)現(xiàn),由于輔助變量的差異,分區(qū)結(jié)果也存在明顯差異。從農(nóng)田的耕作角度出發(fā),基于ECh0.375的分區(qū)結(jié)果,存在相對窄而長的分區(qū),這種分區(qū)與實(shí)際機(jī)械化工作不相符;當(dāng)引入輔助變量后,分區(qū)結(jié)果在每個分區(qū)中最少有2~3 個田塊,與實(shí)際機(jī)械化工作條件相符合。從分割效果出發(fā),基于ECh0.375的分割尺度與基于ECh0.375和深層表觀電導(dǎo)率的分割尺度相同,但分區(qū)個數(shù)有差異,這可能是由深層土壤表觀電導(dǎo)率與表層表觀電導(dǎo)率分布格局不同引起的;基于ECh0.375和植被指數(shù)的分割尺度與基于ECh0.375和土壤指數(shù)的分割尺度相同,但分區(qū)個數(shù)不同。上述對比只是對不同分區(qū)結(jié)果進(jìn)行了定性對比,而未進(jìn)行定量評價。分區(qū)結(jié)果的定量評價,需本著分區(qū)內(nèi)同質(zhì)性最大,且分區(qū)間異質(zhì)性最大原則展開評價,即分區(qū)內(nèi)平均變異系數(shù)越小,分區(qū)間莫蘭指數(shù)數(shù)值越接近于零,代表分區(qū)效果越好。由表4可知,不同數(shù)據(jù)下各分區(qū)內(nèi)平均變異系數(shù)在13.60%~18.65%之間,根據(jù)劉寧[24]對不同土地利用方式下黃河三角洲土壤特性空間變異研究可知,變異系數(shù)在0~20%之間屬于弱變異,在20%~50%屬于中等變異,>50%屬于強(qiáng)變異;說明各分區(qū)內(nèi)均屬于弱變異;且較整個研究區(qū)的變異系數(shù)下降了約60%,說明分區(qū)內(nèi)同質(zhì)性較高。不同數(shù)據(jù)下各分區(qū)間莫蘭指數(shù)均接近于0,表示各分區(qū)間相關(guān)性較低,說明各分區(qū)間異質(zhì)性較高。綜上從農(nóng)田耕作角度、分割效果及分區(qū)評價原則出發(fā),基于ECh0.375和深層表觀電導(dǎo)率的管理分區(qū)為最佳,因為該分區(qū)結(jié)果既符合當(dāng)?shù)剞r(nóng)田的管理和機(jī)械化作業(yè)要求,又可以代表深層鹽分的分區(qū),且分區(qū)間異質(zhì)性也最強(qiáng)?;贓Ch0.375和土壤鹽分指數(shù)的管理分區(qū)次之?;贓Ch0.375的管理分區(qū)最差,是由于其分區(qū)結(jié)果存在著窄而長的分區(qū),與實(shí)際耕作和機(jī)械化作業(yè)不相符。

    圖3 不同變量下的最優(yōu)管理分區(qū)Fig.3 Optimal management zones under different variables

    表4 基于不同數(shù)據(jù)集下的分區(qū)評價Tab.4 Partition evaluation based on different dataset

    4 討論

    新疆南疆地區(qū)是我國棉花主要生產(chǎn)基地,但該地區(qū)鹽漬化造成的生態(tài)環(huán)境惡化、作物生產(chǎn)力下降、土壤肥力降低問題,是制約該地區(qū)作物高產(chǎn)的主要障礙因子,因此對該地區(qū)鹽漬化分區(qū)具有重要意義。本研究選取了阿拉爾墾區(qū)棉花收獲后冬灌前的農(nóng)田為研究對象,此時期鹽分聚集在土壤表層[25],且研究區(qū)土壤質(zhì)地以砂壤土為主,保水保肥能力差,土壤中水分,養(yǎng)分不穩(wěn)定,受環(huán)境影響較大,因此在研究區(qū)進(jìn)行鹽漬化分區(qū)比進(jìn)行水分,養(yǎng)分分區(qū)更有意義。在研究方法上,基于面向?qū)ο蟮亩喑叨确指钏惴ㄊ悄壳斑M(jìn)行農(nóng)田管理分區(qū)研究中最常用的方法之一,其分區(qū)數(shù)是由最優(yōu)分割尺度所決定,不僅消除了人為因素的影響,還能夠精準(zhǔn)地劃分管理分區(qū),但最優(yōu)分割尺度同時受形狀因子、緊致度因子、影像空間分辨率、研究區(qū)面積大小等多個因素影響。本文通過與劉煥軍等[11]利用遙感影像對黑土區(qū)田塊尺度的精準(zhǔn)管理分區(qū)對比發(fā)現(xiàn),影像分辨率不同其最優(yōu)分割尺度也不同,且最優(yōu)分割尺度隨著空間分辨率的增加而減小,即空間分辨率越高,分割尺度越小,這與許盼盼[26]基于高分辨率遙感數(shù)據(jù)的濕地精細(xì)信息提取相一致。但分辨率越高,分割結(jié)果可能較為破碎,其分區(qū)結(jié)果可能不適合新疆這種高度機(jī)械化的農(nóng)事作業(yè)的開展。同時最優(yōu)分割尺度與研究區(qū)面積大小也具有一定關(guān)系,根據(jù)黃萬里[27]對不同面積的森林對象進(jìn)行多尺度分割發(fā)現(xiàn)隨著研究區(qū)面積的增大,最優(yōu)分割尺度也隨之變大。因此在使用面向?qū)ο蠖喑叨确指顣r不僅要考慮研究區(qū)所占面積大小,更要選擇合適的分辨率,從而確定最優(yōu)分割尺度。經(jīng)對比發(fā)現(xiàn)多源數(shù)據(jù)較單一數(shù)據(jù)可較好地劃分管理分區(qū),且在本質(zhì)上更具有說服力[21]。本文基于輔助變量的多源數(shù)據(jù)管理分區(qū)消除了過于破碎化的分區(qū),且分區(qū)間的異質(zhì)性也有明顯提高,但可能存在數(shù)據(jù)冗余和計算量增加等問題,因此在后期研究中需改善此類問題。在分區(qū)管理方面,接下來將根據(jù)鹽漬化程度計算出各個分區(qū)壓鹽所需要的灌溉量,通過測量各個分區(qū)的土壤自然含水量、土壤容重,田間持水量等基礎(chǔ)參數(shù),采用列果斯塔也夫經(jīng)驗公式,進(jìn)行各個分區(qū)灌溉量的計算[23],從而達(dá)到精準(zhǔn)灌溉的目標(biāo)。

    5 結(jié)論

    本研究運(yùn)用以電磁感應(yīng)原理為基礎(chǔ)的EM38-MK2大地電導(dǎo)率儀,獲取棉花收獲后冬灌前表觀電導(dǎo)率數(shù)據(jù)。用半變異函數(shù)進(jìn)行分析,揭示土壤鹽漬化的空間異質(zhì)性。利用空間分辨率為30 m的Landsat8影像不同波段信息,計算出22個光譜指數(shù)。將光譜指數(shù)和地形數(shù)據(jù)分別與表觀電導(dǎo)率(ECh0.375)展開相關(guān)性分析,篩選得出3 個植被指數(shù)和3 個土壤鹽分指數(shù),并將其作為精準(zhǔn)管理分區(qū)的輔助變量,最后進(jìn)行多源數(shù)據(jù)的土壤鹽漬化管理分區(qū)研究,所得結(jié)論如下:

    (1)在半變異函數(shù)分析中,各土層表觀電導(dǎo)率的基臺值與塊金值之比均在50%左右,說明研究區(qū)土壤鹽漬化存在空間變異,且是由隨機(jī)性因素和結(jié)構(gòu)性因素共同導(dǎo)致的變異,因而可以根據(jù)研究區(qū)土壤鹽漬化程度進(jìn)行分區(qū)。

    (2)經(jīng)相關(guān)分析和指數(shù)篩選后,9 月植被指數(shù)(RVI、GRVI、EVI)和11 月裸地的土壤鹽分指數(shù)(NDSI、S5、SI-T)均與ECh0.375達(dá)到極顯著水平,因此將其作為土壤鹽漬化管理分區(qū)的輔助變量。

    (3)分區(qū)后結(jié)果表明,各分區(qū)結(jié)果的變異系數(shù)均小于19%,且較整個研究區(qū)的變異系數(shù)下降約60%,說明面向?qū)ο蠖喑叨确指钏惴捎行У貏澐殖鼍|(zhì)性區(qū)域,從而指導(dǎo)農(nóng)田的差異化管理。

    (4)對比分區(qū)結(jié)果發(fā)現(xiàn),基于單一數(shù)據(jù)的管理分區(qū)其分區(qū)間異質(zhì)性較低,且部分分區(qū)不符合耕作單元,而基于多源數(shù)據(jù)的管理分區(qū),不僅符合耕作單元,且分區(qū)間異質(zhì)性較單一數(shù)據(jù)源的管理分區(qū)均有所提高。

    (5)從農(nóng)田的耕作角度、分割效果及分區(qū)評價原則出發(fā)可得,基于ECh0.375和深層表觀電導(dǎo)率的管理分區(qū)為最佳,該分區(qū)結(jié)果既符合當(dāng)?shù)剞r(nóng)田管理又符合機(jī)械化作業(yè);而基于ECh0.375和土壤鹽分指數(shù)的管理分區(qū)次之,效果最差的為僅基于ECh0.375單一數(shù)據(jù)源的管理分區(qū)。

    猜你喜歡
    鹽漬化表觀電導(dǎo)率
    蔬菜大棚土壤鹽漬化成因及防治措施
    土地質(zhì)量地球化學(xué)調(diào)查成果在判定土壤鹽漬化、沙化中的應(yīng)用
    綠盲蝽為害與赤霞珠葡萄防御互作中的表觀響應(yīng)
    河北果樹(2021年4期)2021-12-02 01:14:50
    甘肅蘇干湖濕地土壤鹽漬化、地下水位埋深及其對生態(tài)環(huán)境的影響
    鋼結(jié)構(gòu)表觀裂紋監(jiān)測技術(shù)對比與展望
    上海公路(2019年3期)2019-11-25 07:39:28
    例析對高中表觀遺傳學(xué)的認(rèn)識
    基于比較測量法的冷卻循環(huán)水系統(tǒng)電導(dǎo)率檢測儀研究
    低溫脅迫葡萄新梢電導(dǎo)率和LT50值的研究
    瑪納斯河流域土壤鹽漬化時空動態(tài)變化
    高電導(dǎo)率改性聚苯胺的合成新工藝
    亚洲精品一卡2卡三卡4卡5卡| 一级毛片久久久久久久久女| 国产精华一区二区三区| 特大巨黑吊av在线直播| 日日摸夜夜添夜夜添小说| 成人特级黄色片久久久久久久| 国产一区二区在线观看日韩| 中国国产av一级| 国产精品一区www在线观看| 国产在视频线在精品| 白带黄色成豆腐渣| 久久精品国产亚洲av香蕉五月| 精品免费久久久久久久清纯| 如何舔出高潮| 免费看光身美女| 在线观看av片永久免费下载| 亚洲一级一片aⅴ在线观看| 国产精品久久久久久av不卡| 成人亚洲精品av一区二区| 国产午夜精品论理片| 一区二区三区高清视频在线| 亚洲国产欧美人成| 性色avwww在线观看| 国产久久久一区二区三区| 久久久久久久久久成人| 免费在线观看影片大全网站| 深夜精品福利| 好男人在线观看高清免费视频| 久久久久久久久久久丰满| 日韩精品中文字幕看吧| 久久99热这里只有精品18| 亚洲aⅴ乱码一区二区在线播放| 欧美中文日本在线观看视频| 黄色配什么色好看| www日本黄色视频网| 超碰av人人做人人爽久久| 99热全是精品| 欧美人与善性xxx| 免费看a级黄色片| 最近视频中文字幕2019在线8| 国产黄色视频一区二区在线观看 | 最后的刺客免费高清国语| 蜜桃久久精品国产亚洲av| 中文字幕av在线有码专区| 国产精品综合久久久久久久免费| 久久久精品94久久精品| 色视频www国产| 国产精品人妻久久久久久| 久久精品人妻少妇| a级毛片免费高清观看在线播放| 亚洲国产精品成人久久小说 | 寂寞人妻少妇视频99o| 久久久久久久久久黄片| 成人av在线播放网站| 日韩大尺度精品在线看网址| 嫩草影视91久久| a级毛片免费高清观看在线播放| 可以在线观看的亚洲视频| 国产av麻豆久久久久久久| 91久久精品电影网| 久久综合国产亚洲精品| 亚洲成a人片在线一区二区| 日韩成人伦理影院| 男人舔女人下体高潮全视频| 最近手机中文字幕大全| 99精品在免费线老司机午夜| 国产男靠女视频免费网站| 俄罗斯特黄特色一大片| 欧美极品一区二区三区四区| 精品福利观看| 免费电影在线观看免费观看| 国产真实乱freesex| 99热精品在线国产| 一区二区三区高清视频在线| 一级毛片久久久久久久久女| 晚上一个人看的免费电影| 大香蕉久久网| 免费不卡的大黄色大毛片视频在线观看 | 美女被艹到高潮喷水动态| 欧美一区二区精品小视频在线| 变态另类丝袜制服| 久久99热6这里只有精品| 国产午夜精品论理片| 日韩欧美精品v在线| 亚洲av一区综合| 91精品国产九色| 婷婷亚洲欧美| 99热只有精品国产| 午夜激情福利司机影院| 尾随美女入室| 看非洲黑人一级黄片| 一个人观看的视频www高清免费观看| 欧美另类亚洲清纯唯美| 一级黄色大片毛片| 国产精品人妻久久久影院| 精品一区二区三区视频在线| 日韩一区二区视频免费看| 国产三级中文精品| 在线观看美女被高潮喷水网站| 蜜桃亚洲精品一区二区三区| 高清午夜精品一区二区三区 | 好男人在线观看高清免费视频| 国产欧美日韩精品一区二区| 精品人妻视频免费看| 亚洲av一区综合| 国产91av在线免费观看| 日本成人三级电影网站| 麻豆国产97在线/欧美| 少妇熟女欧美另类| 中文在线观看免费www的网站| 国产成人a区在线观看| 国产成人福利小说| av在线天堂中文字幕| 男人的好看免费观看在线视频| 别揉我奶头~嗯~啊~动态视频| 能在线免费观看的黄片| 免费看av在线观看网站| 精品久久久久久久久久久久久| 亚洲丝袜综合中文字幕| 日本熟妇午夜| 久久精品91蜜桃| 2021天堂中文幕一二区在线观| 亚洲综合色惰| 色哟哟·www| 两个人视频免费观看高清| 午夜爱爱视频在线播放| 极品教师在线视频| 免费一级毛片在线播放高清视频| 熟妇人妻久久中文字幕3abv| 国内精品美女久久久久久| 国产一区二区三区av在线 | 一a级毛片在线观看| 亚洲婷婷狠狠爱综合网| 成人综合一区亚洲| av国产免费在线观看| 亚洲五月天丁香| 精品久久久久久久久久免费视频| 午夜影院日韩av| 国产成人a∨麻豆精品| 国产综合懂色| 日本免费一区二区三区高清不卡| 狂野欧美激情性xxxx在线观看| 蜜桃亚洲精品一区二区三区| 小说图片视频综合网站| 亚洲精品粉嫩美女一区| 夜夜夜夜夜久久久久| 老熟妇乱子伦视频在线观看| 精品久久国产蜜桃| 干丝袜人妻中文字幕| 精品福利观看| 十八禁网站免费在线| 亚洲欧美成人综合另类久久久 | 欧美一区二区亚洲| 观看免费一级毛片| 亚洲乱码一区二区免费版| 亚洲国产精品成人综合色| 一本精品99久久精品77| 国产三级在线视频| 成人无遮挡网站| 99热精品在线国产| 男女做爰动态图高潮gif福利片| 国产精华一区二区三区| 国产在线男女| 丰满的人妻完整版| 中文字幕精品亚洲无线码一区| 免费看光身美女| 国产私拍福利视频在线观看| 欧美丝袜亚洲另类| 日本黄色视频三级网站网址| 男女之事视频高清在线观看| 国产大屁股一区二区在线视频| 精品人妻偷拍中文字幕| 18禁在线无遮挡免费观看视频 | 两个人视频免费观看高清| 国产色爽女视频免费观看| 国产精品久久久久久亚洲av鲁大| 一本久久中文字幕| 伊人久久精品亚洲午夜| 床上黄色一级片| 免费看日本二区| 在线天堂最新版资源| a级毛片a级免费在线| 欧美中文日本在线观看视频| 久久亚洲国产成人精品v| 在线观看免费视频日本深夜| eeuss影院久久| АⅤ资源中文在线天堂| 亚洲欧美日韩卡通动漫| 国产精品伦人一区二区| 赤兔流量卡办理| 女同久久另类99精品国产91| 国产一区二区三区av在线 | 日本在线视频免费播放| 村上凉子中文字幕在线| 精品一区二区三区视频在线| av福利片在线观看| 偷拍熟女少妇极品色| 亚洲天堂国产精品一区在线| 五月伊人婷婷丁香| 一本一本综合久久| 网址你懂的国产日韩在线| 哪里可以看免费的av片| 亚洲不卡免费看| 悠悠久久av| 狠狠狠狠99中文字幕| 亚洲在线观看片| 99国产极品粉嫩在线观看| 成人无遮挡网站| 国产亚洲精品综合一区在线观看| 亚洲av免费在线观看| 桃色一区二区三区在线观看| 两个人的视频大全免费| 美女内射精品一级片tv| 国产探花极品一区二区| 国产一区二区亚洲精品在线观看| 国产91av在线免费观看| 村上凉子中文字幕在线| 亚洲四区av| 欧美激情国产日韩精品一区| 亚洲熟妇中文字幕五十中出| 成熟少妇高潮喷水视频| 人人妻,人人澡人人爽秒播| 成人无遮挡网站| 久久久久免费精品人妻一区二区| 亚洲熟妇熟女久久| 又爽又黄a免费视频| 国产av在哪里看| 久久久久九九精品影院| 少妇高潮的动态图| 日本免费a在线| 99热这里只有是精品50| 男女做爰动态图高潮gif福利片| 在线免费观看不下载黄p国产| av天堂中文字幕网| 日日干狠狠操夜夜爽| 99久久精品一区二区三区| 热99在线观看视频| 晚上一个人看的免费电影| 12—13女人毛片做爰片一| 久久韩国三级中文字幕| 国产精品亚洲美女久久久| 国产高清视频在线观看网站| 国产乱人偷精品视频| 免费看a级黄色片| 久久午夜亚洲精品久久| 国产精品国产高清国产av| 男女边吃奶边做爰视频| 黄色欧美视频在线观看| 色噜噜av男人的天堂激情| 国产美女午夜福利| 日本黄大片高清| 久久韩国三级中文字幕| 国产精品乱码一区二三区的特点| 乱人视频在线观看| 少妇的逼水好多| 久久久久久伊人网av| 女人被狂操c到高潮| 国产免费一级a男人的天堂| 女人十人毛片免费观看3o分钟| 真人做人爱边吃奶动态| 免费高清视频大片| 亚洲图色成人| 国产不卡一卡二| 国产中年淑女户外野战色| 亚洲电影在线观看av| 日日摸夜夜添夜夜爱| 免费av观看视频| 日韩三级伦理在线观看| 亚洲人成网站在线播| 欧美又色又爽又黄视频| 久久久色成人| 成年女人看的毛片在线观看| 中国美女看黄片| 国产伦一二天堂av在线观看| 久久久精品欧美日韩精品| 日韩大尺度精品在线看网址| 亚州av有码| 欧美另类亚洲清纯唯美| 婷婷色综合大香蕉| 成年免费大片在线观看| 99国产极品粉嫩在线观看| 国产欧美日韩一区二区精品| 免费观看在线日韩| 一边摸一边抽搐一进一小说| h日本视频在线播放| 深爱激情五月婷婷| 啦啦啦观看免费观看视频高清| 精品一区二区免费观看| 国产精品永久免费网站| 最新中文字幕久久久久| 特大巨黑吊av在线直播| 日本色播在线视频| 免费看日本二区| 亚洲精品日韩在线中文字幕 | 成人综合一区亚洲| 五月玫瑰六月丁香| 国产成人a区在线观看| 久久婷婷人人爽人人干人人爱| a级毛色黄片| 免费在线观看影片大全网站| 好男人在线观看高清免费视频| 在线观看66精品国产| 成人午夜高清在线视频| 黄色配什么色好看| 亚洲欧美日韩卡通动漫| 日韩精品中文字幕看吧| 亚洲av熟女| 久久久久免费精品人妻一区二区| 听说在线观看完整版免费高清| 啦啦啦韩国在线观看视频| 日日撸夜夜添| 69av精品久久久久久| 亚洲精华国产精华液的使用体验 | 午夜精品在线福利| 寂寞人妻少妇视频99o| 成人无遮挡网站| 婷婷亚洲欧美| 久久6这里有精品| 色5月婷婷丁香| 成人特级av手机在线观看| 久久久精品大字幕| 99热网站在线观看| 别揉我奶头 嗯啊视频| 3wmmmm亚洲av在线观看| 色av中文字幕| 国产一区二区在线观看日韩| 成人亚洲精品av一区二区| 深夜a级毛片| 久久人人爽人人片av| 久久热精品热| 精品无人区乱码1区二区| 国产高清视频在线播放一区| 成人av在线播放网站| 身体一侧抽搐| 久久久国产成人免费| 精品一区二区三区视频在线观看免费| 国内精品宾馆在线| 两个人的视频大全免费| 国产精品一区www在线观看| 国产精品一区二区性色av| 禁无遮挡网站| 国产精品一区二区性色av| 男人的好看免费观看在线视频| 在线观看66精品国产| 六月丁香七月| 女同久久另类99精品国产91| 午夜福利视频1000在线观看| 色哟哟哟哟哟哟| 国产精华一区二区三区| 日本黄大片高清| 久久九九热精品免费| 亚洲人成网站在线播放欧美日韩| 性色avwww在线观看| 三级男女做爰猛烈吃奶摸视频| 69av精品久久久久久| 国内精品一区二区在线观看| 夜夜爽天天搞| 国产视频一区二区在线看| av卡一久久| 又黄又爽又免费观看的视频| 一个人观看的视频www高清免费观看| 日韩精品中文字幕看吧| 亚洲高清免费不卡视频| 国产精品久久久久久久电影| 亚洲欧美日韩高清专用| 老熟妇仑乱视频hdxx| 国产真实乱freesex| 日韩欧美 国产精品| 国产欧美日韩一区二区精品| 国产乱人视频| 在线免费十八禁| 免费av不卡在线播放| 成人特级av手机在线观看| 秋霞在线观看毛片| 午夜视频国产福利| 亚洲精品在线观看二区| 毛片一级片免费看久久久久| 日本一本二区三区精品| 精品人妻视频免费看| 色视频www国产| 村上凉子中文字幕在线| 成熟少妇高潮喷水视频| 色综合亚洲欧美另类图片| 级片在线观看| 欧美精品国产亚洲| 麻豆一二三区av精品| 九九热线精品视视频播放| 中文字幕av在线有码专区| 淫秽高清视频在线观看| 精品少妇黑人巨大在线播放 | 亚洲国产欧美人成| av在线亚洲专区| 亚洲av一区综合| 精品少妇黑人巨大在线播放 | 激情 狠狠 欧美| 亚洲无线观看免费| av视频在线观看入口| 成人永久免费在线观看视频| 久久久久久久久久黄片| 少妇裸体淫交视频免费看高清| 久久精品影院6| 国产69精品久久久久777片| 一个人看的www免费观看视频| 欧美一区二区亚洲| 精品国内亚洲2022精品成人| 亚洲av.av天堂| 男人狂女人下面高潮的视频| 一区二区三区高清视频在线| 午夜福利18| 久久精品国产亚洲网站| 不卡一级毛片| 一本久久中文字幕| 亚洲av二区三区四区| 午夜亚洲福利在线播放| 日韩成人伦理影院| 国国产精品蜜臀av免费| 99久久成人亚洲精品观看| 亚洲av中文av极速乱| 国产精品一区www在线观看| 国产老妇女一区| 色播亚洲综合网| 久久精品国产亚洲av涩爱 | 亚洲欧美精品自产自拍| 免费人成视频x8x8入口观看| 赤兔流量卡办理| 国产一区亚洲一区在线观看| 直男gayav资源| 少妇人妻一区二区三区视频| 可以在线观看的亚洲视频| 在线观看66精品国产| 日本免费a在线| 在线播放无遮挡| 精品久久久久久久久av| 久久久精品大字幕| 久久精品国产清高在天天线| 亚洲欧美精品综合久久99| 亚洲国产高清在线一区二区三| 一本久久中文字幕| 午夜福利在线观看免费完整高清在 | 久久久久久国产a免费观看| 黄色欧美视频在线观看| 久久精品影院6| 黄色欧美视频在线观看| 日韩三级伦理在线观看| 中文字幕免费在线视频6| 一级a爱片免费观看的视频| 精品无人区乱码1区二区| 久久欧美精品欧美久久欧美| 一进一出抽搐动态| 欧美绝顶高潮抽搐喷水| 露出奶头的视频| 日韩成人伦理影院| 长腿黑丝高跟| 午夜精品一区二区三区免费看| 日韩三级伦理在线观看| 直男gayav资源| 国产伦精品一区二区三区视频9| 亚洲图色成人| 婷婷六月久久综合丁香| 国产av一区在线观看免费| 91麻豆精品激情在线观看国产| 男人的好看免费观看在线视频| 国产欧美日韩一区二区精品| 91狼人影院| 国产真实伦视频高清在线观看| 日本一二三区视频观看| 在线观看免费视频日本深夜| 午夜福利在线观看吧| 91狼人影院| 乱人视频在线观看| 国产伦在线观看视频一区| 在线观看美女被高潮喷水网站| 精品国内亚洲2022精品成人| 久久午夜福利片| av天堂中文字幕网| 日本一二三区视频观看| 婷婷精品国产亚洲av| 美女内射精品一级片tv| 亚洲专区国产一区二区| 免费搜索国产男女视频| 白带黄色成豆腐渣| av卡一久久| 白带黄色成豆腐渣| 久久韩国三级中文字幕| 在线观看免费视频日本深夜| 亚洲va在线va天堂va国产| 在线观看66精品国产| 卡戴珊不雅视频在线播放| 伦精品一区二区三区| 国产精品一及| 成年av动漫网址| 国产 一区精品| 免费看av在线观看网站| 日产精品乱码卡一卡2卡三| 又爽又黄无遮挡网站| 五月伊人婷婷丁香| 精品日产1卡2卡| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲中文字幕一区二区三区有码在线看| 国产高清有码在线观看视频| 少妇熟女aⅴ在线视频| 少妇裸体淫交视频免费看高清| 久久久久久久久久成人| 免费看美女性在线毛片视频| 搞女人的毛片| 欧美在线一区亚洲| .国产精品久久| 最近手机中文字幕大全| 久久鲁丝午夜福利片| 色哟哟·www| 欧美潮喷喷水| 五月玫瑰六月丁香| 中文字幕人妻熟人妻熟丝袜美| 日日撸夜夜添| 国产色爽女视频免费观看| 深夜精品福利| 真实男女啪啪啪动态图| 亚州av有码| 成人精品一区二区免费| 国产精品精品国产色婷婷| 久久精品国产亚洲av涩爱 | 久久午夜福利片| 在线观看66精品国产| 色播亚洲综合网| .国产精品久久| av在线天堂中文字幕| 国产午夜精品论理片| 欧美在线一区亚洲| 99视频精品全部免费 在线| 美女cb高潮喷水在线观看| 三级国产精品欧美在线观看| 国产黄色视频一区二区在线观看 | 国产视频内射| 色噜噜av男人的天堂激情| 别揉我奶头 嗯啊视频| 我的女老师完整版在线观看| 国产在线精品亚洲第一网站| 亚洲精品色激情综合| 午夜日韩欧美国产| 日本a在线网址| 一级黄色大片毛片| av女优亚洲男人天堂| 亚洲精品在线观看二区| 久久草成人影院| 日韩精品有码人妻一区| 精品久久久久久久久久久久久| 国产亚洲91精品色在线| 大又大粗又爽又黄少妇毛片口| 亚洲精品国产成人久久av| 老司机午夜福利在线观看视频| 久久久成人免费电影| 青春草视频在线免费观看| 黄色视频,在线免费观看| 黄色一级大片看看| 99riav亚洲国产免费| 午夜福利高清视频| 亚洲欧美日韩高清专用| 国内少妇人妻偷人精品xxx网站| 秋霞在线观看毛片| 欧美日韩在线观看h| 国产91av在线免费观看| 亚洲综合色惰| 久久久久久久亚洲中文字幕| 国产国拍精品亚洲av在线观看| 天堂av国产一区二区熟女人妻| 中国国产av一级| 网址你懂的国产日韩在线| 精品欧美国产一区二区三| 亚洲精品日韩在线中文字幕 | av在线老鸭窝| 久久精品综合一区二区三区| 青春草视频在线免费观看| a级毛片免费高清观看在线播放| 欧美人与善性xxx| 亚洲av第一区精品v没综合| 黄色一级大片看看| 性插视频无遮挡在线免费观看| 不卡一级毛片| 夜夜夜夜夜久久久久| 久久这里只有精品中国| 免费人成在线观看视频色| 一级毛片久久久久久久久女| 国产精品乱码一区二三区的特点| 美女cb高潮喷水在线观看| 精品不卡国产一区二区三区| 国产淫片久久久久久久久| 久久热精品热| 欧美中文日本在线观看视频| 草草在线视频免费看| 老女人水多毛片| 麻豆国产97在线/欧美| 久久久久免费精品人妻一区二区| 国产一区亚洲一区在线观看| av在线亚洲专区| 国产精品久久视频播放| 国产高潮美女av| 老司机午夜福利在线观看视频| 国产精品,欧美在线| 小蜜桃在线观看免费完整版高清| 18禁在线播放成人免费| 国产伦在线观看视频一区| 少妇猛男粗大的猛烈进出视频 | 人妻制服诱惑在线中文字幕| 最近手机中文字幕大全| 日产精品乱码卡一卡2卡三| 桃色一区二区三区在线观看| 香蕉av资源在线| 欧美日本视频| 一区二区三区四区激情视频 | 欧美最新免费一区二区三区| 国产精品久久视频播放| 欧美一区二区精品小视频在线| 欧美成人精品欧美一级黄| 国产私拍福利视频在线观看| 免费av观看视频|