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

    興凱湖平原表層土壤有機碳空間變異的主控因素

    2022-10-28 05:15:50楊澤張一鶴戴慧敏劉國棟劉凱許江
    物探與化探 2022年5期
    關(guān)鍵詞:成土母質(zhì)土壤有機

    楊澤,張一鶴,戴慧敏,劉國棟,劉凱,許江

    (1.中國地質(zhì)調(diào)查局 沈陽地質(zhì)調(diào)查中心,遼寧 沈陽 110034;2.自然資源部 黑土地演化與生態(tài)效應(yīng)重點實驗室,遼寧 沈陽 110034;3.遼寧省黑土地演化與生態(tài)效應(yīng)重點實驗室,遼寧 沈陽 110034)

    0 引言

    土壤作為陸地最大的碳庫[1]和最大的有機碳庫[2],在全球碳循環(huán)中起著重要的作用。土壤有機碳(SOC)含量及其動態(tài)變化直接影響著全球碳循環(huán),也影響著土壤質(zhì)量和植物生長[3-5]。由于受結(jié)構(gòu)性與隨機性因素的影響,土壤有機碳在空間分布上存在一定的異質(zhì)性,即使在相鄰位置,土壤有機碳含量也會因為復(fù)雜的環(huán)境而產(chǎn)生較大差異[6-9]。因此,準確地獲取區(qū)域尺度內(nèi)土壤有機碳含量及空間變異的主控因素,對土壤有機碳的調(diào)控和恢復(fù)、農(nóng)業(yè)可持續(xù)發(fā)展以及全球環(huán)境保護都具有重要意義。

    目前,國內(nèi)外許多學者[6-13]針對不同尺度上的土壤有機碳空間變異的主控因素開展了研究,結(jié)果表明:在眾多的影響因素中,定性因素如成土母質(zhì)、植被類型或土地利用方式以及土壤類型等對土壤有機碳空間分布的影響顯著[14-17]。不同尺度上土壤有機碳空間分布的定性主控因素也存在明顯的區(qū)域差異,如:河北省省域尺度上土地利用方式對土壤有機碳空間變異的解釋能力大于土壤類型[11];徐淮黃泛平原上土地利用方式對土壤有機碳空間分布的影響大于成土母質(zhì)和土壤類型[10];江蘇省域尺度上成土母質(zhì)的影響大于土地利用方式[18];江西省域尺度上秸稈還田的影響大于成土母質(zhì)和土地利用方式[6];縣級尺度上,福建龍海市[19]的土壤類型較土地利用方式對土壤有機碳的空間分布影響更大,川中丘陵縣域[7]土地利用方式較土壤類型對土壤有機碳含量的空間分布影響更大。以上研究表明,成土母質(zhì)、土壤類型以及土地利用方式等定性因素對土壤有機碳空間分布的影響隨著研究區(qū)域和尺度的不同而發(fā)生變化,存在明顯的區(qū)域特征和尺度效應(yīng)[7]。此外,以往研究多集中于結(jié)構(gòu)性因素對土壤有機碳空間變異的研究,對隨機性因素考慮較少。因此,充分地考慮研究區(qū)內(nèi)土壤的結(jié)構(gòu)性(地形地貌、成土母質(zhì)、土壤類型等)與隨機性(土地利用方式、土壤開墾年限等)因素對準確把握土壤有機碳空間變異的主控因素、掌握調(diào)控土壤有機碳的關(guān)鍵因子等方面具有非常重要的現(xiàn)實意義。

    興凱湖平原位于黑龍江省東部偏南,是我國東北黑土區(qū)重要的商品糧基地之一,但多年來一直未開展過系統(tǒng)的多目標地球化學調(diào)查,有關(guān)整個平原區(qū)土壤有機碳的研究也鮮有報道。本研究基于2019年中國地質(zhì)調(diào)查局沈陽地質(zhì)調(diào)查中心在興凱湖平原開展的1∶25萬土地質(zhì)量地球化學調(diào)查表層土壤數(shù)據(jù),運用地統(tǒng)計分析、半變異函數(shù)分析和回歸分析等方法,結(jié)合成土母質(zhì)、土壤質(zhì)地、土壤類型、土地利用方式及開墾年限等影響因素,對興凱湖平原土壤有機碳空間變異的主控因素進行定量研究,以期為興凱湖平原土壤有機碳調(diào)控、恢復(fù)及農(nóng)業(yè)可持續(xù)發(fā)展提供科學依據(jù)。

    1 研究區(qū)概況

    興凱湖平原位于黑龍江省東部偏南,東與俄羅斯接壤,西至張廣才嶺,北依完達山,南抵興凱湖,總面積約16 000 km2,行政區(qū)跨越虎林市、密山市等。研究區(qū)屬中緯度寒溫帶濕潤、半濕潤大陸性季風氣候,冬季漫長,嚴寒而干燥,夏季短促,溫暖且濕潤,年平均溫度2.9~3.1 ℃。區(qū)內(nèi)雨量充沛,多年平均降水量為526~710 mm,每年6~9月降水占全年降水量總量的70%;蒸發(fā)量為降水量的1.5~2倍,每年5~9月蒸發(fā)量占全年總量的70%。

    研究區(qū)地形總體由西向東傾斜,坡降1/6 000~1/10 000,平原區(qū)地勢低平,微有起伏,僅在虎林市和虎頭鎮(zhèn)散布著被地塹切割殘余的孤山。地貌上,西部為侵蝕低山、丘陵區(qū),東部為沖積—湖積平原區(qū)。區(qū)內(nèi)水系比較發(fā)育,主要河流由南往北依次為穆棱河、七虎林河及阿布沁河,均為烏蘇里江左翼支流。

    區(qū)內(nèi)成土母質(zhì)以第四紀沖積、湖積物及不同時期的砂巖風化物為主。土壤類型以白漿土為主,其次為暗棕壤、沼澤土及草甸土。主要土地利用為耕地,以水田為主,其次為林地和沼澤地。糧食作物以水稻為主,兼有大豆和玉米(圖1)。

    圖1 研究區(qū)高程、采樣點(a)及土地利用類型(b)Fig.1 Elevation, sampling sites(a) and land use(b) of the study area

    2 數(shù)據(jù)來源及研究方法

    2.1 樣品采集與測試

    依據(jù)《多目標區(qū)域地球化學調(diào)查規(guī)范(1∶250 000)》(DZ/T 0258—2014)[20]和《土地質(zhì)量地球化學評價規(guī)范》(DZ/T 0295—2016)[21],中國地質(zhì)調(diào)查局沈陽地質(zhì)調(diào)查中心于2019年在興凱湖平原開展1∶25萬土地質(zhì)量地球化學調(diào)查工作,網(wǎng)格化采集土壤表層和深層樣品,表層土壤采樣密度為1個點/km2,采樣深度為0~20 cm(土柱),1個點/4 km2組合分析;深層土壤采樣密度為1個點/4 km2,采樣深度為150~200 cm(土柱),1個點/16 km2組合分析。樣點均采用GPS進行定位,采集后的樣品進行晾干、碎樣、過20目篩,采取四分法,稱取200 g后裝紙袋送至實驗室進行測試。共獲得表層土壤組合樣品4 151件,深層土壤組合樣品1 000余件。本次數(shù)據(jù)源主要采用表層土壤樣品數(shù)據(jù)。

    土壤樣品測試由遼寧省地質(zhì)礦產(chǎn)研究院有限責任公司承擔,土壤有機碳(SOC)采用重鉻酸鉀法。樣品分析過程中采用國家一級標準土樣監(jiān)控分析測試的準確度,采用重復(fù)樣監(jiān)控分析測試的精密度。一級標準物質(zhì)的所有分析指標合格率為100%,重復(fù)樣合格率符合《多目標區(qū)域地球化學調(diào)查規(guī)范(1∶250 000)》(DZ/T 0258—2014)中的樣品分析質(zhì)量控制要求,測試指標報出率均達到100%,異常點重復(fù)性檢驗合格率95.8%。

    2.2 其他數(shù)據(jù)來源

    本文所用其他數(shù)據(jù)來源說明如下:

    土地利用數(shù)據(jù):來源于中國科學院地理科學與資源研究所,中國土地利用現(xiàn)狀遙感監(jiān)測數(shù)據(jù)庫是在國家科技支撐計劃、中國科學院知識創(chuàng)新工程重要方向項目等多項重大科技項目的支持下,經(jīng)過多年的積累而建立的覆蓋全國陸地區(qū)域的多時相土地利用現(xiàn)狀數(shù)據(jù)庫。采用的4期數(shù)據(jù)集包括20世紀70年代末期、1995年、2010年和2018年,數(shù)據(jù)生產(chǎn)制作是以各期Landsat TM/ETM遙感影像為主要數(shù)據(jù)源,通過人工目視解譯生成,數(shù)據(jù)分辨率為1 km柵格數(shù)據(jù)。

    土壤類數(shù)據(jù):來源于中國科學院南京土壤研究所,原數(shù)據(jù)來自全國土壤普查辦公室組織完成的《1∶100萬中華人民共和國土壤圖》,由中國科學院南京土壤研究所完成數(shù)據(jù)庫建設(shè),屬性字段包含了土綱、土類及亞類。

    2.3 數(shù)據(jù)處理方法

    經(jīng)典統(tǒng)計分析、相關(guān)分析以及回歸分析均在SPSS 25軟件中完成;考慮到成土母質(zhì)、土壤類型、土地利用方式和開墾年限為定性分類變量,研究中采用啞變量(虛擬變量)進行賦值,相關(guān)方法參見文獻[6-7]。土壤有機碳數(shù)據(jù)進行對數(shù)變換后,在GS+9.0中進行半變異函數(shù)分析。由于原始數(shù)據(jù)不符合正態(tài)分布,故在ArcGIS 10.2中采用反距離權(quán)重法獲得土壤有機碳空間分布圖。

    2.3.1 化學蝕變指數(shù)

    化學蝕變指數(shù)(chemical index of alteration,CIA)是判別由硅酸鹽礦物組成的沉積物風化程度最常用的化學指標,該指標也經(jīng)常被用作土壤風化程度的指標[22]。其計算公式為

    (1)

    式中氧化物均為分子摩爾數(shù);CaO*為硅酸鹽相,采用McLennan計算方法[23]。CIA反映了長石風化成黏土礦物的程度,CIA值越高,指示氣候溫暖濕潤,風化程度越高;反之,寒冷干燥,風化程度低。CIA值介于50~65為低化學風化程度,CIA值介于65~85為中化學風化程度,CIA值介于85~100為強化學風化[24-27]。

    2.3.2 土地開墾年限

    使用ArcGIS10.2值提取至點功能,將4期遙感解譯土地利用狀況提取至采樣點,在Excel中根據(jù)采樣點不同時期的土地利用情況計算采樣點位置處的土地開墾年限,一直為草地、林地、濕地等自然狀態(tài)的土地,開墾年限為0,其他土地利用類型的開墾年限以實際情況計算為準,截止計算年為2019年。

    3 結(jié)果與討論

    3.1 興凱湖平原表層土壤有機碳統(tǒng)計特征

    研究區(qū)表層土壤有機碳含量變化范圍(表1)為0.35%~14.49%,平均值2.80%,變異系數(shù)為0.44,屬于中等變異水平,表明研究區(qū)表層土壤有機碳空間分布存在一定的不均。從偏度來看,原始數(shù)據(jù)不符合正態(tài)分布,呈現(xiàn)右偏;從峰度來看,數(shù)據(jù)比較陡峭。對表層土壤有機碳含量數(shù)據(jù)進行對數(shù)變換后,數(shù)據(jù)近似服從正態(tài)分布(圖2)。與以往研究結(jié)果對比(表2),研究區(qū)表層土壤有機碳高于東北黑土、東北平原以及中國大陸的平均值,低于黑龍江省的平均值(1986~1990年)。對照《土地質(zhì)量地球化學評價規(guī)范》(DZ/T 0295—2016)中土壤養(yǎng)分指標等級劃分標準的有機質(zhì)標準,興凱湖平原土壤有機質(zhì)(折算后)含量達到豐富、較豐富水平。

    表1 表層土壤有機碳含量統(tǒng)計Table 1 Statistics of SOC content in topsoil

    圖2 研究區(qū)土壤有機碳含量及其對數(shù)轉(zhuǎn)換值的頻率分布Fig.2 Histograms of SOC contents and its log-transformed values in the study area

    表2 研究區(qū)與其他地區(qū)表層土壤有機碳含量對比Table 2 Comparison of SOC content in topsoil between Xingkai Lake Plain and other regions

    3.2 半變異函數(shù)分析

    半方差分析能較好地刻畫土壤有機碳空間分布的隨機性和結(jié)構(gòu)性。運用GS+軟件對研究區(qū)土壤有機碳含量進行半變異函數(shù)擬合(表3、圖3)。根據(jù)決定系數(shù)(R2)、殘差(RSS)等,選擇變異函數(shù)的擬合模型。由表3可知,研究區(qū)土壤有機碳變異函數(shù)比較符合球形模型,R2為0.962,擬合效果較好。塊金值為0.016,說明在當前的采樣尺度范圍內(nèi)存在由采樣誤差、短距離的變異、隨機因素引起的變異。塊金效應(yīng)為47.06%,表明區(qū)域內(nèi)存在中等強度的空間自相關(guān)性,說明研究區(qū)土壤有機碳含量空間分布受結(jié)構(gòu)性因素和隨機性因素的共同影響且更傾向于受結(jié)構(gòu)性因素的影響。其中,結(jié)構(gòu)性因素主要是指土壤在形成過程中的母質(zhì)類型、土壤類型等,隨機性因素則是指諸如土地利用方式等能在一定程度上反映耕作活動、種植制度、投入管理水平差異的人為活動因素;結(jié)構(gòu)性因素使得土壤屬性具有空間自相關(guān)性,而隨機性因素則會減弱土壤屬性的空間自相關(guān)性,增大其異質(zhì)性。變程(A)為135.7 km,說明研究區(qū)土壤有機碳自相關(guān)范圍較大,可能與研究區(qū)平原地貌有關(guān)。

    表3 研究區(qū)土壤有機碳半方差函數(shù)模型及其參數(shù)Table 3 Theoretical semivariogram model for SOC content and parameters values in the study area

    圖3 研究區(qū)土壤有機碳半方差函數(shù)Fig.3 Isotropic semivariogram of SOC in the study area

    3.3 有機碳空間變異特征

    為了更直觀反映研究區(qū)土壤有機碳的空間變異特征,采用反距離加權(quán)空間插值方法獲得興凱湖平原表層土壤有機碳含量的空間分布。如圖4所示,研究區(qū)表層土壤有機碳呈現(xiàn)“中、西部低,東、北部高”的分布格局,體現(xiàn)了興凱湖平原土壤有機碳空間分布不均的特點。高值區(qū)(>3%)主要分布在研究區(qū)東部的興凱湖農(nóng)場—忠誠鄉(xiāng)、北部的珍寶島鄉(xiāng)—東方紅鎮(zhèn)—云山農(nóng)場及西北部的向陽村附近;低值區(qū)(<2%)主要分布在研究區(qū)中、西部地區(qū)。土壤有機碳含量以2%~3%為主,占研究區(qū)總面積的49.35%;其次為3%~5%,占研究區(qū)總面積的30.64%。有機碳總體水平較高,與統(tǒng)計結(jié)果一致。大部分高值區(qū)均在海拔較低處。

    圖4 研究區(qū)土壤有機碳空間分布Fig.4 Spatial distribution of SOC in the studyarea

    3.4 土壤有機碳空間變異的影響因素

    3.4.1 成土母質(zhì)

    成土母質(zhì)對土壤有機碳含量的影響主要表現(xiàn)在對土壤礦物質(zhì)組成、土壤風化程度與土壤質(zhì)地等方面的影響。對研究區(qū)不同母質(zhì)發(fā)育的土壤有機碳含量進行統(tǒng)計(表4),結(jié)果表明,不同成土母質(zhì)間土壤有機碳均值差異明顯,其中三疊紀變質(zhì)巖風化物發(fā)育的土壤其有機碳含量最高,平均值為4.13%;其次為三疊紀板巖及粉砂巖風化物(3.20%)、元古宙變質(zhì)巖風化物(3.15%)、中更新世松散堆積物(3.08%)及太古宙變質(zhì)巖風化物(3.03%);新近紀砂巖、砂礫巖風化物發(fā)育的土壤有機碳含量最低,其平均值僅為1.81%。同種巖性不同時代的風化物,土壤有機碳含量也有差異。以玄武巖為例,全新世玄武巖風化物發(fā)育而成的土壤有機碳含量明顯低于新近系玄武巖風化物發(fā)育的土壤。從變異情況來看,各成土母質(zhì)下土壤有機碳含量變異系數(shù)為0.18~0.65,除新近紀砂巖、砂礫巖母質(zhì)為低變異外,其余成土母質(zhì)變異系數(shù)均為中等變異。

    表4 研究區(qū)不同成土母質(zhì)土壤有機碳含量統(tǒng)計特征Table 4 Descriptive statistics characteristics of SOC content with different soil parent materials in the study area

    3.4.2 土壤質(zhì)地

    土壤質(zhì)地是影響土壤有機碳含量的重要因素。已有研究表明,土壤機械組成影響著土壤的耕性、通透性、保水保肥性能[32-33],特別是砂粒和黏粒含量,可以顯著影響土壤有機質(zhì)的累積和礦物組成[34]。本次工作未測試土壤的質(zhì)地,但可采用獲得的地球化學指標側(cè)面反映土壤的質(zhì)地,如硅鋁鐵率[35]、土壤化學蝕變指數(shù)(CIA)等。根據(jù)二者與土壤有機碳的相關(guān)系數(shù)大小,本次采用相關(guān)系數(shù)較大的土壤化學蝕變指數(shù)(CIA)。CIA值與長石風化成黏土礦物的程度成正比,CIA值越大,風化強度越大,相應(yīng)的土壤的粒級會越細。研究區(qū)土壤風化指數(shù)平均值為66.21,范圍為52.62~80.82,變異系數(shù)為0.05,為低變異水平。由表5可知,隨著土壤有機碳含量的增加,CIA平均值呈增大的趨勢,表明土壤有機碳含量與CIA具有良好的相關(guān)關(guān)系。據(jù)相關(guān)性分析,土壤有機碳含量與CIA呈現(xiàn)非常顯著的正相關(guān)關(guān)系(P<0.01),相關(guān)系數(shù)為0.549。以上分析表明,土壤質(zhì)地與有機碳的空間分布具有很好的關(guān)聯(lián)性。

    表5 研究區(qū)不同土壤有機碳含量土壤風化指數(shù)統(tǒng)計特征Table 5 Descriptive statistics characteristics of CIA with different SOC content in the study area

    3.4.3 土壤類型

    研究區(qū)不同土壤類型有機碳含量統(tǒng)計結(jié)果表明(表6),不同土壤類型有機碳含量具有一定差別。沼澤土平均值最高,為3.28%,泥炭土、草甸土、暗棕壤、白漿土及黑土次之,平均值在2.54%~2.88%,水稻土最低,平均值為2.47%。各土類有機碳含量之間差距不大,最大值僅為最小值的1.4倍,反映出土壤類型對有機碳的空間分布影響不大。從變異系數(shù)來看,黑土的變異系數(shù)最小,為0.24,暗棕壤的變異系數(shù)最大,為0.50。7種土類的變異系數(shù)平均值為0.4,與全區(qū)土壤有機碳含量的平均變異系數(shù)接近。

    表6 研究區(qū)不同土壤類型有機碳含量統(tǒng)計特征Table 6 Descriptive statistics characteristics of SOC contents under different soil types in the study area

    3.4.4 土地利用方式

    不同土地利用方式下土壤有機碳含量平均值差異明顯(表7)。草地有機碳平均值最高,達3.33%,其次為濕地、水田、有林地,達2.92~3.08%,灘地最低,為2.16%。耕地中水田土壤有機碳高于旱田,約為其1.2倍。從變異系數(shù)來看,居民用地的變異系數(shù)最小,為0.28;其次為水田(0.33),灌木林最大,為0.54,其余土地利用方式變異系數(shù)在0.4左右。

    表7 研究區(qū)不同土地利用方式下土壤有機碳含量統(tǒng)計特征Table 7 Descriptive statistics characteristics of SOC contents under different land use types in the study area

    3.4.5 開墾年限

    不同開墾年限土壤有機碳含量變化明顯,隨著開墾年限的增加,土壤有機碳基本呈下降趨勢(表8)。開墾5年的土壤有機碳平均值最高,達4.37%,開墾10年的平均值為3.78%,開墾15~20年為2.7%左右,25年略有升高,開墾40年的土壤有機碳平均值最低,為2.45%,為開墾5年土壤有機碳的56%。未開墾的土壤多位于山區(qū),少部分位于濕地,受立地條件差異較大,故土壤有機碳平均值略低,變異系數(shù)較大。

    表8 研究區(qū)不同開墾年限土壤有機碳含量統(tǒng)計特征Table 8 Descriptive statistics characteristics of SOC contents under different reclamation years in the study area

    3.5 主控因素研究

    為定量解釋各影響因素對土壤有機碳的空間變異的獨立解釋能力,對成土母質(zhì)、土壤質(zhì)地、土壤類型、土地利用方式及開墾年限進行回歸分析(表9)。所有自變量均達到極顯著水平(P<0.01)。結(jié)果表明,成土母質(zhì)、土壤質(zhì)地、土壤類型、土地利用方式及開墾年限等因素對研究區(qū)土壤有機碳均有顯著性影響,其中土壤質(zhì)地對研究區(qū)土壤有機碳空間變異影響最大,能夠獨立解釋空間變異的30.1%,顯著高于其他因素;其次為土地利用方式、成土母質(zhì)、土壤類型和開墾年限。這5項影響因素合計能夠解釋41.2%的有機碳空間變異。這表明興凱湖平原土壤有機碳空間變異的主控因素為土壤質(zhì)地,其次為土地利用方式。

    表9 研究區(qū)不同因素對土壤有機碳的回歸分析結(jié)果Table 9 Regression analysis of SOC with different factors in the study area

    4 結(jié)論與討論

    4.1 結(jié)構(gòu)性因素對土壤有機碳的影響

    結(jié)構(gòu)性因素中土壤質(zhì)地的獨立解釋能力最高,遠高于成土母質(zhì)和土壤類型,是興凱湖平原土壤有機碳空間變異的主控因素。該結(jié)果與趙明松等[18]的研究結(jié)果一致,獨立解釋能力均在30%左右。這主要是因為土壤中黏粒含量較高會增加物理和化學保護機制,即黏粒通過與有機物質(zhì)結(jié)合形成有機—無機復(fù)合體,吸附穩(wěn)定土壤有機質(zhì)(SOM),降低了SOM的礦化速度,有利于有機質(zhì)的積累;相反,砂粒含量較高的土壤中則減少或缺少這種保護機制,有機質(zhì)礦化分解速率加快,導(dǎo)致有機質(zhì)含量低[36-37]。工作區(qū)受地形地貌及河流沉積作用的影響,土壤質(zhì)地空間分布不均(見圖5)。工作區(qū)西部地區(qū)屬低山丘陵區(qū),CIA值普遍小于65,個別點位CIA值略大,整體反映西部土壤風化程度較低;東部和北部地區(qū)CIA值普遍大于65,在北部的七虎林河以北地區(qū)普遍大于68,東部的興凱湖以東地區(qū)也存在大于68的集中分布地區(qū),表明東部平原區(qū)土壤整體風化程度較西部地區(qū)高。土壤質(zhì)地的分布于有機碳空間分布高度吻合,表明土壤質(zhì)地對土壤有機碳的空間變異影響較大。

    圖5 研究區(qū)表層土壤CIA分布Fig.5 CIA distribution of topsoil in the study area

    成土母質(zhì)對土壤的物理性狀和化學組成均有極其重要的影響。在一定的地質(zhì)構(gòu)造環(huán)境內(nèi),土壤元素含量(如有機碳)和土壤性狀與母質(zhì)有著緊密的關(guān)系[38],但在成土母質(zhì)的劃分上,無法劃分到粒級的級別,如第四紀松散堆積物僅劃分為全新世松散堆積物、晚更新世松散堆積物、中更新世松散堆積物,不能很好地對應(yīng)質(zhì)地的分級,因此成土母質(zhì)的獨立解釋能力要低。土壤類型亦是如此。有研究表明[6-7,39],土壤類型對土壤有機碳空間變異性的獨立解釋能力與土壤分類級別呈反比,即土壤分類級別越低解釋土壤有機碳空間變異性的能力越強。這是因為低級別的土壤類型分類是在高級別的分類基礎(chǔ)上疊加更加精細的指標進一步分類的結(jié)果。本次未搜集到土屬一級的分類,故未開展相關(guān)研究。

    4.2 隨機性因素對土壤有機碳的影響

    本次研究的隨機性因素中,土地利用方式是僅次于土壤質(zhì)地的影響土壤有機碳空間變異的因素,高于土壤開墾年限。土地利用方式作為自然條件和人為活動的綜合反映,通過有機物的輸入以及控制土壤有機碳積累和釋放速度來影響土壤有機碳含量[11,16]。已有研究表明,不同土地利用方式下土壤有機碳含量差異顯著[12,16],同一研究區(qū)內(nèi)土地利用方式對土壤有機碳空間變異的獨立解釋能力高于土類和亞類,但低于土屬[10-11]。本次研究結(jié)果與已有研究結(jié)果基本一致。研究區(qū)土壤有機碳平均值呈現(xiàn)草地>濕地>水田>有林地>灌木林>旱田>居民用地>灘地。耕地中水田土壤有機碳含量明顯高于旱田,是旱田的1.2倍。這是因為水田根茬殘留量大,相當于部分秸稈還田,有機質(zhì)歸還量比旱田高,大部分被微生物稍經(jīng)分解即轉(zhuǎn)化為腐殖質(zhì)儲存于土壤中,且水田長期處于淹水狀態(tài),有機質(zhì)分解緩慢,易積累,因而土壤有機碳含量較旱田高。土地利用方式對有機碳空間變異的獨立解釋能力為9.2%,大于土壤類型的3.8%。以上分析表明,土地利用方式也是研究區(qū)土壤有機碳空間變異的重要因素之一,可通過改變工作區(qū)土地利用方式進行土壤碳庫的調(diào)控,比如地勢低洼的土地退耕還濕、還草;也可通過秸稈還田等保護性耕作制度達到提高土壤碳庫的目標。

    土壤開墾年限對于土壤有機碳空間變異在4種定性因素中解釋能力最低,與土壤類型基本相當??赡茉蛴袃煞矫妫环矫媸潜敬嗡鸭耐恋乩闷诖屋^少,對于開墾年限的確定不夠精確,導(dǎo)致絕大部分采樣點位于0、25和40年,其他年限采樣點數(shù)量太少;另一方面,土壤開墾年限對有機碳的影響可能主要體現(xiàn)在土壤有機碳的變化,因此其在解釋土壤有機碳變化方面的效果應(yīng)該要好于土壤有機碳的空間變異。

    4.3 結(jié)論

    1)興凱湖平原表層土壤有機碳含量在0.35%~14.49%,平均值2.80%,這一含量水平高于東北黑土、東北平原以及中國大陸的平均值,低于黑龍江省的平均值(1986~1990年)。變異系數(shù)分別為0.44,屬于等程度的空間變異性。土壤有機質(zhì)(折算后)含量達到豐富、較豐富水平。

    2)土壤有機碳空間自相關(guān)范圍較大,其空間變異更趨向于受結(jié)構(gòu)性因素的影響,在全區(qū)表現(xiàn)為“中、西部低,東、北部高”的分布格局。

    3)在各影響因素中,成土母質(zhì)對土壤有機碳空間變異的獨立解釋能力為6.8%,土壤質(zhì)地的獨立解釋能力為30.1%,土壤類型能獨立解釋3.8%的有機碳空間變異,土地利用方式能獨立解釋9.2%,開墾年限能獨立解釋3.3%的有機碳空間變異。土壤質(zhì)地是影響興凱湖平原表層土壤有機碳空間變異的主控因子。

    致謝:文中使用的土壤有機碳來源于多目標區(qū)域地球化學調(diào)查,該項目由自然資源部中國地質(zhì)調(diào)查局沈陽地質(zhì)調(diào)查中心組織實施,由于參與人員眾多,無法一一列舉。筆者感謝該項目參加人員付出的勞動及基礎(chǔ)資料的支持。在此特別感謝江葉楓老師、羅由林老師給予的幫助,以及審稿老師提出的寶貴意見。

    猜你喜歡
    成土母質(zhì)土壤有機
    281例兒童眼部毛母質(zhì)瘤的臨床及病理分析
    試論礦山生態(tài)修復(fù)的地質(zhì)成土
    土壤侵蝕作用是雛形土和新成土廣泛分布的成因研究
    貴州成土母巖類型及其與耕地土壤關(guān)系探討
    西雙版納橡膠林土壤有機碳分布特征研究
    貴州母巖(母質(zhì))對土壤類型及分布的影響
    秸稈還田的土壤有機碳周轉(zhuǎn)特征
    土壤有機碳轉(zhuǎn)化研究及其進展
    AMDIS在土壤有機污染物鑒別中表征性統(tǒng)計量的探究
    兒童多發(fā)性毛母質(zhì)瘤1例
    美女 人体艺术 gogo| 制服人妻中文乱码| 成人鲁丝片一二三区免费| 国产v大片淫在线免费观看| 99热这里只有精品一区 | 曰老女人黄片| 一级黄色大片毛片| 午夜福利高清视频| 国产精品影院久久| 亚洲美女黄片视频| 亚洲片人在线观看| 美女黄网站色视频| 国产一区二区在线av高清观看| 舔av片在线| 成人午夜高清在线视频| 国产精品久久久久久人妻精品电影| 日本 欧美在线| 香蕉国产在线看| 最近在线观看免费完整版| 很黄的视频免费| 蜜桃久久精品国产亚洲av| xxx96com| 人人妻人人看人人澡| 色老头精品视频在线观看| 国产极品精品免费视频能看的| 久久久久久九九精品二区国产| 波多野结衣高清无吗| 精华霜和精华液先用哪个| 黄色女人牲交| 长腿黑丝高跟| 成人特级av手机在线观看| 国产伦人伦偷精品视频| www.www免费av| 国产成人一区二区三区免费视频网站| 成人三级做爰电影| 性色avwww在线观看| 91久久精品国产一区二区成人 | 香蕉av资源在线| 日韩三级视频一区二区三区| 亚洲av第一区精品v没综合| 又大又爽又粗| 久久久国产成人免费| 免费在线观看亚洲国产| 国产成人欧美在线观看| 欧美极品一区二区三区四区| 久久人妻av系列| 国产欧美日韩精品一区二区| 看黄色毛片网站| 亚洲国产精品合色在线| 国产免费男女视频| 99视频精品全部免费 在线 | bbb黄色大片| 成熟少妇高潮喷水视频| 观看免费一级毛片| 国产精品99久久久久久久久| 亚洲精品一卡2卡三卡4卡5卡| 成人av一区二区三区在线看| 亚洲真实伦在线观看| 午夜a级毛片| 老汉色av国产亚洲站长工具| 亚洲无线在线观看| 最近视频中文字幕2019在线8| 偷拍熟女少妇极品色| cao死你这个sao货| 麻豆成人午夜福利视频| 欧美成人免费av一区二区三区| 一本精品99久久精品77| 国产成人精品久久二区二区免费| 婷婷亚洲欧美| 全区人妻精品视频| 男女午夜视频在线观看| 丰满的人妻完整版| 国产成年人精品一区二区| 欧美一区二区国产精品久久精品| 亚洲av熟女| 在线观看免费午夜福利视频| 国产v大片淫在线免费观看| 国内精品久久久久精免费| 天天一区二区日本电影三级| 日本黄大片高清| 男人舔女人下体高潮全视频| 亚洲熟妇熟女久久| 国产成人精品无人区| 成人av在线播放网站| 欧美另类亚洲清纯唯美| 国产熟女xx| 国内精品美女久久久久久| 欧美av亚洲av综合av国产av| 老汉色av国产亚洲站长工具| 两个人视频免费观看高清| 男人的好看免费观看在线视频| e午夜精品久久久久久久| 最近最新中文字幕大全电影3| 国产亚洲精品av在线| 亚洲国产高清在线一区二区三| 亚洲自拍偷在线| 免费高清视频大片| 性色avwww在线观看| 熟女少妇亚洲综合色aaa.| 中文字幕高清在线视频| 女同久久另类99精品国产91| 精华霜和精华液先用哪个| 色综合站精品国产| 一本久久中文字幕| 国产单亲对白刺激| 一区二区三区国产精品乱码| 欧美日韩精品网址| 男女午夜视频在线观看| 亚洲国产欧美人成| www.www免费av| 亚洲国产中文字幕在线视频| 少妇的逼水好多| 三级男女做爰猛烈吃奶摸视频| 村上凉子中文字幕在线| 亚洲精品美女久久久久99蜜臀| 男女视频在线观看网站免费| 1024手机看黄色片| 老司机福利观看| 精品国内亚洲2022精品成人| 国产黄片美女视频| av中文乱码字幕在线| 哪里可以看免费的av片| 国模一区二区三区四区视频 | av国产免费在线观看| 舔av片在线| 国产美女午夜福利| 日本在线视频免费播放| 国产亚洲精品一区二区www| 1000部很黄的大片| 欧美黄色片欧美黄色片| 母亲3免费完整高清在线观看| 国产精品亚洲av一区麻豆| 国产1区2区3区精品| 日本黄色视频三级网站网址| 男人舔女人下体高潮全视频| 90打野战视频偷拍视频| 欧美黑人巨大hd| 熟女人妻精品中文字幕| 久久午夜综合久久蜜桃| 波多野结衣巨乳人妻| 午夜两性在线视频| 国产又色又爽无遮挡免费看| 极品教师在线免费播放| 麻豆国产av国片精品| 欧洲精品卡2卡3卡4卡5卡区| 真人做人爱边吃奶动态| 精品不卡国产一区二区三区| 一a级毛片在线观看| 男人的好看免费观看在线视频| 色尼玛亚洲综合影院| 国语自产精品视频在线第100页| 91av网一区二区| 俄罗斯特黄特色一大片| 国产一区二区三区在线臀色熟女| 欧美性猛交╳xxx乱大交人| 国模一区二区三区四区视频 | 亚洲在线自拍视频| 成人永久免费在线观看视频| av在线蜜桃| 久久久久久久久久黄片| 十八禁人妻一区二区| 99热这里只有精品一区 | 美女午夜性视频免费| 精品久久久久久久毛片微露脸| 欧美日韩一级在线毛片| 国产视频内射| 精品久久久久久久人妻蜜臀av| 免费在线观看影片大全网站| 岛国视频午夜一区免费看| 国产精品久久视频播放| 一卡2卡三卡四卡精品乱码亚洲| 国产伦一二天堂av在线观看| 国产乱人视频| 床上黄色一级片| 午夜精品在线福利| 久久人妻av系列| 亚洲成人免费电影在线观看| 热99在线观看视频| 99在线人妻在线中文字幕| 国模一区二区三区四区视频 | 中文字幕人妻丝袜一区二区| 日本a在线网址| 亚洲av成人av| 国产伦精品一区二区三区视频9 | 在线观看日韩欧美| 亚洲精品中文字幕一二三四区| 美女 人体艺术 gogo| 一级毛片高清免费大全| 人妻夜夜爽99麻豆av| 国产免费av片在线观看野外av| 男女之事视频高清在线观看| 亚洲成人久久爱视频| 亚洲欧美日韩高清在线视频| 久久精品91蜜桃| 国产精品一及| 成人无遮挡网站| 久久伊人香网站| 成年版毛片免费区| 国产高清激情床上av| а√天堂www在线а√下载| 性欧美人与动物交配| 在线观看舔阴道视频| 日本成人三级电影网站| 国产高潮美女av| 国产欧美日韩精品亚洲av| 观看美女的网站| 欧美zozozo另类| aaaaa片日本免费| 淫妇啪啪啪对白视频| 国产亚洲av高清不卡| 欧美午夜高清在线| 我要搜黄色片| 欧美中文日本在线观看视频| 免费在线观看日本一区| 欧美一级a爱片免费观看看| 九九在线视频观看精品| 久久亚洲精品不卡| 黄色视频,在线免费观看| 久久性视频一级片| 亚洲在线自拍视频| 国产精品亚洲一级av第二区| 亚洲欧美日韩无卡精品| 国产伦人伦偷精品视频| 日本a在线网址| 国产成人啪精品午夜网站| 亚洲色图av天堂| 熟妇人妻久久中文字幕3abv| 国产三级在线视频| 九色成人免费人妻av| 黄片小视频在线播放| 久久精品91无色码中文字幕| 偷拍熟女少妇极品色| 精品无人区乱码1区二区| 国产精品电影一区二区三区| 日本免费一区二区三区高清不卡| 国产黄a三级三级三级人| 给我免费播放毛片高清在线观看| 国产成人一区二区三区免费视频网站| 免费看十八禁软件| 狠狠狠狠99中文字幕| 特大巨黑吊av在线直播| 成人国产综合亚洲| 久久这里只有精品中国| 亚洲精品美女久久久久99蜜臀| 这个男人来自地球电影免费观看| 婷婷丁香在线五月| 亚洲精品一卡2卡三卡4卡5卡| 国内精品美女久久久久久| 日本黄大片高清| 波多野结衣高清作品| АⅤ资源中文在线天堂| 欧美精品啪啪一区二区三区| 日韩免费av在线播放| 亚洲性夜色夜夜综合| 欧美黄色淫秽网站| 在线十欧美十亚洲十日本专区| 午夜福利欧美成人| 久久天躁狠狠躁夜夜2o2o| 国产又黄又爽又无遮挡在线| 成年女人看的毛片在线观看| 亚洲av成人av| 精品一区二区三区四区五区乱码| 中文字幕最新亚洲高清| 欧美激情久久久久久爽电影| 亚洲精品中文字幕一二三四区| 男插女下体视频免费在线播放| 日本一二三区视频观看| 看免费av毛片| 国产 一区 欧美 日韩| 亚洲一区二区三区色噜噜| 一二三四社区在线视频社区8| 禁无遮挡网站| 中文字幕最新亚洲高清| 国产亚洲精品av在线| 亚洲片人在线观看| 色播亚洲综合网| 成年女人毛片免费观看观看9| 日韩大尺度精品在线看网址| 日日摸夜夜添夜夜添小说| 在线视频色国产色| 色吧在线观看| 欧美性猛交╳xxx乱大交人| 精品99又大又爽又粗少妇毛片 | 日韩欧美在线二视频| 国产精品亚洲美女久久久| 18禁国产床啪视频网站| 最近视频中文字幕2019在线8| 成人亚洲精品av一区二区| 欧美日韩瑟瑟在线播放| 免费电影在线观看免费观看| 男人的好看免费观看在线视频| 亚洲在线自拍视频| 亚洲精品一卡2卡三卡4卡5卡| 国产精品一及| 国产亚洲精品久久久com| 午夜福利成人在线免费观看| 不卡av一区二区三区| 国产亚洲精品久久久久久毛片| 国产成年人精品一区二区| 亚洲欧美一区二区三区黑人| 精品午夜福利视频在线观看一区| 麻豆一二三区av精品| 在线观看66精品国产| 亚洲,欧美精品.| 中文字幕熟女人妻在线| 色av中文字幕| 99国产精品一区二区三区| 999久久久精品免费观看国产| 亚洲av电影在线进入| 日本在线视频免费播放| 国产亚洲av高清不卡| 国产高清视频在线观看网站| 亚洲成a人片在线一区二区| 在线国产一区二区在线| 国内揄拍国产精品人妻在线| 色播亚洲综合网| 欧美zozozo另类| 国产探花在线观看一区二区| 日本免费一区二区三区高清不卡| 男女之事视频高清在线观看| 青草久久国产| 亚洲av中文字字幕乱码综合| 人妻夜夜爽99麻豆av| 欧美日韩乱码在线| 国产又黄又爽又无遮挡在线| 波多野结衣巨乳人妻| 国产不卡一卡二| 精品欧美国产一区二区三| 狂野欧美白嫩少妇大欣赏| 757午夜福利合集在线观看| 欧美精品啪啪一区二区三区| 美女午夜性视频免费| 色播亚洲综合网| 99久久精品一区二区三区| 嫩草影视91久久| 亚洲黑人精品在线| 丰满人妻熟妇乱又伦精品不卡| 99热这里只有是精品50| 久久久国产欧美日韩av| 亚洲黑人精品在线| 国产欧美日韩精品一区二区| 婷婷精品国产亚洲av在线| 欧美xxxx黑人xx丫x性爽| 最新美女视频免费是黄的| 亚洲美女黄片视频| 51午夜福利影视在线观看| 一卡2卡三卡四卡精品乱码亚洲| 一二三四在线观看免费中文在| 亚洲国产欧美人成| 麻豆国产av国片精品| 久久久久久久久中文| 美女午夜性视频免费| 听说在线观看完整版免费高清| 国产av一区在线观看免费| 两个人的视频大全免费| 十八禁网站免费在线| 成年女人毛片免费观看观看9| 国产极品精品免费视频能看的| 床上黄色一级片| 亚洲国产中文字幕在线视频| ponron亚洲| 九九久久精品国产亚洲av麻豆 | 99久久精品一区二区三区| 老司机在亚洲福利影院| 久久久久久大精品| 99久久成人亚洲精品观看| 又粗又爽又猛毛片免费看| 国内毛片毛片毛片毛片毛片| 黄色丝袜av网址大全| 两个人看的免费小视频| 欧美一区二区精品小视频在线| 一个人观看的视频www高清免费观看 | 午夜福利免费观看在线| 国产人伦9x9x在线观看| 无遮挡黄片免费观看| 欧美性猛交╳xxx乱大交人| 久久国产精品人妻蜜桃| 午夜福利在线在线| 毛片女人毛片| 日日夜夜操网爽| 亚洲片人在线观看| 十八禁人妻一区二区| 岛国在线观看网站| 男女做爰动态图高潮gif福利片| 欧美zozozo另类| 99热这里只有是精品50| 国内精品美女久久久久久| 在线永久观看黄色视频| e午夜精品久久久久久久| 久久久精品大字幕| 一级a爱片免费观看的视频| 又大又爽又粗| 精品日产1卡2卡| 久久久久亚洲av毛片大全| 国产野战对白在线观看| 人人妻,人人澡人人爽秒播| 亚洲五月婷婷丁香| 99国产精品99久久久久| 九九久久精品国产亚洲av麻豆 | 无限看片的www在线观看| 午夜福利欧美成人| 三级毛片av免费| 999久久久精品免费观看国产| 日本五十路高清| 性色avwww在线观看| 亚洲av中文字字幕乱码综合| 精品乱码久久久久久99久播| 国内精品久久久久久久电影| 亚洲av片天天在线观看| 精品国产美女av久久久久小说| 国产麻豆成人av免费视频| 国产在线精品亚洲第一网站| 精华霜和精华液先用哪个| 亚洲aⅴ乱码一区二区在线播放| 欧美乱码精品一区二区三区| 最新中文字幕久久久久 | tocl精华| 日韩精品中文字幕看吧| 欧美在线黄色| 91在线精品国自产拍蜜月 | 欧美一级a爱片免费观看看| 午夜福利在线在线| 久9热在线精品视频| 亚洲av成人精品一区久久| 色综合亚洲欧美另类图片| 婷婷丁香在线五月| 一个人看的www免费观看视频| 国产伦人伦偷精品视频| 在线观看舔阴道视频| 国产一区在线观看成人免费| 日韩有码中文字幕| 亚洲av五月六月丁香网| 日韩欧美免费精品| АⅤ资源中文在线天堂| 99热这里只有精品一区 | 男人舔女人下体高潮全视频| aaaaa片日本免费| 女同久久另类99精品国产91| 免费电影在线观看免费观看| 欧美日韩亚洲国产一区二区在线观看| 久久人人精品亚洲av| 91久久精品国产一区二区成人 | 真人一进一出gif抽搐免费| 免费观看精品视频网站| 亚洲在线观看片| 婷婷六月久久综合丁香| e午夜精品久久久久久久| 午夜福利高清视频| 亚洲av第一区精品v没综合| 美女被艹到高潮喷水动态| 黄片小视频在线播放| 国产精品 国内视频| 欧美激情久久久久久爽电影| 欧美日韩乱码在线| 99精品久久久久人妻精品| 国产成人精品久久二区二区91| 国内久久婷婷六月综合欲色啪| 日本五十路高清| 欧美乱妇无乱码| 亚洲中文字幕日韩| 亚洲av美国av| 99久久成人亚洲精品观看| 亚洲精品在线美女| 国产精品久久久人人做人人爽| 国产视频内射| 国产精品久久久久久亚洲av鲁大| 日本一二三区视频观看| 国产又黄又爽又无遮挡在线| 久久99热这里只有精品18| 法律面前人人平等表现在哪些方面| 午夜日韩欧美国产| 90打野战视频偷拍视频| 午夜激情欧美在线| 大型黄色视频在线免费观看| 成人性生交大片免费视频hd| 在线十欧美十亚洲十日本专区| 观看美女的网站| 欧美高清成人免费视频www| 精品国产超薄肉色丝袜足j| 18禁裸乳无遮挡免费网站照片| 欧美绝顶高潮抽搐喷水| 男女床上黄色一级片免费看| 黄色丝袜av网址大全| 搡老岳熟女国产| 在线观看日韩欧美| 夜夜看夜夜爽夜夜摸| 欧美在线一区亚洲| 国产成+人综合+亚洲专区| 香蕉丝袜av| 亚洲成人久久爱视频| 老熟妇仑乱视频hdxx| 欧美成狂野欧美在线观看| 白带黄色成豆腐渣| 制服丝袜大香蕉在线| 亚洲中文日韩欧美视频| 国产精品日韩av在线免费观看| 久久热在线av| 2021天堂中文幕一二区在线观| 99久久国产精品久久久| 麻豆成人av在线观看| 哪里可以看免费的av片| 真实男女啪啪啪动态图| 国产私拍福利视频在线观看| 法律面前人人平等表现在哪些方面| 熟女电影av网| 免费看十八禁软件| 国产精品98久久久久久宅男小说| cao死你这个sao货| 国内精品久久久久精免费| 九色成人免费人妻av| 99久久精品一区二区三区| 视频区欧美日本亚洲| 亚洲国产欧美网| 亚洲av美国av| 两个人看的免费小视频| 中文字幕人妻丝袜一区二区| 国产成人啪精品午夜网站| 国产精品久久久人人做人人爽| 国产高清videossex| 嫩草影视91久久| 亚洲成a人片在线一区二区| 无遮挡黄片免费观看| 夜夜夜夜夜久久久久| 亚洲国产高清在线一区二区三| 母亲3免费完整高清在线观看| 国产精品亚洲一级av第二区| a级毛片a级免费在线| 欧美乱码精品一区二区三区| av欧美777| 色尼玛亚洲综合影院| 午夜福利欧美成人| 久久精品国产99精品国产亚洲性色| 精品国产超薄肉色丝袜足j| 岛国在线免费视频观看| 少妇熟女aⅴ在线视频| 岛国视频午夜一区免费看| 高潮久久久久久久久久久不卡| 老司机午夜十八禁免费视频| 日韩中文字幕欧美一区二区| 国产高清有码在线观看视频| 亚洲欧美精品综合一区二区三区| 听说在线观看完整版免费高清| 亚洲av成人一区二区三| 男人舔奶头视频| 性欧美人与动物交配| 少妇丰满av| 亚洲avbb在线观看| 久久久久亚洲av毛片大全| 久久精品国产99精品国产亚洲性色| 美女午夜性视频免费| 欧美黄色淫秽网站| 亚洲片人在线观看| 国产激情久久老熟女| 亚洲中文字幕一区二区三区有码在线看 | 午夜激情欧美在线| 日韩免费av在线播放| 夜夜夜夜夜久久久久| 黄色成人免费大全| 亚洲欧美日韩高清专用| 日韩中文字幕欧美一区二区| 在线观看免费午夜福利视频| 最新美女视频免费是黄的| 国产一区二区在线av高清观看| 色综合婷婷激情| 在线免费观看的www视频| 真人做人爱边吃奶动态| 制服丝袜大香蕉在线| 国产精品99久久99久久久不卡| 天堂动漫精品| 啦啦啦观看免费观看视频高清| 国产亚洲欧美在线一区二区| 精品国产乱子伦一区二区三区| 国产三级黄色录像| 麻豆一二三区av精品| 91av网站免费观看| 免费av不卡在线播放| 国产成年人精品一区二区| 亚洲中文日韩欧美视频| 91在线精品国自产拍蜜月 | 国产亚洲av高清不卡| 欧美乱妇无乱码| 久久欧美精品欧美久久欧美| 亚洲第一欧美日韩一区二区三区| 久久婷婷人人爽人人干人人爱| 欧美大码av| 琪琪午夜伦伦电影理论片6080| 国产精品乱码一区二三区的特点| 午夜福利在线观看吧| 三级国产精品欧美在线观看 | 亚洲午夜精品一区,二区,三区| 男女视频在线观看网站免费| 久久国产精品影院| 成人性生交大片免费视频hd| 日韩中文字幕欧美一区二区| 精品一区二区三区四区五区乱码| 中文字幕人成人乱码亚洲影| 国产高清videossex| 国产一级毛片七仙女欲春2| 日本撒尿小便嘘嘘汇集6| 国产不卡一卡二| 禁无遮挡网站| 午夜免费成人在线视频| 午夜成年电影在线免费观看| 国产91精品成人一区二区三区| 色综合欧美亚洲国产小说| 美女黄网站色视频| 亚洲午夜理论影院| 每晚都被弄得嗷嗷叫到高潮| 久久久国产成人精品二区| 久久婷婷人人爽人人干人人爱| 级片在线观看| 亚洲av免费在线观看| 看黄色毛片网站| cao死你这个sao货| 免费在线观看影片大全网站|