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

    黃河三角洲引黃灌區(qū)土壤鹽分遙感反演及空間分布特征

    2025-06-28 00:00:00齊翟碩傅新張浩然黃瑾龐桂斌徐征和
    關(guān)鍵詞:模型研究

    Remote Sensing Inversion and Spatial Distribution Characteristics of Soil Salinity in the Yellow River Irrigation Area of the Yellow River Delta

    QI Zhaishuo,F(xiàn)U Xin, ZHANG Haoran,HUANG Jin,PANG Guibin,XU Zhenghe(Schoolof Water Conservancyand Environment,Universityof Jinan,Jinan 25oo22,Shandong,China)

    Abstract:To accurately graspthespatialdistributioncharacteristicsofsoilsalinityintheYellowRiverDelta rrigationarea, Kendong Irrigation District in Binzhou City,Shandong Provincewas chosenas thestudyarea,basedonthe soilsalinitydata fromfield surveyand Sentinel-2B multispectralremotesensing data,thespectral indices thathada strong explanatory powerof soil salinitywereselected using geodetic detector.Thepartial least square methodwas thenusedto construct remote sensing inversion modelsof soil salinityand explorethe spatialdistributioncharacteristicsof soilsalinity in the study area.The results show that the optimal spectral indices with high explanatory power ( qgt;0.05 )andhighly significant level ( P lt;0.01)for soil salinity are visible-infrared salinity index,extended diffrence vegetation index,normalized diference waterindex,nonlinear index,S6salinityindex,brightness index,reflectivityof sateliteB8Aband,INT2 salinity idex,optimized soil adjusted vegetation index,soil adjusted vegetation index,enhancedvegetation index,SI3 salinityindex,extendedratiovegetationindex,salinityindex,atmosphericalyresistantvegetationindex,S5salinityindex. The determinationcoefficient and meanrelative erorof theremote sensing inversionmodel of soil salinityareO.77and 0.36,respectively,the model can predict the salinitycondition in the study area more accurately.The studyarea consists of non-saline soil,mildly saline soil and moderately saline soil. The area with a salt mass ratio less than 0.2g/kg in nonsaline soil accounts for the largest proportion,which is 50.43% . The mildly saline soil and moderately saline soil show a fragmenteddistribution,mainlyoccuring in the non-planting area,which is inline with theresultsoffieldsurvey.The average soil saltcontentof differentlanduse-covertypes is significantly different,andtheorder ofland use-cover types withtheaverage saltcontentfrom large to smallis saline wasteland,abandoned land-summer corn field,forestland, winter wheat-summer corn field,and paddy field.

    Keywords : soil salinity;remote sensing retrieval; spatial distribution; geodetector; the Yellow River Delta

    土壤鹽漬化是土壤表層集聚可溶性鹽的一種現(xiàn)象,極大地限制了農(nóng)業(yè)發(fā)展和土地資源的可持續(xù)利用,因此土壤鹽漬化已經(jīng)被列為中國(guó)生態(tài)文明建設(shè)和全球變化議題中的重要研究方向。黃河三角洲地區(qū)作為我國(guó)重要的后備土地資源區(qū),蘊(yùn)藏著豐富的自然資源,開(kāi)發(fā)前景極為廣闊;但是,由于長(zhǎng)期海水入滲以及相對(duì)干旱的氣候環(huán)境,因此王壤鹽漬化問(wèn)題日益嚴(yán)重,尤其是在以灌溉為主要方式的農(nóng)業(yè)灌區(qū),重灌溉輕排放的農(nóng)耕形式對(duì)土壤環(huán)境造成了嚴(yán)重的威脅,使得黃河下游引黃灌區(qū)土壤鹽漬化治理變得更加困難。

    在灌溉農(nóng)業(yè)區(qū),監(jiān)測(cè)土壤鹽漬化分布特征可以很好地掌握灌溉區(qū)內(nèi)鹽分的變化情況,制定科學(xué)、合理的水資源管理方案,為灌溉農(nóng)業(yè)區(qū)土壤鹽漬化治理和綜合利用提供科學(xué)依據(jù)。目前,通過(guò)定量遙感手段監(jiān)測(cè)土壤鹽分信息是土壤鹽漬化監(jiān)測(cè)研究的熱點(diǎn)之一。楊小虎等2基于Landsat-8號(hào)衛(wèi)星攜帶的陸地成像儀(Landsat-8OLI)多光譜遙感影像和野外實(shí)測(cè)土壤鹽分?jǐn)?shù)據(jù),通過(guò)歸一化差異水體指數(shù)(INDW)等光譜指數(shù)與土壤鹽分構(gòu)建遙感預(yù)測(cè)模型,研究土壤鹽分空間分布特征,結(jié)果表明,模型精度較高且預(yù)測(cè)的土壤鹽分的空間分布規(guī)律與實(shí)地調(diào)查基本一致。陳紅艷等3在傳統(tǒng)植被指數(shù)的基礎(chǔ)上提出了改進(jìn)植被指數(shù),基于改進(jìn)前、后2種對(duì)應(yīng)的植被指數(shù)構(gòu)建土壤鹽分反演模型,獲得了較好的土壤鹽分空間分布反演結(jié)果。周曉紅等[4在傳統(tǒng)型鹽分指數(shù)和植被指數(shù)的基礎(chǔ)上,引入增強(qiáng)型比值植被指數(shù)與研究區(qū)土壤鹽分建立遙感估算模型,用于定量監(jiān)測(cè)研究區(qū)土壤鹽分的動(dòng)態(tài)變化規(guī)律。張同瑞等[5]采用多光譜衛(wèi)星遙感數(shù)據(jù),篩選敏感光譜指數(shù),最終遴選出土壤鹽分最佳監(jiān)測(cè)模型。張素銘等[通過(guò)相關(guān)性分析等手段篩選出最優(yōu)的土壤鹽分敏感波段和光譜參量,構(gòu)建最佳的土壤鹽分遙感反演模型,分析黃河三角洲地區(qū)土壤鹽分動(dòng)態(tài)分布特征,結(jié)果表明,遙感波段之間相互組合可以提高其與土壤鹽分的相關(guān)程度,并且在組合波段中加入敏感波段的建模方法表現(xiàn)更優(yōu)。劉恩等7以黃河三角洲小開(kāi)河引黃灌區(qū)為例,基于Landsat-8OLI多光譜影像,選取診斷指數(shù)較大的光譜反射率和土壤鹽分?jǐn)?shù)據(jù)構(gòu)建土壤鹽分反演模型。上述研究大多數(shù)是利用單一的光譜指數(shù)(如光譜反射率、植被指數(shù)等)或少數(shù)幾個(gè)光譜指數(shù)組合構(gòu)建土壤鹽分反演模型,用于監(jiān)測(cè)土壤鹽漬化信息,但估算精度仍然較低[8]。光譜指數(shù)選擇不夠全面,容易忽略指數(shù)間的耦合作用,使得模型失真和不穩(wěn)定,因此,不局限于少數(shù)光譜指數(shù),嘗試集成應(yīng)用多種光譜指數(shù)來(lái)定量反演土壤鹽分是一種值得探究的方式。對(duì)于濱海典型環(huán)境脆弱區(qū)黃河三角洲來(lái)說(shuō),以哨兵-2號(hào)衛(wèi)星攜帶的多光譜成像儀(Sentinel-2MSI)影像數(shù)據(jù)和地面實(shí)測(cè)數(shù)據(jù)為基礎(chǔ),利用系列光譜指數(shù)反演土壤鹽分的研究較少。Wang等利用Cubist模型對(duì)比分析Landsat-8OLI與Sentinel-2MSI這2種傳感器之間的差異,結(jié)果表明,空間分辨率和光譜分辨率越高的衛(wèi)星遙感影像越能夠有效提高模型反演精度。

    本文中以黃河三角洲引黃灌區(qū)中的墾東灌區(qū)為研究區(qū),基于野外土壤采樣數(shù)據(jù)和Sentinel-2BMSI遙感數(shù)據(jù),獲取不同遙感波段光譜反射率、所有的鹽分指數(shù)和植被指數(shù)等系列光譜指數(shù),利用地理探測(cè)器篩選對(duì)王壤鹽分解釋力較強(qiáng)的光譜指數(shù)因子,構(gòu)建最優(yōu)土壤鹽分反演模型,分析黃河三角洲引黃灌區(qū)土壤鹽分空間分布規(guī)律,以期為黃河三角洲引黃灌區(qū)土壤鹽漬化有效治理以及區(qū)域土地資源可持續(xù)利用提供理論依據(jù)

    數(shù)據(jù)來(lái)源與研究方法

    1. 1 研究區(qū)概況

    墾東灌區(qū)位于山東省東營(yíng)市墾利區(qū),地理位置為東經(jīng) 118°5455′′ 、北緯 37°4040′′ ,是現(xiàn)代黃河三角洲具有代表性的引黃灌區(qū)之一。研究區(qū)四季氣候變化較為明顯,冬季干冷,夏季濕熱,處于溫帶大陸性季風(fēng)氣候區(qū),多年平均氣溫為 11.7°C ,多年平均風(fēng)速為 3.85m/s ,平均年蒸發(fā)量為 2049.4mm ,年平均降水量為 578mm ,降雨多集中在每年6—9月份,一般占全年降水量的 70% 以上[10]。灌區(qū)內(nèi)土壤類型主要為粉土和粉壤土,地下水位和礦化度較高,土壤毛細(xì)作用劇烈導(dǎo)致土壤反鹽、聚鹽現(xiàn)象嚴(yán)重。

    1. 2 土壤樣品采集與指標(biāo)測(cè)定

    于2023年5月18—20日開(kāi)展墾東灌區(qū)的野外調(diào)查和土樣采集等工作。通過(guò)綜合分析研究灌區(qū)土壤表層的鹽漬化空間分布狀況、土壤采集點(diǎn)整體均勻程度和典型程度布設(shè)土壤樣品采集點(diǎn)位。野外用五點(diǎn)法采集土壤樣品,每個(gè)點(diǎn)位采用專業(yè)土鉆采集0~20cm 土層土壤。用手持式全球定位系統(tǒng)(GPS)獲取樣點(diǎn)經(jīng)、緯度,28個(gè)土壤樣點(diǎn)分布如圖1所示。整理好野外采集土樣并帶回實(shí)驗(yàn)室,然后將土壤樣品經(jīng)過(guò)風(fēng)干、磨碎、過(guò)孔徑 2mm 篩等常規(guī)操作后,以土樣與水的質(zhì)量比為1:5混合,經(jīng)振蕩、過(guò)濾后提取測(cè)定溶液,使用哈希多參數(shù)水質(zhì)測(cè)定儀測(cè)定土壤電導(dǎo)率和全鹽量

    圖1黃河三角洲引黃灌區(qū)土壤采樣點(diǎn)分布圖(資料來(lái)源:歐洲航天局哥白尼數(shù)據(jù)中心網(wǎng)站,經(jīng)ArcGIS10.8軟件數(shù)字化處理。)

    根據(jù)陳紅艷等[3]在黃河三角洲針對(duì)土壤鹽漬化的研究,將土壤鹽分含量劃分為5個(gè)等級(jí):鹽分質(zhì)量比小于 2g/kg 的土壤為非鹽漬土,鹽分質(zhì)量比為 2~lt;4g/kg 的土壤為輕度鹽漬土,鹽分質(zhì)量比為 4~lt;6g/kg 的土壤為中度鹽漬土,鹽分質(zhì)量比為6~lt;10g/kg 的土壤為重度鹽漬土,鹽分質(zhì)量比為10g/kg 的土壤為鹽土。

    1.3 遙感數(shù)據(jù)與預(yù)處理

    哨兵-2(Sentinel-2)號(hào)衛(wèi)星包括A、B小衛(wèi)星,2顆衛(wèi)星之間的協(xié)同作用極大地提高了若干環(huán)境和自然資源的應(yīng)用能力,如植被活力、土壤鹽分動(dòng)態(tài)、水質(zhì)以及地方、區(qū)域和全球范圍的氣候變化影響分析。多光譜成像儀(MSI)以空間分辨率為10、20、60m 分別在光譜帶中以大視場(chǎng)對(duì)地球表面反射率進(jìn)行成像。每個(gè)Sentinel-2MSI場(chǎng)景的掃描范圍為290km ,允許每 10d 對(duì)地球表面進(jìn)行一次全球覆蓋。MSI輻射性能以12位編碼,能夠以4095個(gè)數(shù)字采集圖像,確保輻射精度小于 5% 和出色的信噪比。依據(jù)遙感影像數(shù)據(jù)與地表土壤數(shù)據(jù)時(shí)間盡量保持一致的原則以及野外采樣期間衛(wèi)星遙感影像在研究區(qū)內(nèi)的成像質(zhì)量,本文中獲取2023年5月16日Sentinel-2B衛(wèi)星遙感影像作為遙感數(shù)據(jù),選擇經(jīng)過(guò)輻射定標(biāo)和大氣校正處理后的大氣底層反射率影像數(shù)據(jù),即Level-2A級(jí)別產(chǎn)品,利用ENVI5.3軟件統(tǒng)一設(shè)置遙感影像投影坐標(biāo),執(zhí)行數(shù)據(jù)重采樣、影像邊界裁剪、多波段組合,輸出多波段影像等常規(guī)操作。

    選取與土壤鹽分密切相關(guān)的光譜反射率、鹽分指數(shù)和植被指數(shù)等48種光譜指數(shù)。具體鹽分指數(shù)、植被指數(shù)及其計(jì)算公式如表1、2所示。

    表1鹽分指數(shù)及其計(jì)算公式
    注: λB2 、 λB3 、 λB4 、 λB8 、 λB11 、 λB12 分別為哨兵-2B衛(wèi)星B2、B3、B4、B8、B11、B12波段的反射率。
    表2植被指數(shù)及其計(jì)算公式
    續(xù)表
    注: λB2?λB3?λB4?λB8?λB11?λB12 分別為哨兵 -2B 衛(wèi)星B2、B3、B4、B8、B11、B12波段的反射率; L 為土壤調(diào)節(jié)因素,一般取 0.5

    1.4 研究方法

    1.4.1 地理探測(cè)器

    地理探測(cè)器是由王勁峰等[25提出的探測(cè)空間分異性和揭示驅(qū)動(dòng)性因素的統(tǒng)計(jì)方法。本文中主要采用因子探測(cè)進(jìn)行土壤鹽分遙感反演模型因子的篩選。根據(jù)不同影響因素對(duì)土壤鹽分的解釋能力 q 值篩選出影響土壤鹽分的因子, q 值越大說(shuō)明因子對(duì)土壤鹽分的影響程度越高。

    式中: L 為土壤鹽分或光譜指數(shù)的分層數(shù)或分區(qū)數(shù); Nh : σh2 分別為 h 層的單元數(shù)、土壤鹽分的方差; σ2 為土壤鹽分在地理探測(cè)器探測(cè)過(guò)程中的所有分層或分區(qū)的方差

    1.4.2 數(shù)據(jù)分析與檢驗(yàn)

    多元統(tǒng)計(jì)分析方法偏最小二乘法(PLSR)能夠在不排除多重共線變量的前提下有效反映多個(gè)自變量或因變量之間的相關(guān)性,且適用于小容量樣本的多元統(tǒng)計(jì)分析,是定量分析土壤鹽分對(duì)多種光譜指數(shù)組合響應(yīng)的可行方法。本文中將28個(gè)有效采樣點(diǎn)按照7:3的比例隨機(jī)劃分為建模集和驗(yàn)證集,在建模集和驗(yàn)證集中,以土壤鹽分含量為因變量,光譜指數(shù)為自變量,建立光譜指數(shù)組合與土壤鹽分的PLSR回歸模型

    模型穩(wěn)定程度和精度高低的評(píng)價(jià)指標(biāo)分別采用決定系數(shù) R2 和平均相對(duì)誤差 Emr 來(lái)評(píng)估, R2 越大、Emr 越小的模型的預(yù)測(cè)能力越強(qiáng),穩(wěn)定性越好。評(píng)價(jià)指標(biāo)數(shù)學(xué)表達(dá)式為

    式中: Si?Mi 分別為采樣點(diǎn) i 的預(yù)測(cè)值和實(shí)測(cè)值; n 為采集土壤樣點(diǎn)的數(shù)量; 為所有采樣點(diǎn)的實(shí)測(cè)值的平均值。

    2 結(jié)果與分析

    2.1多光譜指數(shù)與土壤鹽分的地理探測(cè)器分析

    表3為土壤鹽分與光譜指數(shù)的因子探測(cè)結(jié)果從表中可以看出, q 大于0.5的因子依次為可見(jiàn)-紅外鹽分指數(shù) 、擴(kuò)展差值植被指數(shù) IEDV 、歸一化差異水體指數(shù) INDW 、非線性指數(shù) INL 、衛(wèi)星 B8 波段的反射率 λB8 、S6鹽分指數(shù) Is6 、衛(wèi)星B6波段的反射率 λB6 和亮度指數(shù) IB 。S2鹽分指數(shù) Is2 、SI11鹽分指數(shù) IsI11 、綠波段比值植被指數(shù) IGRV 、衛(wèi)星B1 波段的反射率 λBl 的 q 均小于0.05,表明這些光譜指數(shù)對(duì)黃河三角洲引黃灌區(qū)土壤鹽分雖然具體一定影響,但相比于其他光譜指數(shù)影響較為微弱。本文中選取 q 大于0.05的因子用于構(gòu)建土壤鹽分反演模型。

    表3黃河三角洲引黃灌區(qū)土壤鹽分與光譜指數(shù)的因子探測(cè)結(jié)果
    續(xù)表

    2.2 模型構(gòu)建及精度分析

    利用統(tǒng)計(jì)與數(shù)據(jù)分析軟件Minitab19對(duì)研究區(qū)土壤鹽分以及 q 大于0.05的光譜指數(shù)進(jìn)行PLSR建模,在建模過(guò)程中,剔除 p 值大于0.01的光譜指數(shù),以確保只保留對(duì)土壤鹽分有顯著影響的變量,從而提高模型的擬合效果和預(yù)測(cè)準(zhǔn)確性。最終優(yōu)選出最佳的光譜指數(shù)組合為 IIS-VIR 、 IEDV INDW ! INL ) Is6 、IB 、衛(wèi)星B8A 波段的反射率 λB8A 、 INT2 鹽分指數(shù)IINT2 、優(yōu)化土壤調(diào)整植被指數(shù) IosAV 、土壤調(diào)節(jié)植被指數(shù) ISRV 、增強(qiáng)型植被指數(shù) IEV 、SI3 鹽分指數(shù) ISI3 、擴(kuò)展比值植被指數(shù) IERV 、鹽分指數(shù) Is 、大氣阻抗植被指數(shù) IARV 、S5鹽分指數(shù) Is5 ?;趯?shí)測(cè)土壤鹽分?jǐn)?shù)據(jù)和最優(yōu)光譜指數(shù)組合,構(gòu)建PLSR土壤鹽分反演模型為

    y=-3.53-27.21IIS-VIR+2.30IEDV-3.73INDW-3.58INL-

    39.73Is6-0.39IB-2.37λBsA+20.44IINT2+12.61IoSAV+ 24.14IsRv-10.01IEV-20.39IsI3-1.13IERV+16.42Is- 8.36IARV+37.56Is5 (4)

    式中 y 為土壤鹽分含量。

    通過(guò)評(píng)價(jià)指標(biāo)對(duì)建模集土壤鹽分預(yù)測(cè)值與實(shí)測(cè)值進(jìn)行對(duì)比,所構(gòu)建的模型 R2 為0.79, Emr 為0.28,表明模型預(yù)測(cè)效果較好,準(zhǔn)確度較高,對(duì)研究區(qū)有較好的適用性。

    對(duì)最優(yōu)光譜指數(shù)組合所構(gòu)建的土壤鹽分遙感反演模型進(jìn)行驗(yàn)證,結(jié)果如圖2所示。由圖可知,土壤鹽分預(yù)測(cè)值和實(shí)測(cè)值之間具有較強(qiáng)的線性相關(guān),并且通過(guò)評(píng)價(jià)指標(biāo)對(duì)驗(yàn)證集土壤鹽分預(yù)測(cè)值與實(shí)測(cè)值進(jìn)行評(píng)價(jià),得出驗(yàn)證集 R2 為0.77, Emr 為0.36,表明本文中構(gòu)建的土壤鹽分遙感反演模型可以有效地預(yù)測(cè)研究灌區(qū)的土壤鹽漬化狀況。

    圖2黃河三角洲引黃灌區(qū)土壤鹽分含量實(shí)測(cè)值和預(yù)測(cè)值比較

    2.3 土壤鹽漬化空間分布分析

    根據(jù)最優(yōu)土壤鹽分遙感反演模型預(yù)測(cè)研究區(qū)土壤鹽分含量空間分布結(jié)果如圖3所示。由圖可見(jiàn),研究區(qū)土壤鹽漬化程度總體較輕,土壤鹽分含量較低,鹽分質(zhì)量比均小于 6g/kg ,土壤鹽漬化等級(jí)包括非鹽漬土、輕度鹽漬土和中度鹽漬土,無(wú)重度鹽漬土和鹽土。大部分土壤的鹽分質(zhì)量比小于 2g/kg ,屬于非鹽漬土范疇。由于研究區(qū)大部分為作物種植區(qū),為了作物更好地生長(zhǎng)會(huì)進(jìn)行土壤淋洗和灌溉等操作,因此土壤鹽漬化程度較低。研究區(qū)內(nèi)不同等級(jí)土壤鹽漬化面積及所占比例統(tǒng)計(jì)結(jié)果如表4所示。從表中可以看出,非鹽漬土、輕度鹽漬土、中度鹽漬土面積分別為 ,非鹽漬土面積所占比例達(dá)到99. 66%

    圖3黃河三角洲引黃灌區(qū)土壤鹽分空間分布(資料來(lái)源:歐洲航天局哥白尼數(shù)據(jù)中心網(wǎng)站,經(jīng)ArcGIS10.8軟件數(shù)字化處理。)
    表4黃河三角洲引黃灌區(qū)土壤鹽漬化面積及所占比例

    為了將研究區(qū)土壤鹽分的空間分布特征更好的匹配各種土地利用-覆被類型,將非鹽漬土進(jìn)行細(xì)分,以非鹽漬土、輕度鹽漬土和中度鹽漬土劃分等級(jí)的 1/10 對(duì)研究區(qū)非鹽漬土等級(jí)進(jìn)行細(xì)分。就非鹽漬土而言,土壤鹽分質(zhì)量比小于 0.2,0.2~lt;0.4 0.4~lt;0.6,0.6~lt;2.0g/kg 的面積分別為94.82、54.30,21.19,17.08km2 。其中,鹽分質(zhì)量比小于0.2g/kg 鹽漬土的面積占比最大,為 50.43% ;其次是鹽分質(zhì)量比為 0.2~lt;0.4g/kg 的鹽漬土,面積占比為 28.88% ;鹽分質(zhì)量比 0.6~lt;2.0g/kg 鹽漬土面積所占比例最小,為 9.08% 。研究區(qū)內(nèi)輕度鹽漬土和中度鹽漬呈現(xiàn)零散的空間分布。因?yàn)檠芯繀^(qū)內(nèi)大片的田塊均進(jìn)行農(nóng)作物種植,所以在分散的無(wú)種植地帶會(huì)有較少的輕度和中度鹽漬土分布。

    在不同土地利用方式下,由于灌溉、入滲和蒸發(fā)等引起土壤水分在時(shí)間和空間上變化,因此導(dǎo)致區(qū)域土壤鹽分含量與運(yùn)移規(guī)律不同。為了更好地分析研究區(qū)土壤鹽分空間分布特征,本文中根據(jù)當(dāng)?shù)胤N植方式、種植結(jié)構(gòu)、休耕期等實(shí)際情況將研究區(qū)劃分為冬小麥-夏玉米地、擢荒地-夏玉米地、林地、鹽荒地和水稻田5種土地利用-覆被類型,探究土壤鹽分和不同土地利用-覆被類型之間的關(guān)系,結(jié)果如圖4所示,土地利用-覆被類型土壤鹽分含量統(tǒng)計(jì)結(jié)果如圖5所示。由圖可以看出:不同土地利用-覆被類型的土壤鹽分含量的均值差異顯著。鹽荒地土壤鹽分質(zhì)量比均值最大,為 0.37g/kg ,水稻田土壤鹽分質(zhì)量比均值最小,為 0.07g/kg 。土壤鹽分含量均值由大到小的土地利用-覆被類型的排序?yàn)辂}荒地、荒地-夏玉米地、林地、冬小麥-夏玉米地、水稻田,這主要是由土地利用方式、地表覆被類型、灌區(qū)種植結(jié)構(gòu)、灌溉制度所導(dǎo)致。研究區(qū)內(nèi)鹽荒地主要分布一些鹽生植物,相關(guān)研究[26]表明,鹽生植物的脫鹽功能能改良鹽漬化土壤,原因是不同鹽生植物對(duì)鹽分吸收具有選擇性,有些鹽生植物甚至可以將鹽分直接移出土壤,例如在研究區(qū)鹽荒地普遍生長(zhǎng)的拒鹽植物——蘆葦以及聚鹽植物—一鹽地堿蓬等,因此研究區(qū)內(nèi)鹽荒地的土壤鹽分含量較其他土地利用-覆被類型的高,但仍低于中度鹽漬土的鹽分含量。由于擢荒地-夏玉米地、林地、冬小麥-夏玉米地和水稻田長(zhǎng)期進(jìn)行種植管理,因此土壤鹽分含量相對(duì)較低。Wang等[27]研究表明,由于植被生長(zhǎng)隨季節(jié)發(fā)生變化,因此不同土地利用-覆被類型的土壤鹽分特征也會(huì)隨著季節(jié)性變化有所差異。本研究是在5月份對(duì)研究區(qū)進(jìn)行采樣及反演等工作,此時(shí)林地和冬小麥-夏玉米地均有較密集的植被覆被,且冬小麥和水稻均經(jīng)過(guò)了春灌等農(nóng)業(yè)處理。由于擢荒地-夏玉米地的土壤鹽分常年處于較高水平,不適宜冬小麥生長(zhǎng),在冬小麥種植期設(shè)置為擢荒期,于翌年直接種植夏玉米,因此在研究期內(nèi)畧荒地-夏玉米地作物處于玉米苗期或未種植期,導(dǎo)致土壤鹽分含量高于林地和冬小麥-夏玉米地。這種空間分布特征與實(shí)地調(diào)查和實(shí)驗(yàn)室分析的不同土地利用-覆被類型土壤鹽分含量相一致。

    圖4黃河三角洲引黃灌區(qū)土地利用-覆被空間分布(資料來(lái)源:歐洲航天局哥白尼數(shù)據(jù)中心網(wǎng)站,經(jīng)ArcGIS10.8軟件數(shù)字化處理。)
    圖5 黃河三角洲引黃灌不同土地利用-覆被類型土壤鹽分含量統(tǒng)計(jì)結(jié)果

    3結(jié)論

    本文中采用野外取樣結(jié)合哨兵-2B衛(wèi)星遙感影像獲取土壤鹽分?jǐn)?shù)據(jù),計(jì)算鹽分指數(shù)、植被指數(shù)并通過(guò)地理探測(cè)器篩選光譜指數(shù)組合,構(gòu)建PLSR土壤鹽分遙感反演模型,監(jiān)測(cè)研究區(qū)土壤鹽分含量,探討土壤鹽分的空間分布特征,得到如下主要結(jié)論:

    1)優(yōu)選出對(duì)土壤鹽分解釋力較高( qgt;0. 05) 且達(dá)到極顯著水平( p 值小于0.01)的最優(yōu)光譜指數(shù)組合是 (20osAv?IsRv?IEv?IsI3?IERv?Is?IARv?Is? 基于PLSR構(gòu)建的土壤鹽分反演模型精度較高,其中驗(yàn)證 R2 為0.77,Emr 為0.36。

    2)研究區(qū)土壤鹽漬化等級(jí)包括非鹽漬土、輕度鹽漬土和中度鹽漬土,其中以非鹽漬土為主,在非鹽漬土中鹽分質(zhì)量比小于 0.2g/kg 的面積占比最大,為 50.43% 。輕度鹽漬土和中度鹽漬土在空間上呈現(xiàn)零散分布,主要分布在非種植區(qū)域。

    3)研究區(qū)不同土地利用-覆被類型的土壤鹽分含量均值差異顯著,鹽分質(zhì)量比由大到小的土地利用-覆被類型排序?yàn)辂}荒地、荒地-夏玉米地、林地、冬小麥-夏玉米地、水稻田,其中鹽荒地土壤鹽分質(zhì)量比均值為 0.37g/kg ,水稻田土壤鹽分均值最小為 0.07g/kg

    本文中的研究成果可為濱海土壤鹽漬化快速監(jiān)測(cè)以及土地利用結(jié)構(gòu)空間優(yōu)化與可持續(xù)發(fā)展提供技術(shù)支撐。

    參考文獻(xiàn):

    [1]ZHUW,YANGJS,YAORJ,etal.Soilwater-saltcontrol andyield improvement under the effect of compound control insalinesoilof the YellowRiverDelta,China[J].Agricultural WaterManagement,2022,263:107455.

    [2] 楊小虎,羅艷琴,楊海昌,等.瑪納斯河流域綠洲農(nóng)田土壤鹽分反演及空間分布特征[J].干旱區(qū)資源與環(huán)境,2021,35(2):156.

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

    [4]周曉紅,張飛,張海威,等.艾比湖濕地自然保護(hù)區(qū)土壤鹽分多光譜遙感反演模型[J].光譜學(xué)與光譜分析,2019,39(4):1229.

    [5]張同瑞,趙庚星,高明秀,等.基于近地面多光譜的黃河三角洲典型地區(qū)土壤含鹽量估算研究[J].光譜學(xué)與光譜分析,2016, 36(1) : 248.

    [6]張素銘,趙庚星,王卓然,等.濱海鹽漬區(qū)土壤鹽分遙感反演及動(dòng)態(tài)監(jiān)測(cè)[J].農(nóng)業(yè)資源與環(huán)境學(xué)報(bào),2018,35(4):349.

    [7]劉恩,王軍濤,常步輝,等.小開(kāi)河引黃灌區(qū)土壤鹽漬化定量遙感反演[J].中國(guó)農(nóng)村水利水電,2019(12):20.

    [8]劉偉,史海濱,苗慶豐,等.基于Landsat8的節(jié)水改造背景下鹽漬化土壤含鹽量反演[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2024,55(1):294.

    [9]WANG JZ,DING JL,YU D L,et al. Machine learning-baseddetection of soil salinity in an arid desert region,Northwest China:a comparison between Landsat-8 OLI and Sentinel-2 MSI[J].Science of the Total Environment,2020,707:136092.

    [10]袁瑞強(qiáng),劉貫群,宋獻(xiàn)方.現(xiàn)代黃河三角洲淺層地下水對(duì)降水的響應(yīng)[J].資源科學(xué),2009,31(9):1514.

    [11]ZHANG JR, ZHANG Z T,CHEN JY,et al. Estimating soilsalinity with different fractional vegetation cover using remotesensing[J]. Land Degradation and Development, 2020,32(2):597.

    [12]姚志華,陳俊英,張智韜,等.覆膜對(duì)無(wú)人機(jī)多光譜遙感反演土壤含鹽量精度的影響[J].農(nóng)業(yè)工程學(xué)報(bào),2019,35(19): 89.

    [13]DOUAOUI A E K,NICOLAS H, WALTER C. Detecting salinityhazards within a semiarid context by means of combining soil andremote-sensing data[J].Geoderma,2006,134(1/2):217.

    [14]TRIKI FOURATI H, BOUAZIZ M,BENZINA M,et al. Modelingof soil salinity withinasemi-arid region using spectral aalysis[J].Arabian Journal of Geosciences,2015,8(12): 11175.

    [15]SCUDIERO E, SKAGGS T H, CORWIN D L. Regional scalesoil salinity evaluation using Landsat 7, western San Joaquin Valley,California,USA[J]. Geoderma Regional,2014,2-3:82.

    [16]張智韜,魏廣飛,姚志華,等.基于無(wú)人機(jī)多光譜遙感的土壤含鹽量反演模型研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2019,50(12):151.

    [17]CHO K H,BEON MS,JEONG JC. Dynamics of soil salinityand vegetation in a reclaimed area in Saemangeum, Republic ofKorea[J]. Geoderma,2018,321: 42.

    [18]馮娟,丁建麗,楊愛(ài)霞,等.干旱區(qū)土壤鹽漬化信息遙感建模[J].干旱地區(qū)農(nóng)業(yè)研究,2018,36(1):266.

    [19]毛海穎,馮仲科,鞏垠熙,等.多光譜遙感技術(shù)結(jié)合遺傳算法對(duì)永定河土壤歸一化水體指數(shù)的研究[J].光譜學(xué)與光譜分析,2014,34(6):1649.

    [20]奚雪,趙庚星,高鵬,等.基于Sentinel衛(wèi)星及無(wú)人機(jī)多光譜的濱海冬小麥種植區(qū)土壤鹽分反演研究:以黃三角墾利區(qū)為例[J].中國(guó)農(nóng)業(yè)科學(xué),2020,53(24):5005.

    猜你喜歡
    模型研究
    一半模型
    FMS與YBT相關(guān)性的實(shí)證研究
    2020年國(guó)內(nèi)翻譯研究述評(píng)
    遼代千人邑研究述論
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    視錯(cuò)覺(jué)在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    EMA伺服控制系統(tǒng)研究
    新版C-NCAP側(cè)面碰撞假人損傷研究
    3D打印中的模型分割與打包
    夜夜夜夜夜久久久久| 伦理电影免费视频| 高清av免费在线| 男女床上黄色一级片免费看| 国产三级黄色录像| 在线看a的网站| 汤姆久久久久久久影院中文字幕| 亚洲av男天堂| 啦啦啦在线免费观看视频4| 国产精品熟女久久久久浪| 久久久久网色| 精品人妻一区二区三区麻豆| 五月开心婷婷网| 国产精品 国内视频| 久久久国产成人免费| 色综合欧美亚洲国产小说| 青青草视频在线视频观看| 精品国产国语对白av| 国产日韩欧美视频二区| 一区福利在线观看| 久久中文字幕一级| 亚洲欧美成人综合另类久久久| 欧美黄色淫秽网站| 一级毛片精品| 免费日韩欧美在线观看| 纯流量卡能插随身wifi吗| 国产真人三级小视频在线观看| 人妻一区二区av| 少妇裸体淫交视频免费看高清 | 爱豆传媒免费全集在线观看| 黄色a级毛片大全视频| 性少妇av在线| 国产精品久久久久久人妻精品电影 | 丝袜美足系列| 两个人看的免费小视频| 精品国产乱码久久久久久小说| 国产精品九九99| 999久久久国产精品视频| 久久久久久久精品精品| 午夜老司机福利片| 国产黄频视频在线观看| 美女高潮到喷水免费观看| 97在线人人人人妻| 超碰97精品在线观看| 9191精品国产免费久久| 爱豆传媒免费全集在线观看| 国产成人av教育| 久久精品亚洲av国产电影网| 久久天堂一区二区三区四区| 精品少妇一区二区三区视频日本电影| 亚洲国产精品一区三区| 久久精品亚洲av国产电影网| 久久中文字幕一级| 91精品三级在线观看| 亚洲国产毛片av蜜桃av| 丝袜在线中文字幕| 日韩 亚洲 欧美在线| 人妻人人澡人人爽人人| xxxhd国产人妻xxx| 两人在一起打扑克的视频| 五月开心婷婷网| av在线老鸭窝| 激情视频va一区二区三区| 久久影院123| 免费在线观看黄色视频的| 少妇裸体淫交视频免费看高清 | 人妻人人澡人人爽人人| 欧美一级毛片孕妇| 99国产精品一区二区蜜桃av | 中国国产av一级| 制服诱惑二区| 欧美黄色片欧美黄色片| 国产老妇伦熟女老妇高清| 欧美少妇被猛烈插入视频| 亚洲av欧美aⅴ国产| 97在线人人人人妻| 精品福利观看| 久久国产精品人妻蜜桃| 女人久久www免费人成看片| 婷婷成人精品国产| 日韩制服丝袜自拍偷拍| 久久久久久免费高清国产稀缺| 热99re8久久精品国产| 动漫黄色视频在线观看| av网站在线播放免费| 超色免费av| 18禁国产床啪视频网站| 老司机靠b影院| 大片免费播放器 马上看| 日本av免费视频播放| www.999成人在线观看| 精品久久久精品久久久| 日本av免费视频播放| 日韩欧美一区视频在线观看| 亚洲专区字幕在线| 别揉我奶头~嗯~啊~动态视频 | 久久久久久久久久久久大奶| 欧美日韩亚洲高清精品| 欧美少妇被猛烈插入视频| 免费高清在线观看日韩| 国产精品一区二区精品视频观看| av超薄肉色丝袜交足视频| 乱人伦中国视频| 亚洲欧美日韩高清在线视频 | 国产日韩欧美在线精品| 成人影院久久| 欧美激情极品国产一区二区三区| 一个人免费看片子| tube8黄色片| 黄色视频在线播放观看不卡| xxxhd国产人妻xxx| 狂野欧美激情性xxxx| 69精品国产乱码久久久| 视频在线观看一区二区三区| 国产亚洲精品久久久久5区| 日韩三级视频一区二区三区| 国产av一区二区精品久久| 久久久久久人人人人人| av免费在线观看网站| 久久 成人 亚洲| 人妻久久中文字幕网| 老司机午夜十八禁免费视频| 99国产精品一区二区蜜桃av | 国产成人免费无遮挡视频| 一区二区三区乱码不卡18| 下体分泌物呈黄色| 高清av免费在线| 国产成人精品在线电影| 每晚都被弄得嗷嗷叫到高潮| 制服人妻中文乱码| 性色av乱码一区二区三区2| 国产成人欧美| 一级毛片电影观看| 天天影视国产精品| 高清欧美精品videossex| 老司机午夜十八禁免费视频| 欧美黄色片欧美黄色片| 777久久人妻少妇嫩草av网站| 欧美一级毛片孕妇| 国产成人精品在线电影| 国产精品一二三区在线看| 一级片'在线观看视频| 亚洲国产精品一区二区三区在线| 久久久久久久国产电影| 黄色a级毛片大全视频| 日本vs欧美在线观看视频| 十八禁人妻一区二区| 欧美日韩成人在线一区二区| 亚洲天堂av无毛| 日本vs欧美在线观看视频| 桃花免费在线播放| 欧美另类亚洲清纯唯美| 精品一区在线观看国产| 啦啦啦在线免费观看视频4| 美国免费a级毛片| 免费在线观看日本一区| 色婷婷av一区二区三区视频| 大香蕉久久成人网| 最近中文字幕2019免费版| 深夜精品福利| 欧美97在线视频| 12—13女人毛片做爰片一| 国产精品久久久久久精品电影小说| 免费日韩欧美在线观看| 91精品伊人久久大香线蕉| 热re99久久国产66热| 大片免费播放器 马上看| 大片免费播放器 马上看| 国产成人精品久久二区二区免费| 9191精品国产免费久久| 丝袜美足系列| 久久午夜综合久久蜜桃| 一级毛片女人18水好多| xxxhd国产人妻xxx| www.熟女人妻精品国产| 中文字幕最新亚洲高清| 久久午夜综合久久蜜桃| 人人澡人人妻人| 国产99久久九九免费精品| 丁香六月欧美| 91精品伊人久久大香线蕉| xxxhd国产人妻xxx| 99久久99久久久精品蜜桃| 亚洲av成人一区二区三| 五月天丁香电影| 国产欧美日韩一区二区三区在线| 在线亚洲精品国产二区图片欧美| 午夜精品国产一区二区电影| 不卡一级毛片| 999久久久国产精品视频| www.自偷自拍.com| 久久久国产欧美日韩av| 日本欧美视频一区| 欧美黄色片欧美黄色片| 岛国毛片在线播放| 久久精品人人爽人人爽视色| 日韩视频一区二区在线观看| 国产免费福利视频在线观看| 亚洲av国产av综合av卡| 制服人妻中文乱码| √禁漫天堂资源中文www| 搡老乐熟女国产| 99久久综合免费| 一本久久精品| a 毛片基地| 妹子高潮喷水视频| 午夜精品国产一区二区电影| www.自偷自拍.com| 欧美97在线视频| 亚洲成人国产一区在线观看| 青春草视频在线免费观看| 丝袜喷水一区| 国产精品香港三级国产av潘金莲| 另类亚洲欧美激情| 日韩熟女老妇一区二区性免费视频| 男女边摸边吃奶| 国产麻豆69| 淫妇啪啪啪对白视频 | 亚洲熟女毛片儿| 久久久久久人人人人人| 国产野战对白在线观看| 精品亚洲成国产av| 久久天堂一区二区三区四区| 国产有黄有色有爽视频| 女人精品久久久久毛片| 精品卡一卡二卡四卡免费| 妹子高潮喷水视频| 精品视频人人做人人爽| 国产一区二区三区在线臀色熟女 | 国产精品久久久久久人妻精品电影 | 日本黄色日本黄色录像| 亚洲五月婷婷丁香| 搡老乐熟女国产| 热99国产精品久久久久久7| 黄色 视频免费看| 日韩视频在线欧美| 日韩,欧美,国产一区二区三区| 亚洲精品久久久久久婷婷小说| 亚洲人成电影免费在线| 日本wwww免费看| 国产一区二区 视频在线| 老司机影院成人| 韩国精品一区二区三区| 欧美黑人欧美精品刺激| 亚洲精华国产精华精| 国产成人av教育| av网站在线播放免费| 午夜成年电影在线免费观看| 高潮久久久久久久久久久不卡| 国产老妇伦熟女老妇高清| 国产黄频视频在线观看| 亚洲黑人精品在线| 免费人妻精品一区二区三区视频| 丰满迷人的少妇在线观看| 精品国产乱子伦一区二区三区 | 在线观看免费高清a一片| 18禁黄网站禁片午夜丰满| 日本av免费视频播放| 欧美老熟妇乱子伦牲交| 成年人免费黄色播放视频| 男人添女人高潮全过程视频| 少妇精品久久久久久久| 精品人妻一区二区三区麻豆| 日韩电影二区| 日本黄色日本黄色录像| 捣出白浆h1v1| videosex国产| 亚洲国产精品一区二区三区在线| 欧美日韩福利视频一区二区| 亚洲国产中文字幕在线视频| 中文字幕色久视频| www.熟女人妻精品国产| 国产精品1区2区在线观看. | 色老头精品视频在线观看| 老司机深夜福利视频在线观看 | 啦啦啦免费观看视频1| 精品一区二区三卡| 丝袜喷水一区| 天堂俺去俺来也www色官网| 在线观看www视频免费| 一二三四在线观看免费中文在| 曰老女人黄片| 啦啦啦啦在线视频资源| 老熟女久久久| 国产成人系列免费观看| 亚洲精品久久成人aⅴ小说| 中文字幕人妻丝袜一区二区| 美女高潮到喷水免费观看| 欧美日本中文国产一区发布| 日韩,欧美,国产一区二区三区| 久久中文字幕一级| 十八禁人妻一区二区| 中文字幕人妻熟女乱码| 国产亚洲欧美精品永久| 正在播放国产对白刺激| 丰满人妻熟妇乱又伦精品不卡| 日本一区二区免费在线视频| 一个人免费在线观看的高清视频 | 一本一本久久a久久精品综合妖精| 日韩欧美一区视频在线观看| 亚洲第一av免费看| 在线观看免费视频网站a站| 深夜精品福利| 欧美在线一区亚洲| 午夜老司机福利片| 欧美成狂野欧美在线观看| 欧美精品一区二区免费开放| 一本—道久久a久久精品蜜桃钙片| 亚洲av电影在线观看一区二区三区| 黄色怎么调成土黄色| 岛国在线观看网站| 日韩中文字幕欧美一区二区| 国产亚洲一区二区精品| 国产亚洲av高清不卡| 亚洲av成人不卡在线观看播放网 | 超色免费av| 日韩中文字幕视频在线看片| 青春草亚洲视频在线观看| 高潮久久久久久久久久久不卡| 啦啦啦视频在线资源免费观看| 夜夜骑夜夜射夜夜干| www.自偷自拍.com| 日韩欧美一区二区三区在线观看 | 伊人久久大香线蕉亚洲五| 免费黄频网站在线观看国产| www.av在线官网国产| 午夜福利免费观看在线| 狠狠狠狠99中文字幕| 国产成人影院久久av| 久久99一区二区三区| 日韩一卡2卡3卡4卡2021年| 成人影院久久| 国产精品久久久久久精品电影小说| 狠狠精品人妻久久久久久综合| 国产在线一区二区三区精| 一区二区日韩欧美中文字幕| 91精品三级在线观看| 亚洲av电影在线进入| 1024视频免费在线观看| 日本撒尿小便嘘嘘汇集6| 国产一区二区在线观看av| av福利片在线| 国产成人av教育| 女警被强在线播放| 亚洲免费av在线视频| √禁漫天堂资源中文www| 两个人看的免费小视频| 欧美国产精品一级二级三级| 黑人欧美特级aaaaaa片| 十分钟在线观看高清视频www| 日韩一卡2卡3卡4卡2021年| 成年人黄色毛片网站| 久久精品久久久久久噜噜老黄| 色综合欧美亚洲国产小说| 色视频在线一区二区三区| 99九九在线精品视频| 亚洲专区国产一区二区| 久久久国产欧美日韩av| 亚洲成人国产一区在线观看| 久久久国产精品麻豆| 国产成人精品无人区| 后天国语完整版免费观看| 久久精品成人免费网站| 99久久精品国产亚洲精品| 伦理电影免费视频| 高清欧美精品videossex| 黄色 视频免费看| 好男人电影高清在线观看| 亚洲精品久久成人aⅴ小说| 精品国产乱码久久久久久男人| 中国美女看黄片| e午夜精品久久久久久久| 我的亚洲天堂| xxxhd国产人妻xxx| 窝窝影院91人妻| 国产精品一二三区在线看| 女人被躁到高潮嗷嗷叫费观| 成年动漫av网址| 成人免费观看视频高清| 99国产精品一区二区蜜桃av | 搡老乐熟女国产| 中文字幕精品免费在线观看视频| 午夜视频精品福利| 亚洲精品一二三| 国产精品自产拍在线观看55亚洲 | 伊人亚洲综合成人网| 国产精品成人在线| 嫁个100分男人电影在线观看| av又黄又爽大尺度在线免费看| 中文字幕av电影在线播放| 婷婷色av中文字幕| 国产不卡av网站在线观看| 1024视频免费在线观看| 看免费av毛片| 色老头精品视频在线观看| 另类精品久久| 精品乱码久久久久久99久播| 国产精品av久久久久免费| 91国产中文字幕| 男女边摸边吃奶| 伊人亚洲综合成人网| 日韩中文字幕视频在线看片| 中文字幕最新亚洲高清| 亚洲va日本ⅴa欧美va伊人久久 | 老熟女久久久| www.999成人在线观看| 欧美日韩视频精品一区| 中文字幕色久视频| 我的亚洲天堂| 亚洲欧美激情在线| 国产免费av片在线观看野外av| 永久免费av网站大全| 制服诱惑二区| 91大片在线观看| 欧美成人午夜精品| 老汉色∧v一级毛片| av在线app专区| 久久av网站| 老汉色av国产亚洲站长工具| 交换朋友夫妻互换小说| 欧美精品一区二区免费开放| 午夜老司机福利片| av一本久久久久| 一本综合久久免费| 91字幕亚洲| 日韩欧美一区二区三区在线观看 | 国精品久久久久久国模美| 无限看片的www在线观看| 少妇精品久久久久久久| 日韩视频一区二区在线观看| 最黄视频免费看| 亚洲精品av麻豆狂野| 一二三四在线观看免费中文在| 亚洲成人免费电影在线观看| 欧美变态另类bdsm刘玥| 美女大奶头黄色视频| 可以免费在线观看a视频的电影网站| 国产av国产精品国产| 纯流量卡能插随身wifi吗| 亚洲中文av在线| 欧美精品一区二区大全| 热99re8久久精品国产| 黄色a级毛片大全视频| 欧美性长视频在线观看| av天堂在线播放| videosex国产| 中文字幕人妻熟女乱码| 精品少妇内射三级| 伊人久久大香线蕉亚洲五| 久久久久精品人妻al黑| 永久免费av网站大全| 亚洲第一欧美日韩一区二区三区 | 欧美av亚洲av综合av国产av| 精品久久久久久久毛片微露脸 | 日韩欧美免费精品| 国产一区二区三区av在线| 欧美日韩亚洲国产一区二区在线观看 | 国产极品粉嫩免费观看在线| 99香蕉大伊视频| 欧美日韩亚洲国产一区二区在线观看 | 日本五十路高清| 大香蕉久久网| 伊人亚洲综合成人网| 国产一区二区在线观看av| 我的亚洲天堂| 日韩一区二区三区影片| 欧美精品高潮呻吟av久久| 国产av精品麻豆| 国产熟女午夜一区二区三区| 大码成人一级视频| 亚洲自偷自拍图片 自拍| 欧美精品高潮呻吟av久久| 岛国毛片在线播放| 免费人妻精品一区二区三区视频| √禁漫天堂资源中文www| 老司机影院成人| 午夜视频精品福利| 国产又色又爽无遮挡免| cao死你这个sao货| 国产精品亚洲av一区麻豆| 一本大道久久a久久精品| videos熟女内射| 91成年电影在线观看| 国产极品粉嫩免费观看在线| 国产精品一区二区免费欧美 | 黄色视频在线播放观看不卡| 欧美变态另类bdsm刘玥| 美女大奶头黄色视频| 国产成人啪精品午夜网站| 九色亚洲精品在线播放| av在线老鸭窝| 国产成人av激情在线播放| 丰满少妇做爰视频| av天堂在线播放| 满18在线观看网站| 99香蕉大伊视频| 欧美中文综合在线视频| netflix在线观看网站| 久久精品国产a三级三级三级| 一区二区日韩欧美中文字幕| 俄罗斯特黄特色一大片| 狂野欧美激情性xxxx| 日韩 亚洲 欧美在线| 午夜免费观看性视频| 最近最新免费中文字幕在线| 高清在线国产一区| 午夜激情av网站| 国产淫语在线视频| 丝袜人妻中文字幕| 精品第一国产精品| 水蜜桃什么品种好| 国产精品香港三级国产av潘金莲| 满18在线观看网站| 高清av免费在线| 亚洲熟女毛片儿| 久久久精品免费免费高清| 欧美激情极品国产一区二区三区| 麻豆国产av国片精品| 91精品国产国语对白视频| 天天躁狠狠躁夜夜躁狠狠躁| 十八禁网站免费在线| 日本一区二区免费在线视频| xxxhd国产人妻xxx| 国产在线一区二区三区精| 99精国产麻豆久久婷婷| 亚洲人成电影免费在线| 亚洲av片天天在线观看| 日韩视频在线欧美| 久久九九热精品免费| 久久久久久久精品精品| 乱人伦中国视频| 永久免费av网站大全| 黑人猛操日本美女一级片| 日韩一区二区三区影片| 1024香蕉在线观看| 制服诱惑二区| 色婷婷久久久亚洲欧美| 亚洲国产精品999| 欧美人与性动交α欧美精品济南到| 又黄又粗又硬又大视频| 另类亚洲欧美激情| 成人国产av品久久久| 欧美xxⅹ黑人| www.熟女人妻精品国产| 手机成人av网站| 午夜福利,免费看| 性少妇av在线| 侵犯人妻中文字幕一二三四区| 精品国产乱码久久久久久男人| 高清在线国产一区| av在线app专区| 色精品久久人妻99蜜桃| 久久ye,这里只有精品| www.自偷自拍.com| 男女下面插进去视频免费观看| 亚洲精品国产一区二区精华液| 亚洲三区欧美一区| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲欧美清纯卡通| avwww免费| 后天国语完整版免费观看| a级毛片黄视频| 国产91精品成人一区二区三区 | 欧美黑人欧美精品刺激| 国产成人欧美在线观看 | 国产日韩欧美在线精品| 狠狠狠狠99中文字幕| 9色porny在线观看| 成人国产一区最新在线观看| 成人亚洲精品一区在线观看| 久久久久久久久久久久大奶| 首页视频小说图片口味搜索| 欧美成人午夜精品| 久久人妻福利社区极品人妻图片| 91老司机精品| 日韩电影二区| 精品国产乱码久久久久久小说| 国产亚洲av高清不卡| 叶爱在线成人免费视频播放| 美女扒开内裤让男人捅视频| 啦啦啦啦在线视频资源| 国产一区二区激情短视频 | 亚洲av国产av综合av卡| 精品国产乱子伦一区二区三区 | 亚洲精品久久成人aⅴ小说| 波多野结衣一区麻豆| 99九九在线精品视频| 欧美成人午夜精品| 精品久久久久久久毛片微露脸 | 亚洲色图 男人天堂 中文字幕| 精品高清国产在线一区| 亚洲精品中文字幕在线视频| 日本黄色日本黄色录像| 在线av久久热| 亚洲专区字幕在线| 午夜福利视频在线观看免费| 大型av网站在线播放| 黄色a级毛片大全视频| 97在线人人人人妻| 亚洲精华国产精华精| 日本av手机在线免费观看| 老熟女久久久| 亚洲av欧美aⅴ国产| 大片免费播放器 马上看| 91成人精品电影| 视频区欧美日本亚洲| 国产人伦9x9x在线观看| 国产精品久久久久久精品古装| 国产成人精品在线电影| √禁漫天堂资源中文www| 蜜桃国产av成人99| 亚洲七黄色美女视频| 色视频在线一区二区三区| 韩国精品一区二区三区| a在线观看视频网站|