• <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)率改性聚苯胺的合成新工藝
    在线天堂中文资源库| 亚洲精品aⅴ在线观看| 只有这里有精品99| 免费久久久久久久精品成人欧美视频 | 精品久久久久久电影网| 久久人人爽av亚洲精品天堂| 激情视频va一区二区三区| 国产又爽黄色视频| 制服丝袜香蕉在线| 日本vs欧美在线观看视频| 国产高清国产精品国产三级| 麻豆乱淫一区二区| 80岁老熟妇乱子伦牲交| 亚洲精品视频女| 久久人人爽人人爽人人片va| 亚洲国产色片| 久久这里有精品视频免费| 中文字幕另类日韩欧美亚洲嫩草| av免费在线看不卡| 国产精品久久久久成人av| 国产成人欧美| 亚洲久久久国产精品| 日韩av不卡免费在线播放| 国产一区亚洲一区在线观看| 一本久久精品| 涩涩av久久男人的天堂| 免费观看在线日韩| 久久久国产一区二区| 最近手机中文字幕大全| 在线观看人妻少妇| 色网站视频免费| 最新中文字幕久久久久| 日本av手机在线免费观看| 亚洲精品第二区| 国产淫语在线视频| 久久久亚洲精品成人影院| 亚洲国产欧美在线一区| 中文字幕精品免费在线观看视频 | 精品久久蜜臀av无| 男人爽女人下面视频在线观看| 丝袜喷水一区| 天天操日日干夜夜撸| 少妇的丰满在线观看| 成人无遮挡网站| av线在线观看网站| 欧美日韩一区二区视频在线观看视频在线| 99热网站在线观看| 久久精品国产综合久久久 | 久久人人爽人人爽人人片va| 狂野欧美激情性bbbbbb| a级片在线免费高清观看视频| 日韩av在线免费看完整版不卡| av又黄又爽大尺度在线免费看| 午夜老司机福利剧场| 三级国产精品片| 亚洲av综合色区一区| 天天躁夜夜躁狠狠躁躁| 精品亚洲成国产av| 亚洲精品久久成人aⅴ小说| 肉色欧美久久久久久久蜜桃| 日韩av免费高清视频| 免费日韩欧美在线观看| 亚洲国产看品久久| 最新的欧美精品一区二区| 伦理电影免费视频| 如何舔出高潮| 日韩熟女老妇一区二区性免费视频| 久久精品国产自在天天线| 久久久久久久精品精品| 侵犯人妻中文字幕一二三四区| 久久久a久久爽久久v久久| 插逼视频在线观看| 国产男女内射视频| 国产精品国产三级国产专区5o| 少妇的逼水好多| 日韩人妻精品一区2区三区| 九色成人免费人妻av| 男人爽女人下面视频在线观看| 成人综合一区亚洲| 国产精品 国内视频| 国产亚洲午夜精品一区二区久久| 国产 精品1| 视频在线观看一区二区三区| 日韩中字成人| 亚洲欧美清纯卡通| 国产色婷婷99| 日韩av不卡免费在线播放| 亚洲精品国产色婷婷电影| 少妇熟女欧美另类| 欧美激情极品国产一区二区三区 | 九九爱精品视频在线观看| 女性生殖器流出的白浆| 久久久久视频综合| 中文乱码字字幕精品一区二区三区| 日韩电影二区| 街头女战士在线观看网站| 免费大片18禁| 黄色 视频免费看| 日本黄色日本黄色录像| 色婷婷av一区二区三区视频| 夜夜骑夜夜射夜夜干| 97在线人人人人妻| 成人毛片60女人毛片免费| 九色成人免费人妻av| 久久久久国产精品人妻一区二区| 人人妻人人爽人人添夜夜欢视频| 草草在线视频免费看| av国产久精品久网站免费入址| videos熟女内射| 免费观看在线日韩| 多毛熟女@视频| 少妇熟女欧美另类| 亚洲国产日韩一区二区| 又大又黄又爽视频免费| 亚洲精品日本国产第一区| 成年动漫av网址| 亚洲激情五月婷婷啪啪| 国产成人免费无遮挡视频| 日日啪夜夜爽| 91午夜精品亚洲一区二区三区| 青春草亚洲视频在线观看| 免费看av在线观看网站| 国产免费现黄频在线看| 亚洲内射少妇av| 人妻系列 视频| 国产成人一区二区在线| 欧美bdsm另类| 高清欧美精品videossex| 久久国产精品男人的天堂亚洲 | 欧美成人精品欧美一级黄| 在线天堂中文资源库| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美清纯卡通| 尾随美女入室| 精品国产乱码久久久久久小说| 精品少妇内射三级| 熟妇人妻不卡中文字幕| 精品人妻在线不人妻| 欧美精品一区二区免费开放| 午夜老司机福利剧场| 日韩视频在线欧美| 国产熟女午夜一区二区三区| 日本欧美国产在线视频| 久久久久久久久久成人| 成人国产av品久久久| 欧美少妇被猛烈插入视频| 激情视频va一区二区三区| 国产在线免费精品| 国产成人免费观看mmmm| 少妇人妻精品综合一区二区| 久久久国产精品麻豆| 一本色道久久久久久精品综合| 曰老女人黄片| 亚洲在久久综合| 亚洲伊人久久精品综合| 免费黄色在线免费观看| videosex国产| 丰满乱子伦码专区| 十八禁高潮呻吟视频| 久久99一区二区三区| 国产视频首页在线观看| 亚洲av在线观看美女高潮| 日本黄大片高清| 亚洲av电影在线观看一区二区三区| 飞空精品影院首页| 美女脱内裤让男人舔精品视频| 有码 亚洲区| 欧美3d第一页| 日本黄大片高清| 亚洲国产精品一区三区| tube8黄色片| a级毛片在线看网站| 在线观看美女被高潮喷水网站| 99久久人妻综合| 亚洲一码二码三码区别大吗| 一二三四中文在线观看免费高清| 老女人水多毛片| 精品人妻偷拍中文字幕| 咕卡用的链子| 人人妻人人爽人人添夜夜欢视频| www.色视频.com| 97人妻天天添夜夜摸| 在线精品无人区一区二区三| 成人国语在线视频| 国产精品蜜桃在线观看| 亚洲精品乱码久久久久久按摩| av国产精品久久久久影院| 国产亚洲av片在线观看秒播厂| 永久免费av网站大全| 午夜免费鲁丝| 亚洲第一av免费看| 久久久久精品久久久久真实原创| 久久精品夜色国产| 看免费成人av毛片| www.熟女人妻精品国产 | xxx大片免费视频| 性色av一级| 九九在线视频观看精品| 国产精品欧美亚洲77777| 亚洲国产欧美日韩在线播放| 中文天堂在线官网| 亚洲精品久久午夜乱码| 免费人成在线观看视频色| 制服诱惑二区| 丝袜脚勾引网站| 国产欧美亚洲国产| 亚洲精品日韩在线中文字幕| 久久国产亚洲av麻豆专区| 国产日韩欧美亚洲二区| 国产一区二区三区av在线| 日韩免费高清中文字幕av| 成人午夜精彩视频在线观看| av片东京热男人的天堂| 天堂俺去俺来也www色官网| 最近中文字幕2019免费版| www.熟女人妻精品国产 | 少妇的逼好多水| 女的被弄到高潮叫床怎么办| 1024视频免费在线观看| 热re99久久国产66热| 青春草视频在线免费观看| 国产一区二区三区av在线| 大香蕉97超碰在线| 国产精品.久久久| 国产激情久久老熟女| 如何舔出高潮| 免费高清在线观看视频在线观看| 免费日韩欧美在线观看| 夫妻性生交免费视频一级片| www.熟女人妻精品国产 | 亚洲欧洲日产国产| 妹子高潮喷水视频| 午夜av观看不卡| 搡老乐熟女国产| 好男人视频免费观看在线| 水蜜桃什么品种好| 不卡视频在线观看欧美| 99热国产这里只有精品6| 亚洲av福利一区| 青青草视频在线视频观看| 美女福利国产在线| 丰满少妇做爰视频| 久久久久精品性色| 成人综合一区亚洲| 精品国产露脸久久av麻豆| 日韩制服骚丝袜av| 久久精品久久精品一区二区三区| 国产在线一区二区三区精| 日本黄大片高清| 观看美女的网站| 国产精品久久久久久久久免| 免费女性裸体啪啪无遮挡网站| 伦理电影免费视频| 丝袜脚勾引网站| 日本91视频免费播放| 男女下面插进去视频免费观看 | 久久99热6这里只有精品| 亚洲在久久综合| 久久精品久久久久久噜噜老黄| 少妇被粗大的猛进出69影院 | 日韩一本色道免费dvd| 九色成人免费人妻av| 侵犯人妻中文字幕一二三四区| 在线观看免费视频网站a站| 亚洲精品aⅴ在线观看| 中文字幕最新亚洲高清| 韩国av在线不卡| 亚洲国产精品999| 亚洲精品,欧美精品| 久久午夜综合久久蜜桃| 九色亚洲精品在线播放| 久久久欧美国产精品| 婷婷色麻豆天堂久久| 美女xxoo啪啪120秒动态图| 婷婷色av中文字幕| 伦理电影大哥的女人| 女人精品久久久久毛片| 最近手机中文字幕大全| 男女无遮挡免费网站观看| 91午夜精品亚洲一区二区三区| 搡老乐熟女国产| 香蕉精品网在线| 亚洲av电影在线进入| 日韩大片免费观看网站| 久久久久精品性色| 午夜视频国产福利| 欧美精品国产亚洲| tube8黄色片| 视频区图区小说| 久久精品久久久久久久性| 满18在线观看网站| 我要看黄色一级片免费的| 免费看光身美女| 免费久久久久久久精品成人欧美视频 | 国产日韩欧美在线精品| 99久久中文字幕三级久久日本| 国产精品久久久久久久久免| 国产精品无大码| 狂野欧美激情性bbbbbb| 欧美xxⅹ黑人| 大香蕉97超碰在线| 中文字幕另类日韩欧美亚洲嫩草| 丰满少妇做爰视频| 国产男女超爽视频在线观看| 最近中文字幕高清免费大全6| 午夜免费男女啪啪视频观看| 一本色道久久久久久精品综合| 精品99又大又爽又粗少妇毛片| 国产极品天堂在线| 久久久久国产精品人妻一区二区| 久久女婷五月综合色啪小说| 日韩中文字幕视频在线看片| 日日爽夜夜爽网站| 免费观看av网站的网址| 人人妻人人澡人人看| 国产av国产精品国产| 国产亚洲精品第一综合不卡 | 日本vs欧美在线观看视频| 99九九在线精品视频| 国产又色又爽无遮挡免| 国产国语露脸激情在线看| 天堂8中文在线网| 纯流量卡能插随身wifi吗| 精品国产露脸久久av麻豆| 午夜久久久在线观看| 狂野欧美激情性bbbbbb| 97超碰精品成人国产| 国产一区二区三区av在线| 国产成人免费无遮挡视频| 国产亚洲精品第一综合不卡 | 9热在线视频观看99| 99re6热这里在线精品视频| 亚洲综合色惰| 水蜜桃什么品种好| 国产一区二区三区综合在线观看 | 国产精品久久久久久精品电影小说| 午夜视频国产福利| av福利片在线| 欧美 亚洲 国产 日韩一| 午夜免费男女啪啪视频观看| 视频区图区小说| 一边摸一边做爽爽视频免费| 精品国产一区二区三区久久久樱花| 少妇人妻久久综合中文| 又黄又粗又硬又大视频| 秋霞伦理黄片| 观看av在线不卡| 亚洲av免费高清在线观看| 99re6热这里在线精品视频| 日韩不卡一区二区三区视频在线| 亚洲成人手机| 99久久人妻综合| 国产av码专区亚洲av| 少妇人妻久久综合中文| 亚洲欧美清纯卡通| 亚洲三级黄色毛片| 亚洲色图 男人天堂 中文字幕 | 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 你懂的网址亚洲精品在线观看| 丰满迷人的少妇在线观看| 免费少妇av软件| 男女啪啪激烈高潮av片| 国国产精品蜜臀av免费| 亚洲精品国产av蜜桃| 成人毛片60女人毛片免费| 毛片一级片免费看久久久久| 亚洲精品一二三| 99久久中文字幕三级久久日本| 一本大道久久a久久精品| 91精品三级在线观看| xxx大片免费视频| 国产爽快片一区二区三区| 午夜日本视频在线| 91精品伊人久久大香线蕉| 亚洲人成网站在线观看播放| 亚洲精品色激情综合| 日日啪夜夜爽| 国产一区二区三区av在线| 九色亚洲精品在线播放| 欧美日本中文国产一区发布| 欧美少妇被猛烈插入视频| 黑人巨大精品欧美一区二区蜜桃 | 黄色 视频免费看| 天堂中文最新版在线下载| 久久午夜综合久久蜜桃| 丝瓜视频免费看黄片| 国产 一区精品| 国产精品久久久久久精品电影小说| 色网站视频免费| 国产免费视频播放在线视频| 赤兔流量卡办理| xxx大片免费视频| 国产乱人偷精品视频| 视频在线观看一区二区三区| 黄色毛片三级朝国网站| 亚洲欧美色中文字幕在线| 亚洲av国产av综合av卡| 国精品久久久久久国模美| 91成人精品电影| 久久久亚洲精品成人影院| 午夜福利视频精品| 蜜桃在线观看..| 十分钟在线观看高清视频www| 国产69精品久久久久777片| 免费观看av网站的网址| 久久99热这里只频精品6学生| 国产又爽黄色视频| 精品国产乱码久久久久久小说| 曰老女人黄片| 免费av不卡在线播放| 人人澡人人妻人| 国产一区有黄有色的免费视频| 国产一区二区在线观看日韩| 91精品国产国语对白视频| 国产xxxxx性猛交| 26uuu在线亚洲综合色| 熟妇人妻不卡中文字幕| 国产精品国产三级国产av玫瑰| 亚洲第一av免费看| 丰满乱子伦码专区| 久久av网站| 啦啦啦中文免费视频观看日本| 五月开心婷婷网| 9191精品国产免费久久| 永久网站在线| 51国产日韩欧美| 丝袜人妻中文字幕| 国产精品久久久久久久久免| 精品久久久久久电影网| 成人综合一区亚洲| 制服人妻中文乱码| 五月开心婷婷网| 午夜老司机福利剧场| 男女免费视频国产| 一边摸一边做爽爽视频免费| 日韩中字成人| 男女午夜视频在线观看 | 日韩一区二区视频免费看| 搡老乐熟女国产| 国产在线视频一区二区| 成人无遮挡网站| 亚洲av日韩在线播放| 国产成人aa在线观看| 欧美亚洲 丝袜 人妻 在线| 香蕉精品网在线| 欧美日韩精品成人综合77777| 国产精品一二三区在线看| 亚洲av欧美aⅴ国产| 看十八女毛片水多多多| 99久久中文字幕三级久久日本| 欧美成人午夜精品| 久久人妻熟女aⅴ| 成人18禁高潮啪啪吃奶动态图| 亚洲国产看品久久| 国产成人免费观看mmmm| 国产精品蜜桃在线观看| 免费高清在线观看视频在线观看| 国产片特级美女逼逼视频| 国产有黄有色有爽视频| 日日摸夜夜添夜夜爱| 免费看光身美女| 丰满迷人的少妇在线观看| 一区二区日韩欧美中文字幕 | 亚洲情色 制服丝袜| 日本-黄色视频高清免费观看| 各种免费的搞黄视频| 久久精品国产a三级三级三级| 国产 一区精品| 国产亚洲一区二区精品| 亚洲精品视频女| 久久久久精品性色| 亚洲国产精品一区二区三区在线| 精品久久久精品久久久| 国产成人精品在线电影| 卡戴珊不雅视频在线播放| 国产精品偷伦视频观看了| 97人妻天天添夜夜摸| 免费观看av网站的网址| 精品国产国语对白av| 亚洲av电影在线进入| 卡戴珊不雅视频在线播放| 亚洲在久久综合| videossex国产| 成人毛片60女人毛片免费| 建设人人有责人人尽责人人享有的| 美女主播在线视频| 国产成人精品一,二区| 亚洲精品国产色婷婷电影| 九草在线视频观看| 老司机亚洲免费影院| 欧美亚洲 丝袜 人妻 在线| 成年av动漫网址| 777米奇影视久久| 亚洲,欧美精品.| av卡一久久| 制服丝袜香蕉在线| 热re99久久精品国产66热6| 一区二区av电影网| 国产淫语在线视频| 亚洲 欧美一区二区三区| 香蕉丝袜av| 最近2019中文字幕mv第一页| 观看av在线不卡| 午夜福利视频在线观看免费| 伦理电影免费视频| 国产精品一区二区在线观看99| 99久国产av精品国产电影| 欧美日韩成人在线一区二区| 乱人伦中国视频| √禁漫天堂资源中文www| 十八禁网站网址无遮挡| 国产亚洲精品久久久com| 亚洲av国产av综合av卡| 99久久人妻综合| 国产精品三级大全| 少妇高潮的动态图| 久久久久精品性色| 看非洲黑人一级黄片| 国产xxxxx性猛交| 欧美日韩视频高清一区二区三区二| 欧美日韩综合久久久久久| 国产 精品1| 26uuu在线亚洲综合色| 亚洲三级黄色毛片| 啦啦啦中文免费视频观看日本| 亚洲第一区二区三区不卡| 国产亚洲一区二区精品| 99久久中文字幕三级久久日本| 国产不卡av网站在线观看| 美女内射精品一级片tv| 一级片免费观看大全| 亚洲精品国产av蜜桃| 国产亚洲欧美精品永久| 在线天堂中文资源库| www日本在线高清视频| 欧美日韩视频精品一区| 国产成人精品一,二区| 天堂8中文在线网| 美女脱内裤让男人舔精品视频| 韩国av在线不卡| 亚洲综合色惰| 国产伦理片在线播放av一区| 美国免费a级毛片| 2022亚洲国产成人精品| 97在线人人人人妻| 久久影院123| 黑丝袜美女国产一区| 久久人人爽人人片av| 99热网站在线观看| 免费看不卡的av| 国产精品久久久久久av不卡| 男女啪啪激烈高潮av片| 国产乱人偷精品视频| 青春草亚洲视频在线观看| 插逼视频在线观看| 欧美3d第一页| 中文精品一卡2卡3卡4更新| 亚洲国产日韩一区二区| 人妻少妇偷人精品九色| 一区二区av电影网| kizo精华| 又黄又粗又硬又大视频| 亚洲天堂av无毛| 男女无遮挡免费网站观看| 高清不卡的av网站| 久久久亚洲精品成人影院| av线在线观看网站| av天堂久久9| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av综合色区一区| 亚洲国产欧美在线一区| 国产在线免费精品| 一二三四在线观看免费中文在 | 一二三四中文在线观看免费高清| av有码第一页| 美国免费a级毛片| 成人毛片a级毛片在线播放| 国产日韩欧美视频二区| 亚洲国产av影院在线观看| 宅男免费午夜| 99久国产av精品国产电影| 中文字幕人妻熟女乱码| 国产av一区二区精品久久| 嫩草影院入口| 一区二区三区精品91| 天天影视国产精品| 91精品国产国语对白视频| 婷婷色麻豆天堂久久| 午夜福利视频在线观看免费| 久久久久久久精品精品| 伊人亚洲综合成人网| 伊人久久国产一区二区| 两个人看的免费小视频| 国产福利在线免费观看视频| 亚洲美女搞黄在线观看| 青春草亚洲视频在线观看| 国产免费福利视频在线观看| 中文字幕免费在线视频6| 欧美精品av麻豆av| 熟女电影av网| 女的被弄到高潮叫床怎么办| 丰满乱子伦码专区| 国产精品久久久久久久电影| 蜜桃国产av成人99| 夫妻午夜视频| 香蕉丝袜av| 亚洲精品国产av蜜桃| av天堂久久9| 男女免费视频国产| 久久精品国产亚洲av天美| 久久人人爽av亚洲精品天堂| 大话2 男鬼变身卡| 国产深夜福利视频在线观看| 美女xxoo啪啪120秒动态图| 午夜福利网站1000一区二区三区|