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

    渝西北土壤有機(jī)質(zhì)空間變異及影響因素分析

    2017-12-01 09:47:24鮑麗然賈中民
    關(guān)鍵詞:插值坡度土地利用

    鮑麗然,周 皎,李 瑜,賈中民

    (重慶市地質(zhì)礦產(chǎn)勘查開發(fā)局川東南地質(zhì)大隊(duì),重慶 400038)

    渝西北土壤有機(jī)質(zhì)空間變異及影響因素分析

    鮑麗然,周 皎,李 瑜,賈中民

    (重慶市地質(zhì)礦產(chǎn)勘查開發(fā)局川東南地質(zhì)大隊(duì),重慶 400038)

    【目的】深入認(rèn)識(shí)重慶西北部丘陵區(qū)土壤有機(jī)質(zhì)空間變異特征及影響因素,為農(nóng)業(yè)生產(chǎn)、土地利用等提供科學(xué)依據(jù)?!痉椒ā坎捎玫亟y(tǒng)計(jì)學(xué)和GIS相結(jié)合的方法,進(jìn)行有機(jī)質(zhì)空間變異和基于土地利用類型修正的空間插值研究?!窘Y(jié)果】研究區(qū)土壤有機(jī)質(zhì)平均含量為16.26 g/kg,屬中等偏下水平;變異系數(shù)為33.8 %,存在中等程度的空間變異性;土壤有機(jī)質(zhì)空間分布呈現(xiàn)中間高,兩端低的趨勢(shì),合川區(qū)整體偏低。【結(jié)論】土壤有機(jī)質(zhì)具有中等空間自相關(guān)性,空間變異性受結(jié)構(gòu)因素(地形坡度、土壤類型)和隨機(jī)因素(土地利用類型、耕作措施等)影響,土地利用類型是主控因素。

    地統(tǒng)計(jì)學(xué);土壤有機(jī)質(zhì);空間變異;影響因素

    【研究意義】有機(jī)質(zhì)是土壤的重要組成部分,包括各種動(dòng)植物的殘?bào)w、微生物體及其分解和合成的各種有機(jī)質(zhì)。土壤有機(jī)質(zhì)在耕地養(yǎng)分管理、農(nóng)業(yè)發(fā)展等方面具有重要意義,是土壤學(xué)、農(nóng)學(xué)和地球化學(xué)研究的重點(diǎn)問題之一[1-2]。查明土壤有機(jī)質(zhì)的空間變異特征及其影響因素,有助于制定合理的作物種植規(guī)劃及農(nóng)田配方施肥,促進(jìn)耕地的高效利用,保障農(nóng)業(yè)可持續(xù)發(fā)展[3]。【前人研究進(jìn)展】土壤有機(jī)質(zhì)空間變異的研究始于20世紀(jì)70年代[4],隨著GIS技術(shù)的發(fā)展,80年代后多位學(xué)者相繼運(yùn)用地統(tǒng)計(jì)學(xué)結(jié)合GIS方法,探討了農(nóng)田土壤有機(jī)質(zhì)空間變異特征[5-7],如Moorman等[6]對(duì)包括土壤有機(jī)質(zhì)在內(nèi)的多種養(yǎng)分進(jìn)行了空間變異研究,Mishra等[7]曾針對(duì)紅壤地區(qū)的有機(jī)質(zhì)和土壤酸堿度展開空間變異特征分析。近年,國(guó)內(nèi)學(xué)者對(duì)土壤有機(jī)質(zhì)空間變異性在農(nóng)田尺度、區(qū)域尺度上展開了大量研究,研究區(qū)涉及平原區(qū)、干旱荒漠區(qū)、丘陵區(qū)、典型喀斯特區(qū)、高原流域、山區(qū)小流域、城市近郊等[8-13],主要探討了自然因素(地形地貌、氣候和土壤類型等)和人為因素(土地利用和耕作方式等)對(duì)有機(jī)質(zhì)空間變異的影響?!颈狙芯壳腥朦c(diǎn)】2003年以來(lái),重慶市多目標(biāo)區(qū)域地球化學(xué)調(diào)查獲得了大量高精度土壤有機(jī)質(zhì)數(shù)據(jù),為研究其空間變異特征提供了基礎(chǔ)。渝西北地區(qū)地勢(shì)平緩、土壤肥沃,耕地廣泛分布,農(nóng)業(yè)是該區(qū)最重要的經(jīng)濟(jì)支柱,因此土壤有機(jī)質(zhì)的分布特征至關(guān)重要?!緮M解決的關(guān)鍵問題】基于GIS操作平臺(tái),研究了重慶西部丘陵區(qū)土壤有機(jī)質(zhì)空間變異特征及影響因素,為全面、準(zhǔn)確掌握該區(qū)土壤有機(jī)質(zhì)空間變異規(guī)律,更好的進(jìn)行農(nóng)業(yè)生產(chǎn)活動(dòng)、合理規(guī)劃土地利用等提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)位于重慶市西北部,地理坐標(biāo)為東經(jīng)105°17′ ~ 106°39′,北緯29°15′ ~ 30°26′,總面積約7130 km2(圖1)。地勢(shì)起伏變化不大,海拔一般為200 ~ 500 m,廣大平緩淺丘遍布全區(qū)。研究區(qū)出露地層主要為侏羅系沙溪廟組、遂寧組和蓬萊鎮(zhèn)組,其次為三疊系須家河組,第四系零星分布于山間河谷地區(qū)。礦產(chǎn)有煤、天青石、石灰?guī)r、鹽和天然氣等。土地利用類型以耕地為主,占土地總面積的60.1 %。土壤類型主要為廣大丘陵區(qū)的紫色土和水稻土,其次為黃壤、石灰土。該區(qū)屬亞熱帶季風(fēng)氣候,氣候溫和,雨量充沛,年均降雨量1022 mm。區(qū)內(nèi)分布嘉陵江水系和沱江水系兩大長(zhǎng)江一級(jí)支流,嘉陵江水系包括渠江、瓊江和涪江等,沱江水系包括瀨溪河及其支流。糧食作物主要為水稻、玉米、油菜、大豆、薯類,建有大型果蔬基地,是重慶市重要的糧油果蔬生產(chǎn)基地。

    1.2 樣品采集與測(cè)試

    按照1件/km2的密度采集0 ~ 20 cm表層土壤,樣點(diǎn)布設(shè)以代表性為主要原則,同時(shí)兼顧均勻性和合理性,最大限度控制調(diào)查面積。采樣時(shí)除去地表雜物,樣品加工處理后重量不低于800 g。將大格內(nèi)(4 km2)樣品組合為一個(gè)分析樣,全區(qū)共采集表層土壤樣品7135件,組合分析樣1783件。

    樣品分析測(cè)試按照DD2005-01《多目標(biāo)區(qū)域地球化學(xué)調(diào)查規(guī)范》(1∶250000)[14]的相關(guān)要求,由安徽省地質(zhì)實(shí)驗(yàn)研究所完成。實(shí)驗(yàn)室人員采用重鉻酸鉀容量法分析檢測(cè),檢出限0.05 %,報(bào)出率100 %,合格率100 %。

    1.3 數(shù)據(jù)處理

    運(yùn)用SPSS19.0軟件對(duì)土壤有機(jī)質(zhì)進(jìn)行經(jīng)典描述性統(tǒng)計(jì)分析;依據(jù)峰度、偏度值和非參數(shù)檢驗(yàn)的單樣本K-S值判別統(tǒng)計(jì)數(shù)據(jù)是否符合正態(tài)分布;利用單因素方差分析F檢驗(yàn)和回歸分析,分別判斷有機(jī)質(zhì)含量差異顯著性和不同因素對(duì)土壤有機(jī)質(zhì)的影響程度。

    借助地統(tǒng)計(jì)學(xué)軟件GS+9.0進(jìn)行半變異函數(shù)計(jì)算和模型擬合,結(jié)合專業(yè)知識(shí)分析土壤有機(jī)質(zhì)的空間變異性及其影響因素。

    圖1 地理位置Fig.1 Geographical position of the study area

    指標(biāo)Index樣本Samples含量范圍Contentrange中位數(shù)Median均值Mean標(biāo)準(zhǔn)差Std.dev.變異系數(shù)c.v.偏度Skewness峰度KurtosisK-S檢驗(yàn)K-Stest土壤有機(jī)質(zhì)(g/kg)17836.36~63.7115.5216.265.5033.833.2528.810.0對(duì)數(shù)轉(zhuǎn)換0.3941.540.165

    基于Arcgis10.0的Geostatistic analysis模塊,采用Kriging法對(duì)土壤有機(jī)質(zhì)殘差進(jìn)行空間插值,運(yùn)用柵格計(jì)算將土地利用圖層和殘差預(yù)測(cè)圖層相加,得到土壤有機(jī)質(zhì)含量空間分布圖。

    2 結(jié)果與分析

    2.1 土壤有機(jī)質(zhì)經(jīng)典統(tǒng)計(jì)分析

    從一般的經(jīng)典統(tǒng)計(jì)分析結(jié)果(表1)看,土壤有機(jī)質(zhì)含量介于6.36 ~ 63.71 g/kg,平均值為16.26 g/kg,根據(jù)《土地質(zhì)量地球化學(xué)評(píng)價(jià)規(guī)范》屬較缺乏水平[15];變異系數(shù)33.8 %,揭示了有機(jī)質(zhì)含量的離散程度,屬中等變異強(qiáng)度[16]。峰度值為28.8,較陡峭,并且K-S檢驗(yàn)中P值為0,明顯不符合正態(tài)分布。經(jīng)對(duì)數(shù)轉(zhuǎn)換后,偏度和峰度值分別為0.394和1.54,再經(jīng)K-S檢驗(yàn)表明符合正態(tài)分布(Pgt;0.05),滿足地統(tǒng)計(jì)學(xué)分析的基本要求。

    2.2 土壤有機(jī)質(zhì)的空間變異結(jié)構(gòu)特征

    運(yùn)用GS+9.0軟件對(duì)土壤有機(jī)質(zhì)含量進(jìn)行各向同性變異函數(shù)擬合(表2),按照決定系數(shù)(R2)最大、殘差(RSS)最小的原則選擇最優(yōu)半方差函數(shù)擬合模型,最佳理論模型為指數(shù)模型。

    空間變異主要包括隨機(jī)性和結(jié)構(gòu)性兩部分變異。塊金值C0表示區(qū)域化變量?jī)?nèi)部隨機(jī)性的可能程度,由測(cè)量誤差和最小取樣間距內(nèi)土壤性質(zhì)的人為隨機(jī)因子(施肥、耕作措施、土地利用方式等)引起。結(jié)構(gòu)方差C是由自然因素(地質(zhì)背景、地形、土壤類型等)引起土壤性狀的結(jié)構(gòu)變異[17]?;_(tái)值C0+C通常表示系統(tǒng)總變異,包括隨機(jī)變異和結(jié)構(gòu)變異。塊金系數(shù),即塊金值和基臺(tái)值之比C0/C0+C,表示空間異質(zhì)性程度,比值高表明隨機(jī)因素引起的空間變異程度較大;相反,則結(jié)構(gòu)性因素引起的空間變異程度較大[18]。C0/C0+Clt;25 %,說明系統(tǒng)具有強(qiáng)烈的空間相關(guān)性;25 %~ 75 %,說明系統(tǒng)具有中等空間相關(guān)性;C0/C0+Cgt;75 %,說明系統(tǒng)空間相關(guān)性弱[19]。變程反映在一定觀測(cè)尺度下,變量的空間自相關(guān)距離,只有變程以內(nèi)的空間變量才具有空間自相關(guān)性。

    研究區(qū)土壤有機(jī)質(zhì)的塊金系數(shù)C0/C0+C為46.2 %,表現(xiàn)出中等空間相關(guān)性,說明結(jié)構(gòu)性因素和人為隨機(jī)因素土壤有機(jī)質(zhì)的空間變異影響作用相當(dāng)。變程為184 km,說明土壤有機(jī)質(zhì)在此空間范圍內(nèi)存在空間自相關(guān)性。

    2.3 土壤有機(jī)質(zhì)空間變異影響因素分析

    2.3.1 土地利用類型 不同土地利用類型下,人類對(duì)土地的使用、保護(hù)和改造活動(dòng)不同,土壤有機(jī)質(zhì)含量隨之產(chǎn)生了較大差別[20]。研究區(qū)土地利用主要為旱地和水田,其次為林地。林地土壤有機(jī)質(zhì)平均含量最高(表3),為23.82 g/kg,水田為17.15 g/kg,旱地只有14.72 g/kg,相當(dāng)于水田的85 %;旱地的變異系數(shù)為0.35,是水田的1.24倍。研究區(qū)水田更易于土壤有機(jī)質(zhì)的累積,且含量較穩(wěn)定,與許泉[21]等發(fā)現(xiàn)的不同土地利用類型下農(nóng)田土壤有機(jī)碳密度區(qū)域差異規(guī)律相似。方差分析(表4)顯示,組間均方大于組內(nèi)均方,表明組間差異大于組內(nèi)差異,F(xiàn)檢驗(yàn)結(jié)果(Plt;0.05)說明不同土地利用類型土壤有機(jī)質(zhì)含量具有顯著性差異。

    表2 土壤有機(jī)質(zhì)空間變異的理論模型和相應(yīng)參數(shù)

    表3 不同土地利用類型土壤有機(jī)質(zhì)含量

    表4 不同土地利用類型土壤有機(jī)質(zhì)含量方差分析

    注:F檢驗(yàn)Plt;0.05。

    2.3.2 地形(海拔、坡度) 海拔、坡度等地形因子也會(huì)影響土壤有機(jī)質(zhì)含量。高海拔區(qū),人類活動(dòng)較少,易于積累土壤有機(jī)質(zhì)。坡度較大地區(qū),強(qiáng)烈的水土流失會(huì)攜帶走大量的土壤有機(jī)質(zhì)。

    將研究區(qū)1073個(gè)樣點(diǎn)有機(jī)質(zhì)含量和海拔高度進(jìn)行相關(guān)分析,得出相關(guān)系數(shù)為0.170,兩者關(guān)系極弱,認(rèn)為不相關(guān)[22]。進(jìn)一步統(tǒng)計(jì)了不同海拔高度上各等級(jí)有機(jī)質(zhì)含量百分比。由圖2可知,不同海拔高程帶上土壤有機(jī)質(zhì)均以10 ~ 20 g/kg為主,占所在高程帶的80 %左右,其次為20 ~ 30 g/kg。隨著海拔高度的變化,土壤有機(jī)質(zhì)沒有表現(xiàn)出較大的差異性。本區(qū)屬低緩丘陵區(qū),不存在高山或洼地,地勢(shì)起伏、氣候差異較小,土地利用、人類耕作活動(dòng)少受海拔影響。因此,海拔高度對(duì)土壤有機(jī)質(zhì)的空間變異影響不明顯。

    從地形的坡度因素考慮,根據(jù)各采樣點(diǎn)坡度信息,得到不同坡度范圍土壤有機(jī)質(zhì)均值(表5),土壤有機(jī)質(zhì)的空間分布和坡度密切相關(guān)。坡度小于5°時(shí),有機(jī)質(zhì)含量最高,為19.77 g/kg;大于35°時(shí),有機(jī)質(zhì)含量最低,為12.23 g/kg;隨著坡度增大,土壤有機(jī)質(zhì)含量呈逐漸降低趨勢(shì)。這與李婷等[23]對(duì)沱江流域中游土壤有機(jī)質(zhì)含量與坡度關(guān)系研究結(jié)果相符,不同坡度土地利用類型和土地管理水平不同,低緩坡度土地利用主要為水田,比較高陡的坡度土地利用為旱地和林地。另一方面研究區(qū)屬紫色丘陵區(qū),坡度大的區(qū)域水土流失較為嚴(yán)重,大量土壤有機(jī)質(zhì)被攜帶流出。

    圖2 不同海拔土壤有機(jī)質(zhì)分布特征Fig.2 Distributions of soil organic matter in different altitudes

    2.3.3 土壤類型 土壤類型也是有機(jī)質(zhì)含量的重要影響因素之一。研究區(qū)土壤類型以紫色土和水稻土為主,分別占總面積的70 %和21.6 %,黃壤和石灰土分布較少。紫色土主要分布于北部潼南、合川、大足地區(qū),水稻土于銅梁、榮昌地區(qū)分布較集中。統(tǒng)計(jì)分析(表6)表明,土壤有機(jī)質(zhì)平均含量表現(xiàn)為石灰土gt;黃壤gt;水稻土gt;紫色土,最高值石灰土是最低值紫色土的1.79倍,水稻土和黃壤差異較小。這與成都平原不同類型表層土壤中有機(jī)質(zhì)分布規(guī)律相同[24],均為石灰土、黃壤中較高,紫色土最低。由表7可知,不同土壤類型間,組間均方大于組內(nèi)均方,組間差異較大,F(xiàn)檢驗(yàn)結(jié)果(Plt;0.05)說明不同土壤類型中有機(jī)質(zhì)含量具有顯著性差異。

    表5 不同坡度土壤有機(jī)質(zhì)含量

    表6 不同土類有機(jī)質(zhì)含量組成

    表7 不同土類土壤有機(jī)質(zhì)含量方差分析

    注:F檢驗(yàn)Plt;0.05。

    表8 不同影響因子與土壤有機(jī)質(zhì)相關(guān)性

    2.4 土壤有機(jī)質(zhì)的空間分布格局

    土壤有機(jī)質(zhì)空間分布格局可借助空間插值模型實(shí)現(xiàn)從離散的點(diǎn)狀信息向連續(xù)面狀信息的轉(zhuǎn)換,用有限的樣點(diǎn)預(yù)測(cè)整個(gè)研究區(qū)域的土壤屬性空間分布特征。在空間尺度較小、地勢(shì)平坦且土地利用類型單一的區(qū)域,通常用普通克里格插值進(jìn)行預(yù)測(cè)[25],而研究區(qū)空間尺度較大,地形條件復(fù)雜、土地利用類型多樣,應(yīng)利用土壤有機(jī)質(zhì)相關(guān)因子與克里格插值相結(jié)合的方法來(lái)展示土壤有機(jī)質(zhì)的空間分布。目前,結(jié)合相關(guān)因子的空間插值的研究尚不多見。劉靜等[26]在研究梁平地區(qū)土壤有機(jī)質(zhì)空間分布時(shí),結(jié)合地形因素改進(jìn)了空間插值方法。顧成軍、文雯等[27-29]探討了基于土地利用類型修正的殘差法空間插值問題。

    本文首先通過回歸分析定量判斷了土地利用、土壤類型、地形坡度等因素對(duì)土壤有機(jī)質(zhì)空間變異的影響程度(表8)。研究區(qū)土壤有機(jī)質(zhì)含量與不同影響因素的相關(guān)系數(shù)大小依次為土地利用類型gt;土壤類型gt;地形(坡度),土地利用類型是該區(qū)土壤有機(jī)質(zhì)分布最為重要的影響因素。因此,為提高土壤有機(jī)質(zhì)空間分布預(yù)測(cè)精度,本文選擇結(jié)合土地利用類型的克里格插值法展開研究。首先將樣本土壤有機(jī)質(zhì)含量去除土地利用類型均值,得到殘差,并將其作為一個(gè)新的區(qū)域變量行普通克里格插值。然后,將插值過后的殘差值圖層按土地利用類型的不同,分別加上其地類均值,即得出土壤有機(jī)質(zhì)密度最終預(yù)測(cè)值,由此生成該區(qū)土壤有機(jī)質(zhì)含量空間分布圖(圖3)。

    研究區(qū)土壤有機(jī)質(zhì)含量插值結(jié)果范圍為1.96

    圖3 研究區(qū)土壤有機(jī)質(zhì)空間分布Fig.3 Spatial distribution of soil organic matter in the study area

    ~ 71.21 g/kg,空間分布上整體呈中間高,南北兩端低的趨勢(shì)。土壤有機(jī)質(zhì)高值區(qū)主要分布于大足區(qū)東部、銅梁區(qū)西部、潼南區(qū)南部和榮昌區(qū)最南端等地區(qū)。大足區(qū)境內(nèi)高值區(qū)占總面的60 %左右。這些高值區(qū)土地利用類型以水田為主,并有林地分布。合川區(qū)土壤有機(jī)質(zhì)整體偏低,榮昌區(qū)北部含量最低。合川區(qū)和潼南區(qū)北部分布嘉陵江及其支流渠江、涪江,水網(wǎng)交錯(cuò),水土流失相對(duì)嚴(yán)重,這可能是土壤有機(jī)質(zhì)低的主要因素。

    3 討 論

    研究區(qū)土壤有機(jī)質(zhì)的空間變異主要受土地利用、土壤類型的影響,不同土地利用或土壤類型的有機(jī)質(zhì)含量具有顯著差異。

    從土地利用類型看,林地區(qū)植被遍布,源源不斷的樹木殘枝落葉進(jìn)入土壤,有機(jī)質(zhì)來(lái)源充足,因此有機(jī)質(zhì)含量最豐富。耕地中水田土壤有機(jī)質(zhì)高于旱地,主要是由于水田區(qū)海拔較低,地形平緩和土地平整,這些都屬于有機(jī)質(zhì)聚集的有利的條件,并且水田一般為滯水狀態(tài),有機(jī)質(zhì)的分解速率較低,容易積累[23]。另一方面,相比旱地而言,水田作物秸稈殘存量多,大部分被微生物轉(zhuǎn)化為腐殖質(zhì),構(gòu)成了土壤有機(jī)質(zhì)的重要來(lái)源。

    從土壤類型看,該區(qū)石灰土、紫色土、水稻土和成土母巖不同,物質(zhì)來(lái)源不同,加之土壤自身的理化性質(zhì)不同,導(dǎo)致了不同土類中有機(jī)質(zhì)含量的差異。石灰土的成土母巖為三疊系碳酸鹽巖(灰?guī)r和白云巖),這類巖石中生物化石豐富,攜帶了較多的有機(jī)質(zhì)成分[30];同時(shí)石灰土本身的物質(zhì)成分也有利于有機(jī)質(zhì)的貯存[31]。紫色土、水稻土和黃壤由侏羅系或三疊系的紫色砂巖、泥巖及粉砂巖風(fēng)化而成,該類巖石形成于陸相環(huán)境,有機(jī)質(zhì)成分較少。紫色土呈中性或微堿性,易風(fēng)化、墾殖率高、水土流失嚴(yán)重,不利于有機(jī)質(zhì)的保存[32]。

    4 結(jié) 論

    (1)重慶西北部丘陵區(qū)土壤有機(jī)質(zhì)平均含量為16.26 g/kg,屬中等偏下水平。有機(jī)質(zhì)的空間變異系數(shù)為33.8 %,屬中等變異程度??臻g結(jié)構(gòu)采用指數(shù)型模型擬合較好,有機(jī)質(zhì)的塊金系數(shù)C0/C0+C為46.2 %,表明土壤有機(jī)質(zhì)具有中等空間相關(guān)格局,其空間變異性受結(jié)構(gòu)性因素和隨機(jī)因素影響程度相當(dāng)。

    (2)研究區(qū)土壤有機(jī)質(zhì)分布受土地利用類型、地形坡度、土壤類型等因素影響。不同土地利用類型下,林地土壤有機(jī)質(zhì)含量最高,水田次之。不同坡度下,低緩坡度中有機(jī)質(zhì)含量較高;海拔高度對(duì)有機(jī)質(zhì)含量影響不明顯。不同土壤類型,石灰土有機(jī)質(zhì)含量最高,黃壤和水稻土次之,紫色土含量最低。土地利用類型是有機(jī)質(zhì)空間分布最重要的影響因素。

    (3)基于土地利用類型修正的殘差法空間插值預(yù)測(cè)顯示,研究區(qū)土壤有機(jī)質(zhì)含量空間分布呈中間高,南北兩端低的趨勢(shì)。合川區(qū)水土流失較嚴(yán)重,土壤有機(jī)質(zhì)含量整體偏低。

    [1]黃昌勇. 土壤學(xué)[M]. 北京: 中國(guó)農(nóng)業(yè)出版社, 2000: 32.

    [2]韓 丹, 程先富, 謝金紅, 等. 大別山區(qū)江子河流域土壤有機(jī)質(zhì)的空間變異及其影響因素[J]. 土壤學(xué)報(bào), 2012, 49(2): 403-407.

    [3]苑小勇, 黃元仿, 高如泰, 等. 北京市平谷區(qū)農(nóng)用地土壤有機(jī)質(zhì)空間變異特征[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2008, 24(2): 70-75.

    [4]WARRICK A W. Analytical solutions to the one-dimensional linarized moisture flow equation for arbitrary input[J]. Soil Science, 1975, 120: 79-84.

    [5]CANN M D, HUMMEL J W, BROUER B H. Spatial analysis of soil fertility for site specific crop management[J]. Soil Sci Soc Am J, 1994, 58: 1240-1258

    [6]CAMBARDELLA C A, MOORMAN T B, NOVAK J M, et al. Field-scale variability of soil properties in central Iowa soils[J]. Soil Sci Soc Am, 1994, 58: 1501-1511.

    [7]MISHRA T K, BANERJEE S K. Spatial variability of soil pH and organic matter under Shorea robusta in lateritic region[J]. Indian Journal Forest, 1995, 18(2): 144-152.

    [8]趙明松, 張甘霖, 王德彩, 等. 徐淮黃泛平原土壤有機(jī)質(zhì)空間變異特征及主控因素分析[J]. 土壤學(xué)報(bào), 2013, 50(1):1-9.

    [9]楊 東, 劉 強(qiáng). 河西地區(qū)土壤全氮及有機(jī)質(zhì)空間變異特征分析-以張掖市甘州區(qū)為例[J]. 干旱地區(qū)農(nóng)業(yè)研究, 2010, 28(4): 183-186.

    [10]楊 葳. 紫色丘陵區(qū)土壤養(yǎng)分空間變異及肥力評(píng)價(jià)研究——以重慶市銅梁縣為例[D]. 重慶: 西南大學(xué), 2012.

    [11]王建鋒, 謝世友. 喀斯特石漠化地區(qū)土壤養(yǎng)分空問異質(zhì)性研究構(gòu)想[J]. 湖北農(nóng)業(yè)科學(xué), 2010, 49(7): 1605-1608.

    [12]邱 揚(yáng), 傅伯杰, 王 軍, 等. 黃土高原流域土壤養(yǎng)分的時(shí)空變異及其影響因子[J].自然科學(xué)進(jìn)展, 2004, 14: 294-299.

    [13]宋 莎, 李廷軒, 王永東, 等. 縣域農(nóng)田土壤有機(jī)質(zhì)空間變異及其影響因素分析[J]. 土壤, 2011, 43(1): 44-49.

    [14]中國(guó)地質(zhì)調(diào)查局. DD2005-01《多目標(biāo)區(qū)域地球化學(xué)調(diào)查規(guī)范》(1∶250000)[S]. 2005: 10.

    [15]中華人民共和國(guó)國(guó)土資源部. DZ/T 0295-2016,土地質(zhì)量地球化學(xué)評(píng)價(jià)規(guī)范[S]. 2016: 41.

    [16]陳延良, 雷國(guó)平. 松嫩平原南部土壤養(yǎng)分空間變異規(guī)律分析及其等級(jí)評(píng)價(jià)[J]. 揚(yáng)州大學(xué)學(xué)報(bào)(農(nóng)業(yè)與生命科學(xué)版), 2010, 31(2): 43-47.

    [17]蔣勇軍, 袁道先, 謝世友, 等. 典型巖溶流域土壤有機(jī)質(zhì)空間變異——以云南小江流域?yàn)槔齕J]. 生態(tài)學(xué)報(bào), 2007, 27(5):2041-2045.

    [18]劉愛利, 王培法, 丁 園. 地統(tǒng)計(jì)學(xué)概論[M]. 北京: 科學(xué)出版社, 2012: 60-67.

    [19]孔 毅, 劉國(guó)順, 臘貴曉, 等. 小尺度下基于GIS和地統(tǒng)計(jì)學(xué)的煙田土壤養(yǎng)分空間變異[J]. 江蘇農(nóng)業(yè)科學(xué), 2012, 40(4): 317-320.

    [20]SCHUIP C J E, VELDKAMP A. Long-term landscape-land use interactions as explaining factor for soil orgnic matter variability in Dutch agricultural landspace[J]. Geodema, 2008, 146: 457-465.

    [21]許 泉, 芮雯奕, 何 航, 等. 不同土地利用方式下中國(guó)農(nóng)田土壤有機(jī)碳密度特征及區(qū)域差異[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2006, 39(12): 2505-2510.

    [22]孫逸敏. 利用 SPSS軟件分析變量間的相關(guān)性[J]. 新疆教育學(xué)院學(xué)報(bào), 2007, 23(2):120-122.

    [23]李 婷, 張世熔, 劉 潯. 沱江流域中游土壤有機(jī)質(zhì)的空間變異特點(diǎn)及其影響因素[J]. 土壤學(xué)報(bào), 2011,48(4): 864-867.

    [24]奚小環(huán), 楊忠芳, 廖啟林, 等. 中國(guó)典型地區(qū)土壤碳儲(chǔ)量研究[J]. 第四紀(jì)研究, 2010, 30(3):573-583.

    [25]申廣榮, 葛曉燁, 黃秀梅. 上海崇明島表層土壤有機(jī)碳密度的空間分布特征及碳儲(chǔ)量估算[J]. 上海交通大學(xué)學(xué)報(bào)(農(nóng)業(yè)科學(xué)版), 2011, 29(6): 61-66.

    [26]劉 靜, 蔡國(guó)學(xué), 劉洪斌. 西南丘陵地區(qū)土壤有機(jī)質(zhì)含量的空間插值法研究[J]. 西南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2008, 30(3): 107-111.

    [27]文 雯, 周寶同, 汪亞峰, 等. 基于輔助環(huán)境變量的土壤有機(jī)碳空間插值——以黃土丘陵區(qū)小流域?yàn)槔齕J]. 生態(tài)學(xué)報(bào), 2013, 33(19): 6390-6395.

    [28]顧成軍, 史學(xué)正, 于東升. 結(jié)合土地利用與克里格插值的區(qū)域土壤有機(jī)碳空間表征方法[J]. 水土保持研究, 2014, 21(2):40-42.

    [29]范勝龍, 黃 炎, 林金石. 表征土壤有機(jī)碳區(qū)域分布的優(yōu)化空間插值模型研究[J]. 水土保持研究, 2011, 18(6): 1-5.

    [30]杜遠(yuǎn)生,童金南. 古生物地史學(xué)概論[M]. 北京: 中國(guó)地質(zhì)大學(xué)出版社, 1998: 150-194.

    [31]曹建華, 袁道先, 潘根興. 巖溶生態(tài)系統(tǒng)中的土壤[J]. 地球科學(xué)進(jìn)展, 2003, 18(1): 37-44.

    [32]倪九派, 袁道先, 謝德體, 等. 重慶巖溶區(qū)土壤有機(jī)碳庫(kù)的估算及其空間分布特征[J]. 生態(tài)學(xué)報(bào), 2009, 29(11):6292-6300.

    (責(zé)任編輯 李 潔)

    SpatialVariabilityofSoilOrganicMatterandItsInfluenceFactorsofHillyAreainNorthwestChongqing

    BAO Li-ran , ZHOU Jiao, LI Yu , JIA Zhong-min

    (Southeast Sichuan Geological Group, Chongqing Bureau of Geology and Minerals Exploration, Chongqing 400038, China)

    【Objective】The spatial variability of soil organic matter and its influence factors of hilly area in Northwest Chongqing ws discussed for providing scientific basis on agricultural production and land uses. 【Method】The geostatistics and GIS methods were used to study, spatial variability of soil organic matter and spatial interpolation based on land uses amendment.【Result】The soil organic matter average content was 16.26 g/kg, in the lower middle level and the variation coefficient was 33.8 %, in the middle level. The spatial distribution of soil organic matter was higher in middle, lower on both ends, and the whole Hechuan district was lower. 【Conclusion】Geostatistical analysis suggested that soil organic matter had medium spatial correlation; The spatial variability were resulted from structure factors(topography slope and soil types) and random factors(land uses, farming measures, ect.), and land uses was the dominated factor.

    Geostatistics;Soil organic matter; Spatial variability;Influence factors

    P596

    A

    1001-4829(2017)11-2541-07

    10.16213/j.cnki.scjas.2017.11.026

    2016-08-26

    全國(guó)土壤現(xiàn)狀調(diào)查及污染防治 ([2008]GZTR01-11);重慶市國(guó)土房管科技計(jì)劃項(xiàng)目(K J-2015043)

    鮑麗然(1983-),女,河北石家莊人,工程師,碩士,從事環(huán)境地球化學(xué)研究,Email:305946963@qq.com。

    猜你喜歡
    插值坡度土地利用
    關(guān)于公路超高漸變段合成坡度解析與應(yīng)用
    基于Sinc插值與相關(guān)譜的縱橫波速度比掃描方法
    土地利用生態(tài)系統(tǒng)服務(wù)研究進(jìn)展及啟示
    基于圖像處理的定位器坡度計(jì)算
    電氣化鐵道(2016年4期)2016-04-16 05:59:46
    一種改進(jìn)FFT多譜線插值諧波分析方法
    基于四項(xiàng)最低旁瓣Nuttall窗的插值FFT諧波分析
    濱??h土地利用挖潛方向在哪里
    坡度在巖石風(fēng)化層解譯中的應(yīng)用
    河北遙感(2015年2期)2015-07-18 11:11:14
    CT和MR對(duì)人上脛腓關(guān)節(jié)面坡度的比較研究
    Blackman-Harris窗的插值FFT諧波分析與應(yīng)用
    国产伦精品一区二区三区四那| 成年女人看的毛片在线观看| 乱码一卡2卡4卡精品| 免费在线观看成人毛片| 麻豆成人午夜福利视频| 国产欧美日韩一区二区精品| 亚洲精品久久国产高清桃花| 亚洲aⅴ乱码一区二区在线播放| 我要看日韩黄色一级片| 成人国产麻豆网| 婷婷精品国产亚洲av在线| 国产免费一级a男人的天堂| 精品熟女少妇av免费看| 成人亚洲欧美一区二区av| 少妇熟女aⅴ在线视频| 亚洲性久久影院| 日韩欧美在线乱码| 欧美高清成人免费视频www| 日日摸夜夜添夜夜添小说| 国产精品电影一区二区三区| 色av中文字幕| a级毛色黄片| 晚上一个人看的免费电影| 一卡2卡三卡四卡精品乱码亚洲| 99久国产av精品| 国产69精品久久久久777片| 男人狂女人下面高潮的视频| 91精品国产九色| 成熟少妇高潮喷水视频| 一卡2卡三卡四卡精品乱码亚洲| 欧美高清性xxxxhd video| 中文字幕免费在线视频6| 内地一区二区视频在线| 亚洲久久久久久中文字幕| 亚洲精品乱码久久久v下载方式| 久久精品国产亚洲av天美| 久久精品国产鲁丝片午夜精品| 国内久久婷婷六月综合欲色啪| 99久久成人亚洲精品观看| 一卡2卡三卡四卡精品乱码亚洲| 欧美另类亚洲清纯唯美| 亚洲第一电影网av| 亚洲久久久久久中文字幕| 色av中文字幕| 国内久久婷婷六月综合欲色啪| 最近最新中文字幕大全电影3| 久久久a久久爽久久v久久| 精品久久久噜噜| 美女被艹到高潮喷水动态| 亚洲三级黄色毛片| 日韩一区二区视频免费看| 久久久欧美国产精品| 搡老熟女国产l中国老女人| 国内精品一区二区在线观看| 日韩,欧美,国产一区二区三区 | 99在线人妻在线中文字幕| 亚洲性夜色夜夜综合| 日本免费a在线| 嫩草影院新地址| 欧美成人一区二区免费高清观看| 99riav亚洲国产免费| 99riav亚洲国产免费| 午夜免费激情av| 成人性生交大片免费视频hd| 99热6这里只有精品| 麻豆精品久久久久久蜜桃| 99久国产av精品| 精品人妻一区二区三区麻豆 | 日本黄色视频三级网站网址| 人妻夜夜爽99麻豆av| 国产精品久久久久久亚洲av鲁大| 91av网一区二区| 丰满的人妻完整版| 99在线人妻在线中文字幕| 国产精品精品国产色婷婷| 91在线精品国自产拍蜜月| 亚洲婷婷狠狠爱综合网| 麻豆国产av国片精品| 少妇熟女欧美另类| 欧美成人精品欧美一级黄| 久久精品国产亚洲av天美| 亚洲在线观看片| 国产av麻豆久久久久久久| 99精品在免费线老司机午夜| 亚洲在线观看片| 久久久久免费精品人妻一区二区| 黄色一级大片看看| av福利片在线观看| 国内少妇人妻偷人精品xxx网站| 亚洲激情五月婷婷啪啪| 91av网一区二区| 免费av毛片视频| 国产一区二区在线av高清观看| 一夜夜www| 欧美成人a在线观看| 国产亚洲av嫩草精品影院| 午夜福利视频1000在线观看| 亚洲一区高清亚洲精品| av黄色大香蕉| 亚洲av中文字字幕乱码综合| 欧美激情久久久久久爽电影| 12—13女人毛片做爰片一| 搡老岳熟女国产| 少妇的逼好多水| 69人妻影院| 好男人在线观看高清免费视频| 国产精品一二三区在线看| 亚洲av电影不卡..在线观看| 亚洲人与动物交配视频| 成人国产麻豆网| 1000部很黄的大片| 美女大奶头视频| 色综合亚洲欧美另类图片| 亚洲国产欧美人成| 国产成人a∨麻豆精品| 国产精品美女特级片免费视频播放器| 女的被弄到高潮叫床怎么办| 亚洲第一电影网av| 国产精品野战在线观看| 黄色配什么色好看| 不卡一级毛片| 美女cb高潮喷水在线观看| 久久韩国三级中文字幕| 一区二区三区免费毛片| 国产真实乱freesex| 免费看av在线观看网站| 精品人妻一区二区三区麻豆 | 变态另类成人亚洲欧美熟女| 国产高清视频在线播放一区| eeuss影院久久| 亚洲精品粉嫩美女一区| 丝袜美腿在线中文| 99热6这里只有精品| 日韩欧美国产在线观看| 国产精品不卡视频一区二区| 在线观看66精品国产| 最好的美女福利视频网| 91麻豆精品激情在线观看国产| 精品久久国产蜜桃| 国产成人a∨麻豆精品| 最近2019中文字幕mv第一页| 成年av动漫网址| 99热精品在线国产| 亚洲色图av天堂| 日韩一本色道免费dvd| 美女免费视频网站| 国产一区二区在线av高清观看| av卡一久久| 一级黄片播放器| 国产aⅴ精品一区二区三区波| 免费av观看视频| 婷婷亚洲欧美| 久久精品人妻少妇| 亚洲av.av天堂| 激情 狠狠 欧美| 俺也久久电影网| 晚上一个人看的免费电影| 国产精品野战在线观看| 1000部很黄的大片| 亚洲欧美日韩无卡精品| 日韩三级伦理在线观看| 午夜亚洲福利在线播放| 男人的好看免费观看在线视频| 又爽又黄无遮挡网站| 啦啦啦韩国在线观看视频| 搡老熟女国产l中国老女人| 国产精品一区二区性色av| 99国产极品粉嫩在线观看| 国产精品一及| 少妇裸体淫交视频免费看高清| 99久国产av精品国产电影| 成年av动漫网址| 欧美人与善性xxx| 1000部很黄的大片| 草草在线视频免费看| 丰满乱子伦码专区| 高清毛片免费观看视频网站| 啦啦啦啦在线视频资源| 亚洲av成人精品一区久久| 亚洲国产精品合色在线| 小说图片视频综合网站| 亚洲自拍偷在线| 亚洲国产欧洲综合997久久,| 夜夜爽天天搞| 可以在线观看的亚洲视频| 欧洲精品卡2卡3卡4卡5卡区| 精品久久久久久久人妻蜜臀av| 国产视频一区二区在线看| 久久人人精品亚洲av| 日韩欧美三级三区| 在线观看美女被高潮喷水网站| 亚洲欧美中文字幕日韩二区| 成人二区视频| 欧美高清成人免费视频www| 国内精品美女久久久久久| 婷婷六月久久综合丁香| 欧美性猛交╳xxx乱大交人| 男插女下体视频免费在线播放| 干丝袜人妻中文字幕| 99热这里只有是精品50| 国产伦在线观看视频一区| 国产精品综合久久久久久久免费| 婷婷六月久久综合丁香| 69人妻影院| 久久99热6这里只有精品| 美女被艹到高潮喷水动态| 老熟妇乱子伦视频在线观看| 午夜福利18| 哪里可以看免费的av片| av卡一久久| 晚上一个人看的免费电影| 成人欧美大片| 寂寞人妻少妇视频99o| 日韩精品有码人妻一区| av在线天堂中文字幕| 国产美女午夜福利| 乱码一卡2卡4卡精品| 色综合色国产| 午夜福利高清视频| 国产av不卡久久| 18禁裸乳无遮挡免费网站照片| 99久国产av精品| 小蜜桃在线观看免费完整版高清| 日本熟妇午夜| 在线观看av片永久免费下载| 免费av毛片视频| 丝袜美腿在线中文| 韩国av在线不卡| 一级毛片aaaaaa免费看小| 国产精品一区二区三区四区免费观看 | 亚洲无线观看免费| 夜夜看夜夜爽夜夜摸| 国产高清激情床上av| 亚洲七黄色美女视频| 亚洲中文字幕一区二区三区有码在线看| 欧美一区二区精品小视频在线| 91麻豆精品激情在线观看国产| 免费搜索国产男女视频| 如何舔出高潮| 欧美高清成人免费视频www| 日韩欧美 国产精品| 国产高清激情床上av| a级毛片免费高清观看在线播放| 麻豆成人午夜福利视频| 亚洲美女搞黄在线观看 | 成人漫画全彩无遮挡| 简卡轻食公司| 中文字幕久久专区| 精品熟女少妇av免费看| 菩萨蛮人人尽说江南好唐韦庄 | 一级黄片播放器| 可以在线观看毛片的网站| 精品福利观看| 欧美+亚洲+日韩+国产| 亚洲av.av天堂| 欧美xxxx黑人xx丫x性爽| 久久久国产成人精品二区| 亚洲成av人片在线播放无| 国产爱豆传媒在线观看| 亚洲熟妇中文字幕五十中出| 免费观看的影片在线观看| 一本精品99久久精品77| 插逼视频在线观看| 国产蜜桃级精品一区二区三区| 亚洲国产精品成人综合色| 男人狂女人下面高潮的视频| 亚洲,欧美,日韩| 少妇人妻一区二区三区视频| 偷拍熟女少妇极品色| 干丝袜人妻中文字幕| 国产亚洲精品av在线| 国产爱豆传媒在线观看| 亚洲熟妇中文字幕五十中出| 免费观看的影片在线观看| 乱码一卡2卡4卡精品| 国产熟女欧美一区二区| 男女做爰动态图高潮gif福利片| 一级毛片电影观看 | 国产成人一区二区在线| 欧美又色又爽又黄视频| 噜噜噜噜噜久久久久久91| 韩国av在线不卡| 国产精品人妻久久久久久| 午夜福利在线在线| 一级毛片久久久久久久久女| 国产精品一二三区在线看| 亚洲精华国产精华液的使用体验 | 99国产极品粉嫩在线观看| 成人漫画全彩无遮挡| 桃色一区二区三区在线观看| 欧美激情国产日韩精品一区| 亚洲国产精品sss在线观看| 亚洲av熟女| 别揉我奶头 嗯啊视频| 九九在线视频观看精品| 日韩欧美三级三区| 成人性生交大片免费视频hd| 人人妻人人看人人澡| 2021天堂中文幕一二区在线观| 少妇裸体淫交视频免费看高清| 久久99热6这里只有精品| 伦精品一区二区三区| 人人妻人人澡欧美一区二区| 我要看日韩黄色一级片| 久久精品国产99精品国产亚洲性色| 日韩,欧美,国产一区二区三区 | 波多野结衣巨乳人妻| 亚洲美女黄片视频| 成人午夜高清在线视频| 最近最新中文字幕大全电影3| 12—13女人毛片做爰片一| 我要看日韩黄色一级片| 亚洲欧美精品综合久久99| 欧美不卡视频在线免费观看| 菩萨蛮人人尽说江南好唐韦庄 | 国产亚洲精品综合一区在线观看| 亚洲精品一区av在线观看| 久久综合国产亚洲精品| 女的被弄到高潮叫床怎么办| 变态另类成人亚洲欧美熟女| 亚洲18禁久久av| 少妇人妻精品综合一区二区 | 欧美高清性xxxxhd video| 国产熟女欧美一区二区| 久久久久久国产a免费观看| 国产单亲对白刺激| 人人妻人人澡人人爽人人夜夜 | 午夜福利18| av在线观看视频网站免费| 人妻久久中文字幕网| 麻豆久久精品国产亚洲av| 欧美成人精品欧美一级黄| 91精品国产九色| 精品无人区乱码1区二区| 久久久久久久亚洲中文字幕| 亚洲国产精品合色在线| 久久久久久久久中文| 一个人看的www免费观看视频| 99久久中文字幕三级久久日本| 欧美色视频一区免费| 免费av观看视频| 久久国产乱子免费精品| videossex国产| 亚洲三级黄色毛片| 亚洲乱码一区二区免费版| 亚洲av一区综合| 久久精品夜色国产| 亚洲国产精品合色在线| 91麻豆精品激情在线观看国产| 精品一区二区三区视频在线| 欧美激情国产日韩精品一区| 搡女人真爽免费视频火全软件 | 欧美成人精品欧美一级黄| 亚洲精品一卡2卡三卡4卡5卡| 免费不卡的大黄色大毛片视频在线观看 | 久久精品国产清高在天天线| 麻豆国产av国片精品| 国产午夜精品论理片| 美女免费视频网站| 成人无遮挡网站| 十八禁国产超污无遮挡网站| 色哟哟哟哟哟哟| 国产精品精品国产色婷婷| 成熟少妇高潮喷水视频| 九九在线视频观看精品| 国产一区亚洲一区在线观看| 国产精品一二三区在线看| 三级经典国产精品| 免费高清视频大片| 日本撒尿小便嘘嘘汇集6| 在线a可以看的网站| 亚洲天堂国产精品一区在线| 国产精品久久久久久久电影| 国国产精品蜜臀av免费| 成年女人毛片免费观看观看9| 亚洲成人av在线免费| 99在线人妻在线中文字幕| 麻豆av噜噜一区二区三区| 欧美绝顶高潮抽搐喷水| 久久99热这里只有精品18| 国产人妻一区二区三区在| 日本熟妇午夜| 亚洲av免费在线观看| 啦啦啦啦在线视频资源| 欧美三级亚洲精品| 国产精品女同一区二区软件| 国产av一区在线观看免费| 午夜福利在线观看免费完整高清在 | 免费无遮挡裸体视频| 变态另类丝袜制服| 国产精品一二三区在线看| 免费无遮挡裸体视频| 国产精品日韩av在线免费观看| 亚洲美女搞黄在线观看 | 成人无遮挡网站| 免费av观看视频| 国产精品,欧美在线| 麻豆av噜噜一区二区三区| 日本-黄色视频高清免费观看| 欧美性猛交╳xxx乱大交人| 中文在线观看免费www的网站| 自拍偷自拍亚洲精品老妇| 18禁裸乳无遮挡免费网站照片| 天堂av国产一区二区熟女人妻| 国内揄拍国产精品人妻在线| av女优亚洲男人天堂| 久久精品国产亚洲av涩爱 | 亚洲av成人av| 久久精品国产99精品国产亚洲性色| 亚洲第一区二区三区不卡| 久久精品影院6| 如何舔出高潮| 嫩草影视91久久| 中国美女看黄片| 人人妻人人看人人澡| 精品一区二区免费观看| 少妇人妻精品综合一区二区 | 在线播放国产精品三级| 亚洲激情五月婷婷啪啪| 18禁在线无遮挡免费观看视频 | 亚州av有码| 国产v大片淫在线免费观看| 国产女主播在线喷水免费视频网站 | 国产不卡一卡二| 51国产日韩欧美| av福利片在线观看| 俄罗斯特黄特色一大片| 国产在线精品亚洲第一网站| 蜜臀久久99精品久久宅男| 久久天躁狠狠躁夜夜2o2o| 在线国产一区二区在线| 国产精品国产三级国产av玫瑰| 国产精品伦人一区二区| 国内精品美女久久久久久| 69av精品久久久久久| 成人欧美大片| 蜜桃久久精品国产亚洲av| 18+在线观看网站| 久久99热这里只有精品18| 国产日本99.免费观看| 亚洲电影在线观看av| 欧美激情在线99| 日韩一区二区视频免费看| 成人国产麻豆网| 免费搜索国产男女视频| 亚洲久久久久久中文字幕| 久久久午夜欧美精品| 亚洲av免费高清在线观看| 亚洲成人av在线免费| 我的女老师完整版在线观看| 国产色爽女视频免费观看| 国产免费男女视频| 亚洲自偷自拍三级| 亚洲av不卡在线观看| 午夜福利高清视频| 一区二区三区免费毛片| 伦精品一区二区三区| 变态另类丝袜制服| 成熟少妇高潮喷水视频| 久久欧美精品欧美久久欧美| 一本久久中文字幕| 少妇熟女aⅴ在线视频| 露出奶头的视频| 午夜精品一区二区三区免费看| 欧美日韩精品成人综合77777| 精品一区二区三区视频在线| 成年版毛片免费区| 亚洲最大成人av| 国产高清三级在线| 少妇丰满av| 夜夜夜夜夜久久久久| 麻豆久久精品国产亚洲av| 99riav亚洲国产免费| 特级一级黄色大片| 国国产精品蜜臀av免费| 久久久a久久爽久久v久久| 搡女人真爽免费视频火全软件 | 免费av观看视频| 久久久久久国产a免费观看| 国产色爽女视频免费观看| 婷婷六月久久综合丁香| 晚上一个人看的免费电影| 人人妻人人澡欧美一区二区| 久久久国产成人免费| 夜夜看夜夜爽夜夜摸| 免费搜索国产男女视频| 熟妇人妻久久中文字幕3abv| 中文字幕熟女人妻在线| 国产探花在线观看一区二区| 国产精品久久久久久精品电影| 成人二区视频| 男女啪啪激烈高潮av片| 好男人在线观看高清免费视频| 波野结衣二区三区在线| 欧美激情在线99| 欧美+日韩+精品| 亚洲最大成人av| 人人妻人人看人人澡| 精品久久久久久久久av| 久久人妻av系列| 香蕉av资源在线| 色尼玛亚洲综合影院| 热99在线观看视频| 老司机午夜福利在线观看视频| 最好的美女福利视频网| 在线免费十八禁| 国产午夜福利久久久久久| 在线观看午夜福利视频| 亚洲av五月六月丁香网| 成人午夜高清在线视频| 伦精品一区二区三区| 美女xxoo啪啪120秒动态图| 尾随美女入室| 俄罗斯特黄特色一大片| 在线播放无遮挡| 91久久精品国产一区二区三区| 最近最新中文字幕大全电影3| 亚洲欧美成人综合另类久久久 | 国产精品免费一区二区三区在线| 高清毛片免费看| 久久精品夜夜夜夜夜久久蜜豆| 午夜精品国产一区二区电影 | 亚洲精品乱码久久久v下载方式| 18禁黄网站禁片免费观看直播| 成人鲁丝片一二三区免费| 人人妻人人澡欧美一区二区| 国产日本99.免费观看| 国产男靠女视频免费网站| 97超级碰碰碰精品色视频在线观看| 亚洲激情五月婷婷啪啪| 一进一出抽搐动态| 中文字幕av成人在线电影| 最近2019中文字幕mv第一页| 亚洲人成网站在线播放欧美日韩| 男人的好看免费观看在线视频| 女同久久另类99精品国产91| 中国美女看黄片| 三级男女做爰猛烈吃奶摸视频| 久久天躁狠狠躁夜夜2o2o| .国产精品久久| 99riav亚洲国产免费| 亚洲国产高清在线一区二区三| 99久久精品国产国产毛片| 不卡一级毛片| 麻豆成人午夜福利视频| ponron亚洲| 非洲黑人性xxxx精品又粗又长| 神马国产精品三级电影在线观看| a级毛色黄片| 一级a爱片免费观看的视频| 精品99又大又爽又粗少妇毛片| 男人舔奶头视频| 美女黄网站色视频| 又黄又爽又刺激的免费视频.| 成年女人永久免费观看视频| 国产精品一区二区免费欧美| 亚洲熟妇熟女久久| 麻豆一二三区av精品| 色吧在线观看| 久久九九热精品免费| 欧美+日韩+精品| 美女内射精品一级片tv| 我的老师免费观看完整版| 桃色一区二区三区在线观看| 男人舔女人下体高潮全视频| 精品一区二区三区av网在线观看| or卡值多少钱| 精品福利观看| 国产真实乱freesex| 国产精品福利在线免费观看| 国产av一区在线观看免费| 十八禁网站免费在线| a级一级毛片免费在线观看| 美女xxoo啪啪120秒动态图| 99久久九九国产精品国产免费| 禁无遮挡网站| 久久精品人妻少妇| 特大巨黑吊av在线直播| 久久久久久伊人网av| 亚洲三级黄色毛片| 成人特级av手机在线观看| 看免费成人av毛片| 国内精品宾馆在线| 亚洲美女黄片视频| 国产综合懂色| av免费在线看不卡| 免费av观看视频| 六月丁香七月| 国产免费一级a男人的天堂| 狂野欧美白嫩少妇大欣赏| 亚洲精品日韩av片在线观看| 国产高清三级在线| 国产国拍精品亚洲av在线观看| 亚洲无线观看免费| 久久久久九九精品影院| 亚洲av中文av极速乱| 亚洲欧美精品综合久久99| 少妇被粗大猛烈的视频| av在线老鸭窝| 欧美一区二区国产精品久久精品| 欧美最黄视频在线播放免费| 小说图片视频综合网站| 亚洲最大成人手机在线| 日韩欧美 国产精品| 秋霞在线观看毛片| 国产三级在线视频| 一卡2卡三卡四卡精品乱码亚洲| 国产一区亚洲一区在线观看| 真人做人爱边吃奶动态| 日本精品一区二区三区蜜桃| 亚洲天堂国产精品一区在线| 亚洲av免费在线观看| 中文字幕av在线有码专区| 日日干狠狠操夜夜爽| 亚洲国产高清在线一区二区三|