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

    巖溶山區(qū)土壤水的地球化學(xué)特征及其與土壤環(huán)境的關(guān)系

    2017-12-05 06:05:44覃星銘何丙輝吳華英
    關(guān)鍵詞:排序特征分析

    覃星銘 何丙輝 吳華英

    (1.西南大學(xué)資源環(huán)境學(xué)院, 重慶 400716; 2.中國地質(zhì)科學(xué)院巖溶地質(zhì)研究所, 桂林 541004)

    巖溶山區(qū)土壤水的地球化學(xué)特征及其與土壤環(huán)境的關(guān)系

    覃星銘1,2何丙輝1吳華英2

    (1.西南大學(xué)資源環(huán)境學(xué)院, 重慶 400716; 2.中國地質(zhì)科學(xué)院巖溶地質(zhì)研究所, 桂林 541004)

    巖溶峰叢山區(qū); 土壤水; 地球化學(xué)性質(zhì); 土壤環(huán)境因子; 灰度關(guān)聯(lián); 典范對應(yīng)分析

    引言

    作為全球巖溶三大分布中心之一的中國西南地區(qū),境內(nèi)碳酸鹽巖類巖層縱深橫廣,巖溶現(xiàn)象普遍、多樣,類型繁雜、發(fā)育強(qiáng)烈,在相對脆弱的巖溶多重介質(zhì)與不合理的人為活動的疊加影響下,引發(fā)植被破壞、土地質(zhì)量衰退、水土流失、干旱缺水、群眾生活貧困等一系列的問題[1-2]。

    圖1 研究區(qū)采樣點(diǎn)的環(huán)境地質(zhì)簡圖Fig.1 Schematic environmental geology map of monitoring point in research area

    土壤水是大氣降雨經(jīng)巖溶表層帶中的土壤向地下水轉(zhuǎn)化的重要紐帶,也是巖溶關(guān)鍵帶植物、微生物生存發(fā)展的物質(zhì)和水分來源,直接關(guān)系到巖溶關(guān)鍵帶水資源的形成、轉(zhuǎn)化與消耗。土壤作為土壤水的空間存儲載體和物質(zhì)能量來源,可認(rèn)為土壤質(zhì)量條件是響應(yīng)土壤水的重要環(huán)境因子。目前,有部分專家學(xué)者對巖溶峰叢山區(qū)土壤異質(zhì)性及其環(huán)境的影響進(jìn)行了研究[3-9],也取得了很多研究成果[10-11],但對巖溶峰叢山區(qū)土壤水的地球化學(xué)性質(zhì)方面的研究較少[12-13],特別是利用灰度關(guān)聯(lián)分析(Grey relational analysis,GRA)和典范對應(yīng)分析(Canonical correspondence analysis,CCA)等方法討論巖溶山區(qū)土壤水地球化學(xué)性質(zhì)及其環(huán)境響應(yīng)方面的研究鮮有報道[14-15]。因此,本文選擇巖溶峰叢山區(qū)土壤-土壤水為研究對象,研究土壤水地球化學(xué)性質(zhì)的空間變化特征及其對土壤環(huán)境條件的響應(yīng)關(guān)系,以期完善土壤水與環(huán)境關(guān)系的方法研究、指導(dǎo)巖溶峰叢山區(qū)生態(tài)重建和農(nóng)林水土工程調(diào)蓄。

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)位于桂西南右江中游石漠化最為嚴(yán)重的中心地帶之一的廣西平果縣果化鎮(zhèn)布堯村一帶,面積約3 km2,由多組高低錯落的錐狀山峰聯(lián)座組合形成巖溶峰叢洼地地貌,海拔150~500 m,為典型的巖溶峰叢山區(qū)。區(qū)內(nèi)出露的地層主要有二疊系下統(tǒng)棲霞組(P1q)、茅口組(P1m)和石炭系上統(tǒng)(C3)[16],成土母巖和主要含水層為中厚層狀灰?guī)r,地表水系不發(fā)育。受人為活動和巖溶地質(zhì)結(jié)構(gòu)的影響,研究區(qū)內(nèi)石漠化現(xiàn)象嚴(yán)重,土層淺薄,坡地土壤厚度為0~60 cm,土被不連續(xù)分布,土壤主要為棕褐色石灰土。植被類型單一,主要以稀疏林地、低矮灌草叢為主,耕種地主要以火龍果(HylocereusundulatusBritt)為主,種植面積較大,間種黃豆(Glycinemax)和花生(Arachishypogaea)等農(nóng)作物。研究區(qū)所設(shè)置的7個樣地為旱耕地土被條件,均為二疊系下統(tǒng)棲霞組地層,地表為散流狀表層巖溶帶(圖1)。

    1.2 樣品采集與分析

    于2015年4—8月在研究區(qū)內(nèi)對7個樣地進(jìn)行野外觀測和采樣分析,挖開土壤剖面按0~20 cm、20~40 cm分層采集土壤樣品,樣品的采集按剖面深度分別選取4~6個樣點(diǎn),分層采集樣點(diǎn)土壤混合風(fēng)干,經(jīng)研磨過篩后備用。土壤樣品測定由廣西農(nóng)業(yè)科學(xué)院土壤分析檢測中心依據(jù)《土壤分析技術(shù)規(guī)范》及常規(guī)的方法測定[17-18]。土壤水由自行研制的滲流土壤水收集器分層進(jìn)行水樣采集[19],樣品由國土資源部巖溶地質(zhì)資源環(huán)境監(jiān)督檢測中心依據(jù)地下水標(biāo)準(zhǔn)檢驗方法進(jìn)行水化學(xué)性質(zhì)的測定。

    1.3 數(shù)據(jù)處理方法與統(tǒng)計分析

    運(yùn)用Excel 2016、Matlab 7.0和Canoco 5.0等軟件對試驗數(shù)據(jù)結(jié)果進(jìn)行統(tǒng)計、分析和作圖。其中Excel 2016主要用于描述性統(tǒng)計分析,計算數(shù)據(jù)的平均值、標(biāo)準(zhǔn)差和變異系數(shù)等;Matlab 7.0用于灰度關(guān)聯(lián)分析;Canoco 5.0用于降趨對應(yīng)分析(Detrended correspondence analysis,DCA)和典范對應(yīng)分析[20-21]。

    2 結(jié)果與分析

    2.1巖溶峰叢山區(qū)土壤-土壤水地球化學(xué)性質(zhì)的空間變化特征

    表1為巖溶峰叢山區(qū)土壤性質(zhì)的空間描述性統(tǒng)計分析,分析結(jié)果表明:土壤pH值的變異系數(shù)小于10.5%,遠(yuǎn)低于土壤的其他地球化學(xué)指標(biāo),說明其空間變化較?。幌鄬Χ?,土壤速效磷含量和有效鋅含

    量的變異系數(shù)在0~20 cm和20~40 cm均較大,其空間變化較大;全磷含量、pH值、有效鋅含量、有效鐵含量、有效硼含量、交換性鈣含量在0~20 cm的含量的變異系數(shù)比20~40 cm的大,而全氮含量、全鉀含量、速效氮含量、速效磷含量、速效鉀含量、有機(jī)質(zhì)含量、有效銅含量、有效錳含量、交換性鎂含量的變異系數(shù)則表現(xiàn)相反。從巖溶峰叢山區(qū)土壤水的地球化學(xué)指標(biāo)分析結(jié)果可知(表2),土壤水pH值的變異系數(shù)最小,屬于弱變異,而鉀離子含量、鎂離子含量的變異系數(shù)均大于100%;除了鎂離子含量和氯離子含量的變異系數(shù)表現(xiàn)出0~20 cm小于20~40 cm,土壤水的其他地球化學(xué)指標(biāo)含量的變異系數(shù)均表現(xiàn)為0~20 cm大于20~40 cm。由此說明巖溶峰叢山區(qū)不同土地利用條件與不同深度土壤-土壤水地球化學(xué)性質(zhì)均存在不同程度的異質(zhì)性,各地球化學(xué)指標(biāo)值對巖溶峰叢山區(qū)的環(huán)境條件和人為活動影響的響應(yīng)程度不盡相同。

    2.2巖溶峰叢山區(qū)土壤水地球化學(xué)特征的灰色關(guān)聯(lián)分析

    表1 巖溶峰叢山區(qū)土壤地球化學(xué)的空間描述性統(tǒng)計特征Tab.1 Descriptive statistical analysis of soil geochemical indicators in karst peak-cluster area

    表2 巖溶峰叢山區(qū)土壤水地球化學(xué)的空間描述性統(tǒng)計特征Tab.2 Descriptive statistical analysis of soil water geochemical indicators in karst peak-cluster area

    表3 巖溶峰叢山區(qū)土壤水地球化學(xué)指標(biāo)的灰色關(guān)聯(lián)系數(shù)Tab.3 Grey correlation coefficient analysis of soil water geochemical indicators in karst peak-cluster area

    2.3 巖溶峰叢山區(qū)土壤環(huán)境質(zhì)量的灰度關(guān)聯(lián)分析

    以土壤全氮含量(TN)、全鉀含量(TK)、全磷含量(TP)、速效氮含量(AN)、速效磷含量(AP)、速效鉀含量(AK)、有機(jī)質(zhì)含量(OM)、pH值和土壤腐殖質(zhì)中的胡敏酸含量(HA)、富里酸含量(FA)、胡敏素含量(HM)等指標(biāo)作為綜合評價土壤環(huán)境質(zhì)量的指標(biāo)。由表4的灰度關(guān)聯(lián)分析結(jié)果可知,峰叢山區(qū)土壤地球化學(xué)指標(biāo)對土壤環(huán)境質(zhì)量的影響各不相同;0~20 cm淺層土壤各地球化學(xué)指標(biāo)對土壤環(huán)境質(zhì)量的影響能力由大到小依次為TN、HA、AK、FA、TK、pH值、OM、HM、TP、AN、AP,而20~40 cm深度指標(biāo)影響由大到小排序為pH值、TK、TN、HA、TP、AK、FA、AP、HM、AN、OM;0~20 cm淺層土壤各地球化學(xué)指標(biāo)對土壤質(zhì)量的影響程度相差較大,其中TN、HA、AK、FA、TK、pH值對土壤質(zhì)量的影響較大,相對而言,除了pH值和TK含量的影響較大外,20~40 cm深度土壤其他各地球化學(xué)指標(biāo)對土壤質(zhì)量的影響程度較為接近,說明隨著土壤深度的增加,土壤環(huán)境質(zhì)量受單一地球化學(xué)指標(biāo)的影響不突出,而是趨向于受土壤地球化學(xué)指標(biāo)的綜合影響。

    表4 巖溶峰叢山區(qū)土壤環(huán)境質(zhì)量的灰度關(guān)聯(lián)系數(shù)Tab.4 Grey correlation coefficient of soil geochemical indicators in karst peak-cluster area

    2.4土壤水地球化學(xué)特征對土壤環(huán)境因子響應(yīng)程度的典范對應(yīng)分析

    2.4.1排序模型的選擇及其可靠性分析

    先利用Canoco 5.0軟件對巖溶峰叢山區(qū)土壤水地球化學(xué)指標(biāo)進(jìn)行降趨對應(yīng)分析,DCA分析結(jié)果顯示4個排序軸中梯度最大值(即第一排序軸)大于4,說明單峰模型比線性模型更適合用于分析巖溶峰叢山區(qū)土壤水地球化學(xué)特征與土壤環(huán)境的響應(yīng)關(guān)系。分別將土壤養(yǎng)分(全氮、全磷、全鉀、速效氮、速效磷、速效鉀、有機(jī)質(zhì)的含量和pH值)、土壤中微量元素有效態(tài)和交換態(tài)含量(交換性鈣ECa、交換性鎂EMg、有效銅ACu、有效鐵AFe、有效鋅AZn、有效錳AMn、有效硼AB)、土壤機(jī)械和腐殖質(zhì)組成(0.2~2.0 mm粒徑、0.02~0.20 mm粒徑、0.002~0.020 mm粒徑、小于0.002 mm粒徑、土壤含水率、胡敏酸含量、富里酸含量、胡敏素含量)作為3類土壤環(huán)境因子,分別與土壤水地球化學(xué)指標(biāo)進(jìn)行典范對應(yīng)分析。由單峰模型中的約束性排序方法的排序圖判圖規(guī)則進(jìn)行相關(guān)關(guān)系的判斷。

    通過CCA排序分析討論土壤環(huán)境特征和土壤水地球化學(xué)特征之間的關(guān)系,其響應(yīng)關(guān)系研究結(jié)果如表5所示,典范對應(yīng)分析排序軸第一軸和第二軸的特征值變化范圍分別為0.361~0.381、0.092~0.100,土壤水地球化學(xué)指標(biāo)與土壤環(huán)境因子排序軸的相關(guān)系數(shù)均大于或等于0.940,用蒙特卡擬合(Monte Carlo permutation test)對其關(guān)系的顯著性進(jìn)行檢驗的結(jié)果值范圍為0.002~0.008,且第一和第二排序軸對土壤水地球化學(xué)性質(zhì)的解釋量分別大于等于51.50%和13.12%,對土壤與土壤水地球化學(xué)性質(zhì)之間關(guān)系的累計百分比變化率大于62.3%,第一軸和第二軸相關(guān)關(guān)系分別為趨近于零的0.020、0.013和0.040而表現(xiàn)出兩軸幾近于垂直的關(guān)系,以上分析結(jié)果說明典范對應(yīng)分析能可靠地反映土壤水的地球化學(xué)指標(biāo)值隨土壤地球化學(xué)環(huán)境因子的變化趨勢,可通過前2個排序軸的二維排序圖來分析討論土壤水地球化學(xué)特征對土壤環(huán)境因子的響應(yīng)。

    表5 土壤水地球化學(xué)性質(zhì)與土壤環(huán)境因子的CCA分析Tab.5 CCA results of soil water geochemical indicators and soil environmental factors

    2.4.2土壤養(yǎng)分環(huán)境因子對土壤水地球化學(xué)特征的影響

    由圖2中CCA排序分析的結(jié)果說明土壤養(yǎng)分的8個環(huán)境因子對土壤水地球化學(xué)性質(zhì)的空間分布特征都表現(xiàn)出了不同程度的影響,由箭頭連線長度可判斷出土壤養(yǎng)分的速效鉀對土壤水地球化學(xué)性質(zhì)的影響最大,全氮的影響最小。由線段箭頭所在的象限可知,土壤養(yǎng)分環(huán)境因子與土壤水地球化學(xué)特征的第一排序軸正相關(guān)的由大到小排序依次為全鉀含量、pH值、全氮含量、有機(jī)質(zhì)含量、速效氮含量,對應(yīng)的相關(guān)系數(shù)分別為0.45、0.42、0.31、0.30、0.18;負(fù)相關(guān)的由大到小排序依次為速效鉀含量、全磷含量、速效磷含量。由此可以得出土壤養(yǎng)分環(huán)境因子與土壤水地球化學(xué)性質(zhì)第一排序軸的相關(guān)性由大到小依次為速效鉀含量、全鉀含量、pH值、全磷含量、有機(jī)質(zhì)含量、速效氮含量、速效磷含量、全氮含量。土壤養(yǎng)分環(huán)境因子與土壤水地球化學(xué)性質(zhì)第二排序軸的相關(guān)性由大到小排序依次為速效鉀含量、全鉀含量、速效氮含量、速效磷含量、pH值、全氮含量、全磷含量、有機(jī)質(zhì)含量,由此說明速效鉀含量、全鉀含量、pH值和速效氮含量是影響土壤水地球化學(xué)特征的主要土壤養(yǎng)分因子,結(jié)合研究區(qū)施肥情況的走訪調(diào)查分析可知,受人工施用鉀肥影響下的土壤速效鉀和全鉀的含量對土壤水地球化學(xué)特征的影響最大。

    圖2 土壤水地球化學(xué)特征與土壤養(yǎng)分因子的CCA二維排序圖Fig.2 Two-dimensional sorting diagram of CCA for soil water geochemical indicators and soil nutrients environmental factors

    2.4.3土壤中微量元素含量與土壤水地球化學(xué)特征的相關(guān)聯(lián)分析

    由土壤中微量元素有效態(tài)和交換態(tài)含量與土壤水地球化學(xué)特征的CCA排序分析結(jié)果可知(圖3),土壤交換性鎂含量對土壤水地球化學(xué)特征的影響最大,有效鐵、錳、硼含量的影響非常弱,與土壤水地球化學(xué)特征的第一排序軸相關(guān)關(guān)系由大到小排序依次為EMg、AZn、ACu、ECa、AFe、AB、AMn,對應(yīng)的相關(guān)系數(shù)分別為0.83、0.54、-0.29、0.12、0.06、-0.04、-0.01。土壤中微量元素有效態(tài)與交換態(tài)含量與第二排序軸的相關(guān)系數(shù)整體比第一排序軸小,其由大到小排序依次為ACu、AFe、AZn、EMg、AMn、AB、ECa,其中有效硼和交換性鈣含量的影響非常小。根據(jù)以上分析可知土壤中微量元素中的交換性鎂、有效鋅、有效銅含量對土壤水地球化學(xué)特征的影響較大。

    圖3 土壤水地球化學(xué)特征與土壤中微量元素因子的CCA二維排序圖Fig.3 Two-dimensional sorting diagram of CCA for soil water geochemical indicators and soil medium and microelements environmental factors

    2.4.4土壤水地球化學(xué)特征對土壤機(jī)械與腐殖質(zhì)組成因子的響應(yīng)分析

    由于土壤腐殖質(zhì)具有較好的黏結(jié)性,是團(tuán)粒結(jié)構(gòu)形成的良好膠結(jié)劑,對土壤的結(jié)構(gòu)有一定的影響,因此將土壤腐殖質(zhì)的3種重要組成成分與土壤機(jī)械組成、土壤含水率歸類在一起討論它們與土壤水地球化學(xué)性質(zhì)的相關(guān)關(guān)系。分析結(jié)果顯示(圖4),土壤機(jī)械和腐殖質(zhì)組成與土壤水地球化學(xué)特征第一排序軸的相關(guān)關(guān)系由大到小依次為土壤含水率、0.2~2.0 mm粒徑、小于0.002 mm粒徑、0.002~0.020 mm粒徑、FA、0.02~0.20 mm粒徑、HM、HA,對應(yīng)的相關(guān)系數(shù)分別為-0.799、0.705、-0.568、0.441、0.419、0.364、0.244、-0.107,土壤含水率與第一排序軸的相關(guān)關(guān)系最大,呈顯著負(fù)相關(guān),這符合土壤水地球化學(xué)指標(biāo)值濃度隨土壤含水率增加而降低的客觀現(xiàn)象;土壤機(jī)械和腐殖質(zhì)組成與土壤水地球化學(xué)特征第二排序軸的相關(guān)關(guān)系由大到小依次為HA、HM、0.002~0.020 mm粒徑、小于0.002 mm粒徑、0.02~0.20 mm粒徑、FA、土壤含水率、0.2~2.0 mm粒徑,胡敏酸與第二排序軸呈顯著正相關(guān)關(guān)系,相關(guān)系數(shù)為0.736。相關(guān)系數(shù)大小與第一軸的差別不大,但是大小排序的前后指標(biāo)卻呈相反的方向,由此說明土壤機(jī)械和腐殖質(zhì)組成的各項指標(biāo)對土壤水地球化學(xué)特征指標(biāo)的影響比較復(fù)雜且具有多重影響性,因此需要對土壤含水率、0.2~2.0 mm粒徑、HA、HM、0.002~0.020 mm粒徑、小于0.002 mm粒徑6個指標(biāo)與土壤水地球化學(xué)性質(zhì)進(jìn)行適應(yīng)值分析。

    圖4 土壤水地球化學(xué)特征與土壤腐殖質(zhì)和機(jī)械組成因子的CCA二維排序圖Fig.4 Two-dimensional sorting diagram of CCA for soil water geochemical indicators and soil humus and mechanical composition factors

    3 討論

    CCA分析的結(jié)果顯示土壤水中Ca2+與影響土壤水地球化學(xué)特征的主要土壤環(huán)境因子之間的響應(yīng)程度較弱,這可能與富鈣的巖溶地質(zhì)背景條件有關(guān)[27-28],此外本文只綜合性的對土壤水與土壤環(huán)境因子進(jìn)行典范對應(yīng)分析,未細(xì)分樣地的坡度、坡向、土地利用等條件,這些問題有待于進(jìn)一步的研究。

    4 結(jié)論

    (1)巖溶峰叢山區(qū)土壤水地球化學(xué)性質(zhì)與土壤環(huán)境指標(biāo)均存在不同程度的空間異質(zhì)性,其中土壤水和土壤中pH值變異系數(shù)最小,能較好地保持巖溶地球化學(xué)背景條件下中性偏弱堿的酸堿條件;土壤水地球化學(xué)和土壤的其他指標(biāo)則受植被條件、人為活動等因素的影響而表現(xiàn)出不同的異質(zhì)性特征,而且除了鎂離子和氯離子,土壤水其他地球化學(xué)指標(biāo)的變異系數(shù)均表現(xiàn)為0~20 cm大于20~40 cm,說明表層土壤水地球化學(xué)性質(zhì)更易受環(huán)境條件的影響。

    1 蔣忠誠,羅為群,童立強(qiáng),等.21世紀(jì)西南巖溶石漠化演變特點(diǎn)及影響因素[J].中國巖溶,2016,35(5):461-468.

    JIANG Zhongcheng, LUO Weiqun, TONG Liqiang, et al. Evolution features of rocky desertification and influence factors in karst areas of southwest China in the 21St century[J]. Carsologica Sinica, 2016, 35(5):461-468. (in Chinese)

    2 JIANG Zhongcheng, LIAN Yanqing, QIN Xiaoqun. Rocky desertification in southwest China:impacts, causes, and restoration[J]. Earth-Science Reviews, 2014, 132(3):1-12.

    3 岳躍民,王克林,張偉,等.基于典范對應(yīng)分析的喀斯特峰叢洼地土壤-環(huán)境關(guān)系研究[J].環(huán)境科學(xué),2008,29(5):1400-1405.

    YUE Yuemin, WANG Kelin, ZHANG Wei, et al. Relationships between soil and environment in peak-cluster depression areas of karst region based on canonical correspondence analysis[J]. Environmental Science, 2008, 29(5): 1400-1405. (in Chinese)

    4 楊奇勇,張發(fā)旺.西南巖溶盆地土壤干容重協(xié)同克里格分析[J/OL].農(nóng)業(yè)機(jī)械學(xué)報,2015,46(2):126-131.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20150219amp;flag=1.DOI:10.6041/j.issn.1000-1298.2015.02.019.

    YANG Qiyong, ZHANG Fawang.Cokriging analysis on soil bulk density in karst basin of southwest China[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(2): 126-131. (in Chinese)

    5 李文軍,楊奇勇,彭保發(fā),等.西南巖溶區(qū)土壤全氮含量的空間變異分析[J/OL].農(nóng)業(yè)機(jī)械學(xué)報,2014,45(9):150-154.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20140925amp;flag=1.DOI:10.6041/j.issn.1000-1298.2014.09.025.

    LI Wenjun, YANG Qiyong,PENG Baofa, et al. Spatial variability of total soil nitrogen content in karst area of southwest China[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(9): 150-154. (in Chinese)

    6 劉月嬌,張洋,倪九派,等.基于地統(tǒng)計學(xué)煙區(qū)土壤養(yǎng)分及pH值空間異質(zhì)性分析[J].中國巖溶,2014,33(3):319-325.

    LIU Yuejiao, ZHANG Yang, NI Jiupai, et al. Spatial heterogeneity of soil nutrients and pH in tobacco area based on geostatistics: a case study in a karst area of Youyang county, Chongqing[J]. Carsologica Sinica,2014, 33(3): 319-325. (in Chinese)

    7 吳敏,劉淑娟,葉瑩瑩,等.喀斯特地區(qū)坡耕地與退耕地土壤有機(jī)碳空間異質(zhì)性及其影響因素[J].生態(tài)學(xué)報,2016,36(6):1619-1627.

    WU Min, LIU Shujuan, YE Yingying, et al. Spatial variability of surface soil organic carbon and its influencing factors in cultivated slopes and abandoned lands in a karst peak-cluster depression area[J]. Acta Ecologica Sinica, 2016, 36(6): 1619-1627. (in Chinese)

    8 徐慧芳,宋同清,黃國勤,等.喀斯特峰叢洼地區(qū)坡地不同土地利用方式下土壤水分的時空變異特征[J].生態(tài)學(xué)報,2014,34(18):5311-5319.

    XU Huifang, SONG Tongqing, HUANG Guoqin, et al. Spatiotemporal variation of soil moisture under different land use types in a typical karst hill region[J]. Acta Ecologica Sinica, 2014,34(18): 5311-5319. (in Chinese)

    9 張川,陳洪松, 張偉,等.喀斯特坡面表層土壤含水量、容重和飽和導(dǎo)水率的空間變異特征[J].應(yīng)用生態(tài)學(xué)報,2014,25(6):1585-1591.

    ZHANG Chuan, CHEN Hongsong, ZHANG Wei, et al. Spatial variation characteristics of surface soil water content, bulk density and saturated hydraulic conductivity on karst slopes [J]. Chinese Journal of Applied Ecology, 2014, 25(6): 1585-1591. (in Chinese)

    10 鄧艷,曹建華,蔣忠誠,等.西南巖溶石漠化綜合治理水-土-植被關(guān)鍵技術(shù)進(jìn)展與建議[J].中國巖溶,2016,35(5):476-485.

    DENG Yan, CAO Jianhua, JIANG Zhongcheng, et al. Advancement in key technologies for comprehensive treatment of water, soil and vegetation resources in karst rocky desertification areas[J]. Carsologica Sinica, 2016, 35(5): 476-485. (in Chinese)

    11 蔣忠誠,羅為群,鄧艷,等.巖溶峰叢洼地水土漏失及防治研究[J].地球?qū)W報,2014,35(5):535-542.

    JIANG Zhongcheng, LUO Weiqun, DENG Yan, et al. The leakage of water and soil in the karst peak cluster depression and its prevention and treatment[J].Acta Geoscientica Sinica, 2014, 35(5): 535-542. (in Chinese)

    12 覃星銘,蔣忠誠,馬祖陸,等.巖溶峰叢坡面土壤養(yǎng)分與水化學(xué)特征研究[J].中國水土保持,2012(3):26-30.

    QIN Xingming, JIANG Zhongcheng, MA Zulu, et al. Research on the slope soil nutrient and water chemical characteristics of karst peak-cluster[J].Soil and Water Conservation in China, 2012(3): 26-30. (in Chinese)

    13 吳華英,覃星銘,蔣忠誠,等.桂林會仙巖溶作用對石灰土壤水鋅離子濃度的影響[J].農(nóng)業(yè)現(xiàn)代化研究,2009,30(6):761-765.

    WU Huaying, QIN Xingming, JIANG Zhongcheng, et al. Effect of karstification on concentration of Zn2+in lime soil water at Huixian Guilin[J]. Research of Agricultural Modernization, 2009, 30(6):761-765. (in Chinese)

    14 王宇飛,趙秀蘭,何丙輝,等.漢豐湖夏季浮游植物群落與環(huán)境因子的典范對應(yīng)分析[J].環(huán)境科學(xué),2015,36(3): 922-927.

    WANG Yufei, ZHAO Xiulan,HE Binghui, et al. Canonical correspondence analysis of summer phytoplankton community and its environmental factors in Hanfeng Lake[J]. Environmental Science, 2015, 36(3): 922-927. (in Chinese)

    15 張瑩,劉暢,宋昂,等.基于典范對應(yīng)分析的會仙巖溶濕地土壤理化性質(zhì)與土壤酶活性關(guān)系研究[J].中國巖溶,2016,35(1):11-18.

    ZHANG Ying, LIU Chang, SONG Ang, et al. Relationship between soil physicochemical properties and soil enzyme activities in Huixian karst wetland system based on canonical correspondence analysis[J]. Carsologica Sinica, 2016, 35(1): 11-18. (in Chinese)

    16 勞文科,祁曉凡,劉慧敏,等.廣西果化龍何地區(qū)表層帶巖溶水系統(tǒng)及其水資源特征[J].中國巖溶,2008,27(2):122-128.

    LAO Wenke, QI Xiaofan, LIU Huimin, et al. Characters of epikarst water system and water resources in Longhe, Guohua, Guangxi[J]. Carsologica Sinica, 2008,27(2):122-128. (in Chinese)

    17 鮑士旦.土壤農(nóng)化分析[M].3版. 北京:中國農(nóng)業(yè)出版社,2000: 263-282.

    18 魯如坤.土壤化學(xué)分析[M].北京:中國農(nóng)業(yè)科技出版社,1999.

    19 中國地質(zhì)科學(xué)院巖溶地質(zhì)研究所.滲流土壤水收集器:中國,201120118981[P].2011-04-21.

    20 BRAAK T, CAJO J F. Canonical correspondence analysis: a new eigenvector technique for multivariate direct gradient analysis[J]. Ecology, 1986, 67(5):1167-1179.

    21 張金屯.?dāng)?shù)量生態(tài)學(xué)[M].北京:科學(xué)出版社,2004:157-164.

    22 張偉,劉淑娟,葉瑩瑩,等.典型喀斯特林地土壤養(yǎng)分空間變異的影響因素[J].農(nóng)業(yè)工程學(xué)報,2013,29(1):93-101.

    ZHANG Wei, LIU Shujuan, YE Yingying, et al. Spatial variability of soil nutrients and its influencing factors in typical karst virgin forest [J]. Transactions of the CSAE, 2013, 29(1):93-101. (in Chinese)

    23 胡忠良,潘根興,李戀卿,等. 貴州喀斯特山區(qū)不同植被下土壤C、N、P 含量和空間異質(zhì)性[J].生態(tài)學(xué)報,2009,29(8):4187-4195.

    HU Zhongliang, PAN Genxing, LI Lianqing, et al. Changes in pools and heterogeneity of soil organic carbon, nitrogen and phosphorus under different vegetation types in karst mountainous area of central Guizhou Province, China[J]. Acta Ecologica Sinica, 2009, 29(8):4187-4195. (in Chinese)

    24 朱曉鋒,陳洪松,付智勇.喀斯特灌叢坡地土壤-表層巖溶帶產(chǎn)流及氮素流失特征[J].應(yīng)用生態(tài)學(xué)報,2017,28(7):2197-2206.

    ZHU Xiaofeng, CHEN Hongsong, FU Zhiyong, et al. Runoff and nitrogen loss characteristics in soil-epikarst system on a karst shrub hillslope[J]. Chinese Journal of Applied Ecology, 2017, 28(7): 2197-2206. (in Chinese)

    25 袁道先.我國巖溶資源環(huán)境領(lǐng)域的創(chuàng)新問題[J].中國巖溶,2015,34(2):98-100.

    YUAN Daoxian. Scientific innovation in karst resources and environment research field of China[J]. Carsologica Sinica, 2015, 34(2):98-100. (in Chinese)

    26 袁道先,章程.巖溶動力學(xué)的理論探索與實踐[J].地球?qū)W報,2008,29(3):355-365.

    YUAN Daoxian, ZHANG Cheng. Karst dynamics theory in China and its practice[J]. Acta Geoscientica Sinica, 2008, 29(3): 355-365. (in Chinese)

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

    CAO Jianhua, YUAN Daoxian, PAN Genxing. Some soil features in karst ecosystem[J]. Advance in Earth Sciences, 2003, 18(1): 37-44. (in Chinese)

    28 李忠云,魏興琥,李保生,等.粵北巖溶峰叢自然坡面土壤鈣元素的空間分布[J].水土保持通報,2016,36(3):62-68.

    LI Zhongyun, WEI Xinghu, LI Baosheng, et al. Spatial distribution of soil calcium on natural slope of karst peak-cluster in northern Guangdong Province[J]. Bulletin of Soil and Water Conservation, 2016, 36(3): 62-68. (in Chinese)

    CharacteristicsofSoilWaterGeochemistryandItsResponsetoSoilEnvironmentalFactorsinKarstPeak-clusterArea

    QIN Xingming1,2HE Binghui1WU Huaying2

    (1.SchoolofResourcesandEnvironment,SouthwestUniversity,Chongqing400716,China2.InstituteofKarstGeology,ChineseAcademyofGeologicalSciences,Guilin541004,China)

    karst peak-cluster area; soil water; geochemistry; soil environmental factors; grey correlation; canonical correspondence analysis

    10.6041/j.issn.1000-1298.2017.11.034

    S153.5; X824

    A

    1000-1298(2017)11-0280-09

    2017-08-10

    2017-09-07

    國家重點(diǎn)研發(fā)計劃項目(2016YFC0502402)、國家自然科學(xué)基金項目(41502257)、中國地質(zhì)大調(diào)查項目(DD20160324)和中國地質(zhì)科學(xué)院基本科研業(yè)務(wù)費(fèi)項目(2017026)

    覃星銘(1983—),男,博士生,中國地質(zhì)科學(xué)院高級工程師,主要從事巖溶水土工程研究,E-mail: qxm212@karst.ac.cn

    何丙輝(1966—),男,教授,博士生導(dǎo)師,主要從事農(nóng)業(yè)水土工程研究,E-mail: hebinghui@swu.edu.cn

    猜你喜歡
    排序特征分析
    排序不等式
    隱蔽失效適航要求符合性驗證分析
    恐怖排序
    如何表達(dá)“特征”
    不忠誠的四個特征
    節(jié)日排序
    電力系統(tǒng)不平衡分析
    電子制作(2018年18期)2018-11-14 01:48:24
    刻舟求劍
    兒童繪本(2018年5期)2018-04-12 16:45:32
    抓住特征巧觀察
    電力系統(tǒng)及其自動化發(fā)展趨勢分析
    久久精品人人爽人人爽视色| 日韩欧美一区二区三区在线观看| 男女床上黄色一级片免费看| 18禁美女被吸乳视频| 精品熟女少妇八av免费久了| 99国产精品99久久久久| 亚洲成人久久性| 国产蜜桃级精品一区二区三区| 国产av精品麻豆| 欧美日韩福利视频一区二区| av国产精品久久久久影院| 老司机亚洲免费影院| 亚洲少妇的诱惑av| 俄罗斯特黄特色一大片| 操美女的视频在线观看| 人人澡人人妻人| 每晚都被弄得嗷嗷叫到高潮| 88av欧美| 欧美成狂野欧美在线观看| 国产人伦9x9x在线观看| 久久香蕉精品热| 日韩欧美一区视频在线观看| 成在线人永久免费视频| 国产亚洲欧美98| 黄色怎么调成土黄色| 在线观看66精品国产| 亚洲精品久久成人aⅴ小说| 久久精品国产亚洲av高清一级| 两人在一起打扑克的视频| 久久久久久久久中文| 99久久久亚洲精品蜜臀av| 在线免费观看的www视频| 最近最新中文字幕大全免费视频| 久久久国产一区二区| 免费在线观看视频国产中文字幕亚洲| 亚洲自拍偷在线| 大陆偷拍与自拍| 人人妻人人澡人人看| 啦啦啦在线免费观看视频4| 动漫黄色视频在线观看| 19禁男女啪啪无遮挡网站| 国内毛片毛片毛片毛片毛片| 一边摸一边抽搐一进一小说| 亚洲精品国产色婷婷电影| xxx96com| 色哟哟哟哟哟哟| 88av欧美| 视频区图区小说| 可以在线观看毛片的网站| 免费在线观看影片大全网站| 国产一区二区三区在线臀色熟女 | 欧美激情久久久久久爽电影 | а√天堂www在线а√下载| 成人国语在线视频| 97人妻天天添夜夜摸| 色在线成人网| 亚洲国产精品一区二区三区在线| 欧洲精品卡2卡3卡4卡5卡区| 在线观看免费视频日本深夜| 欧美黑人欧美精品刺激| 美女 人体艺术 gogo| 亚洲 欧美一区二区三区| 亚洲熟女毛片儿| 91成人精品电影| 国产av精品麻豆| 亚洲一区高清亚洲精品| 亚洲欧美精品综合一区二区三区| 国产精品一区二区三区四区久久 | 国产精品二区激情视频| 欧美色视频一区免费| 亚洲在线自拍视频| 少妇粗大呻吟视频| 黑丝袜美女国产一区| 很黄的视频免费| 91九色精品人成在线观看| 色在线成人网| 日日干狠狠操夜夜爽| 中出人妻视频一区二区| 欧美不卡视频在线免费观看 | 国产激情欧美一区二区| 乱人伦中国视频| 俄罗斯特黄特色一大片| 国内久久婷婷六月综合欲色啪| 在线观看66精品国产| 国产91精品成人一区二区三区| 嫁个100分男人电影在线观看| 欧美久久黑人一区二区| 嫩草影视91久久| 黑丝袜美女国产一区| 久久久久久亚洲精品国产蜜桃av| a级毛片黄视频| 丰满的人妻完整版| 中文字幕人妻熟女乱码| 精品国产国语对白av| 欧美日本亚洲视频在线播放| 18美女黄网站色大片免费观看| 久久久久精品国产欧美久久久| 黑人巨大精品欧美一区二区mp4| av视频免费观看在线观看| www.www免费av| 我的亚洲天堂| 男人舔女人的私密视频| 女性生殖器流出的白浆| 久久久国产一区二区| a级毛片黄视频| av欧美777| 国产精品 国内视频| 国产精品 欧美亚洲| 国产精品亚洲一级av第二区| 一个人免费在线观看的高清视频| 日本wwww免费看| 国产成人精品久久二区二区免费| 亚洲伊人色综图| 精品国产乱码久久久久久男人| 国产亚洲欧美精品永久| 亚洲av五月六月丁香网| 亚洲av美国av| 久久精品亚洲av国产电影网| 日韩免费高清中文字幕av| av天堂在线播放| 欧美激情高清一区二区三区| 天堂影院成人在线观看| 99国产极品粉嫩在线观看| 天堂中文最新版在线下载| www.精华液| 夜夜看夜夜爽夜夜摸 | www国产在线视频色| 午夜激情av网站| 欧美激情 高清一区二区三区| 别揉我奶头~嗯~啊~动态视频| 欧美日韩国产mv在线观看视频| 巨乳人妻的诱惑在线观看| 国产精品 国内视频| 久久久久久久久免费视频了| 自线自在国产av| 中文字幕另类日韩欧美亚洲嫩草| 制服人妻中文乱码| 亚洲成人久久性| 黑丝袜美女国产一区| 久久久国产一区二区| 亚洲中文字幕日韩| 亚洲九九香蕉| 国产高清国产精品国产三级| 国产欧美日韩一区二区三区在线| 99国产精品一区二区蜜桃av| 国产av一区在线观看免费| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利一区二区在线看| 成在线人永久免费视频| 丰满人妻熟妇乱又伦精品不卡| 国产成人欧美| 国产成人精品久久二区二区91| 国产人伦9x9x在线观看| 欧美日本亚洲视频在线播放| 嫁个100分男人电影在线观看| 中文字幕高清在线视频| 两个人免费观看高清视频| 青草久久国产| 国产亚洲精品综合一区在线观看 | 人人澡人人妻人| 久久久久亚洲av毛片大全| 大型av网站在线播放| 亚洲情色 制服丝袜| 欧美成人午夜精品| 国产亚洲精品久久久久5区| 国产高清videossex| 12—13女人毛片做爰片一| 国内毛片毛片毛片毛片毛片| 成人18禁在线播放| 亚洲少妇的诱惑av| 如日韩欧美国产精品一区二区三区| 一级毛片女人18水好多| 欧美在线黄色| 久久久精品国产亚洲av高清涩受| 一级a爱视频在线免费观看| 成年人免费黄色播放视频| 99久久国产精品久久久| 国产亚洲av高清不卡| 女性生殖器流出的白浆| 国产激情久久老熟女| 国产成人欧美| 色精品久久人妻99蜜桃| 免费在线观看日本一区| 中文亚洲av片在线观看爽| 国产91精品成人一区二区三区| 十八禁人妻一区二区| 亚洲黑人精品在线| 涩涩av久久男人的天堂| 香蕉丝袜av| 欧美日韩福利视频一区二区| 久久国产亚洲av麻豆专区| 久久精品91蜜桃| 欧美老熟妇乱子伦牲交| 18禁裸乳无遮挡免费网站照片 | 夜夜看夜夜爽夜夜摸 | 欧美丝袜亚洲另类 | 免费日韩欧美在线观看| 超碰97精品在线观看| 悠悠久久av| 男女午夜视频在线观看| 99国产精品免费福利视频| 一级作爱视频免费观看| 日日爽夜夜爽网站| 身体一侧抽搐| 国产欧美日韩精品亚洲av| 嫁个100分男人电影在线观看| 精品久久蜜臀av无| 精品国产国语对白av| 18禁观看日本| 啦啦啦 在线观看视频| 桃红色精品国产亚洲av| 级片在线观看| 女同久久另类99精品国产91| 国产精品久久久人人做人人爽| 国产成人精品无人区| 色综合欧美亚洲国产小说| 国产精品久久久久成人av| 精品久久久久久,| 俄罗斯特黄特色一大片| 夜夜夜夜夜久久久久| 好看av亚洲va欧美ⅴa在| 亚洲少妇的诱惑av| 少妇裸体淫交视频免费看高清 | 亚洲avbb在线观看| av天堂在线播放| 午夜福利在线观看吧| 日日摸夜夜添夜夜添小说| 精品一区二区三卡| 很黄的视频免费| 女同久久另类99精品国产91| 亚洲av第一区精品v没综合| 精品一区二区三区av网在线观看| 国产精品香港三级国产av潘金莲| 女警被强在线播放| 午夜两性在线视频| 欧美日韩亚洲综合一区二区三区_| 国产高清videossex| 午夜免费成人在线视频| 69精品国产乱码久久久| 9热在线视频观看99| 国产午夜精品久久久久久| 法律面前人人平等表现在哪些方面| 大型黄色视频在线免费观看| 免费看十八禁软件| 国产精品久久久久成人av| 99精品欧美一区二区三区四区| 欧美成人免费av一区二区三区| 欧美日韩亚洲高清精品| 最近最新免费中文字幕在线| 欧美人与性动交α欧美软件| 久久精品人人爽人人爽视色| 9色porny在线观看| 两个人免费观看高清视频| bbb黄色大片| 免费日韩欧美在线观看| 91成年电影在线观看| 国产极品粉嫩免费观看在线| 一级毛片高清免费大全| 精品人妻在线不人妻| 欧美在线一区亚洲| xxxhd国产人妻xxx| 国产精品国产av在线观看| 精品国产一区二区久久| 国产无遮挡羞羞视频在线观看| 日本a在线网址| 可以在线观看毛片的网站| 啦啦啦免费观看视频1| 在线观看日韩欧美| 天堂俺去俺来也www色官网| 免费女性裸体啪啪无遮挡网站| 亚洲欧美日韩高清在线视频| 一区二区日韩欧美中文字幕| 色综合欧美亚洲国产小说| 神马国产精品三级电影在线观看 | 精品电影一区二区在线| 久久精品亚洲av国产电影网| 欧美日韩av久久| 亚洲人成77777在线视频| 亚洲欧美日韩高清在线视频| 18禁国产床啪视频网站| 精品午夜福利视频在线观看一区| 国产视频一区二区在线看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲午夜理论影院| 视频区欧美日本亚洲| 国产av一区在线观看免费| av天堂在线播放| 一区在线观看完整版| 国产亚洲av高清不卡| 在线看a的网站| 午夜免费鲁丝| 免费在线观看完整版高清| 一级片'在线观看视频| 嫩草影视91久久| 91国产中文字幕| 国产精品1区2区在线观看.| 香蕉丝袜av| 成人精品一区二区免费| 女性生殖器流出的白浆| 看免费av毛片| 天堂√8在线中文| 欧美大码av| 欧美黑人精品巨大| 亚洲男人的天堂狠狠| 最好的美女福利视频网| 日本撒尿小便嘘嘘汇集6| 亚洲国产欧美日韩在线播放| 国产av又大| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品中文字幕一二三四区| 成人精品一区二区免费| 啦啦啦 在线观看视频| 久久人人爽av亚洲精品天堂| 国产蜜桃级精品一区二区三区| 国产精品亚洲av一区麻豆| 伊人久久大香线蕉亚洲五| 久久午夜亚洲精品久久| 欧美日韩黄片免| 在线观看午夜福利视频| 亚洲欧美日韩另类电影网站| 久久精品国产亚洲av高清一级| 欧美日韩国产mv在线观看视频| 亚洲熟妇熟女久久| 99精品在免费线老司机午夜| 宅男免费午夜| 久久热在线av| 美女高潮喷水抽搐中文字幕| 精品第一国产精品| 午夜福利,免费看| 午夜成年电影在线免费观看| 一区二区日韩欧美中文字幕| 在线播放国产精品三级| 欧美午夜高清在线| 国产高清国产精品国产三级| 欧美激情 高清一区二区三区| 国产片内射在线| 国产真人三级小视频在线观看| 国产成人一区二区三区免费视频网站| 国产深夜福利视频在线观看| 久久这里只有精品19| 精品电影一区二区在线| 精品一区二区三区四区五区乱码| 日韩欧美国产一区二区入口| 两人在一起打扑克的视频| 神马国产精品三级电影在线观看 | 欧美一区二区精品小视频在线| 高清毛片免费观看视频网站 | 一个人免费在线观看的高清视频| 日日摸夜夜添夜夜添小说| 天堂√8在线中文| 亚洲欧美日韩高清在线视频| 19禁男女啪啪无遮挡网站| 免费观看精品视频网站| 19禁男女啪啪无遮挡网站| 看免费av毛片| 欧美激情极品国产一区二区三区| 脱女人内裤的视频| 视频区欧美日本亚洲| 美女福利国产在线| 国产精品自产拍在线观看55亚洲| 咕卡用的链子| 亚洲精品中文字幕一二三四区| 免费在线观看视频国产中文字幕亚洲| 高清黄色对白视频在线免费看| 热99国产精品久久久久久7| 国产伦人伦偷精品视频| 亚洲自偷自拍图片 自拍| 久久中文字幕一级| 老司机亚洲免费影院| 日日干狠狠操夜夜爽| 欧美精品啪啪一区二区三区| cao死你这个sao货| 一二三四在线观看免费中文在| 69av精品久久久久久| 久久久久国产精品人妻aⅴ院| 欧美成人性av电影在线观看| 国产精品1区2区在线观看.| av超薄肉色丝袜交足视频| 午夜激情av网站| 久久精品91无色码中文字幕| 国产亚洲精品久久久久久毛片| 桃红色精品国产亚洲av| 免费在线观看影片大全网站| 日本wwww免费看| 亚洲少妇的诱惑av| 精品午夜福利视频在线观看一区| 久久精品国产亚洲av香蕉五月| 法律面前人人平等表现在哪些方面| 两个人免费观看高清视频| 国产精品永久免费网站| 亚洲成人国产一区在线观看| 9色porny在线观看| 亚洲精品国产精品久久久不卡| 久久国产亚洲av麻豆专区| 国产成人精品在线电影| 亚洲自偷自拍图片 自拍| а√天堂www在线а√下载| 色老头精品视频在线观看| 一级片免费观看大全| 麻豆成人av在线观看| 男人的好看免费观看在线视频 | 热99re8久久精品国产| 国产精品影院久久| 精品欧美一区二区三区在线| 国产av一区二区精品久久| 成人影院久久| 国产黄a三级三级三级人| 久久 成人 亚洲| 久久久精品欧美日韩精品| 成人永久免费在线观看视频| 人成视频在线观看免费观看| 一二三四在线观看免费中文在| 国产成人影院久久av| 成人特级黄色片久久久久久久| 真人一进一出gif抽搐免费| 女人被躁到高潮嗷嗷叫费观| 操美女的视频在线观看| 男女下面进入的视频免费午夜 | 成年人黄色毛片网站| 老司机靠b影院| 99久久国产精品久久久| 美女大奶头视频| 法律面前人人平等表现在哪些方面| 午夜免费激情av| 亚洲精品一卡2卡三卡4卡5卡| www.www免费av| 999精品在线视频| 中文字幕人妻熟女乱码| 免费看a级黄色片| 十分钟在线观看高清视频www| 日本欧美视频一区| av天堂久久9| 亚洲熟妇中文字幕五十中出 | 国产伦人伦偷精品视频| 变态另类成人亚洲欧美熟女 | 啦啦啦在线免费观看视频4| videosex国产| 亚洲va日本ⅴa欧美va伊人久久| 欧美日韩亚洲综合一区二区三区_| 看片在线看免费视频| 亚洲性夜色夜夜综合| 757午夜福利合集在线观看| 1024视频免费在线观看| 午夜视频精品福利| 两人在一起打扑克的视频| 午夜福利欧美成人| 久久久久久久久免费视频了| 一个人免费在线观看的高清视频| 色老头精品视频在线观看| 丝袜人妻中文字幕| 亚洲 国产 在线| 欧美黄色淫秽网站| 少妇被粗大的猛进出69影院| 日韩大尺度精品在线看网址 | 免费在线观看亚洲国产| 精品国产一区二区三区四区第35| 天堂影院成人在线观看| 国产一区二区在线av高清观看| 女人爽到高潮嗷嗷叫在线视频| 亚洲aⅴ乱码一区二区在线播放 | 日本wwww免费看| 国产精品久久久av美女十八| 国产高清videossex| 99在线人妻在线中文字幕| 免费在线观看日本一区| 国产精品av久久久久免费| 亚洲一区中文字幕在线| 国产色视频综合| 久久欧美精品欧美久久欧美| 国产精品久久电影中文字幕| 亚洲精品美女久久av网站| 99久久99久久久精品蜜桃| 69av精品久久久久久| 久久精品91蜜桃| 精品久久久久久电影网| 日韩欧美一区视频在线观看| 99精品在免费线老司机午夜| 51午夜福利影视在线观看| 国产亚洲av高清不卡| 国产欧美日韩综合在线一区二区| 国产xxxxx性猛交| 亚洲精品国产色婷婷电影| 欧美激情极品国产一区二区三区| 欧美国产精品va在线观看不卡| 久久久久九九精品影院| www.自偷自拍.com| 欧美中文综合在线视频| 国产一区二区在线av高清观看| 亚洲自偷自拍图片 自拍| 日韩高清综合在线| 日韩精品青青久久久久久| 精品午夜福利视频在线观看一区| 中文字幕人妻丝袜一区二区| 亚洲中文日韩欧美视频| 99国产精品一区二区蜜桃av| 免费在线观看影片大全网站| 久久久久国产一级毛片高清牌| 欧美激情 高清一区二区三区| 丝袜美足系列| 久久亚洲真实| 一进一出抽搐gif免费好疼 | 一a级毛片在线观看| 久久久久久久久免费视频了| 日韩欧美免费精品| 亚洲三区欧美一区| av片东京热男人的天堂| 午夜免费观看网址| 精品卡一卡二卡四卡免费| 电影成人av| 纯流量卡能插随身wifi吗| 欧美成人免费av一区二区三区| 88av欧美| 免费在线观看视频国产中文字幕亚洲| 中文欧美无线码| 女人高潮潮喷娇喘18禁视频| 搡老乐熟女国产| 露出奶头的视频| 国产精品1区2区在线观看.| 国产欧美日韩综合在线一区二区| 亚洲精品久久成人aⅴ小说| 日日夜夜操网爽| 久久人妻熟女aⅴ| 日本黄色视频三级网站网址| 首页视频小说图片口味搜索| 精品乱码久久久久久99久播| 欧美乱色亚洲激情| 久久午夜亚洲精品久久| 亚洲第一av免费看| 国产99久久九九免费精品| 欧美成人性av电影在线观看| 欧美日韩亚洲高清精品| 美女扒开内裤让男人捅视频| 国产一区二区在线av高清观看| 亚洲一区中文字幕在线| 日韩有码中文字幕| 一区二区日韩欧美中文字幕| 国产成人精品久久二区二区免费| 久久精品91无色码中文字幕| 午夜免费成人在线视频| 91大片在线观看| 欧美日韩亚洲综合一区二区三区_| 黄网站色视频无遮挡免费观看| 女人精品久久久久毛片| 在线免费观看的www视频| 久久精品亚洲精品国产色婷小说| 欧美精品一区二区免费开放| 国产精品久久久av美女十八| 免费观看精品视频网站| 亚洲黑人精品在线| 亚洲精品国产一区二区精华液| 亚洲熟妇熟女久久| 午夜激情av网站| 日本wwww免费看| 久久精品国产清高在天天线| av电影中文网址| 欧美成狂野欧美在线观看| 在线观看一区二区三区激情| 国产精品电影一区二区三区| cao死你这个sao货| 韩国精品一区二区三区| 国产人伦9x9x在线观看| 少妇裸体淫交视频免费看高清 | 十八禁人妻一区二区| 亚洲成人免费av在线播放| 50天的宝宝边吃奶边哭怎么回事| 国产成人精品久久二区二区91| 大型黄色视频在线免费观看| 级片在线观看| 国产亚洲精品久久久久5区| 黄色 视频免费看| 国产激情久久老熟女| 午夜福利在线免费观看网站| 久久精品91蜜桃| 欧美成人免费av一区二区三区| 国产三级黄色录像| 久久国产精品影院| 性欧美人与动物交配| 99在线人妻在线中文字幕| 国产精品99久久99久久久不卡| 国产精品秋霞免费鲁丝片| 国产精品免费一区二区三区在线| 热99国产精品久久久久久7| 久99久视频精品免费| av视频免费观看在线观看| 国产亚洲精品一区二区www| 后天国语完整版免费观看| 国产一区在线观看成人免费| 国产1区2区3区精品| 又大又爽又粗| 欧美av亚洲av综合av国产av| 日日爽夜夜爽网站| 久久久久国内视频| 亚洲国产毛片av蜜桃av| 亚洲 欧美一区二区三区| 亚洲一区二区三区欧美精品| 天天躁狠狠躁夜夜躁狠狠躁| 久久 成人 亚洲| 超色免费av| 亚洲国产看品久久| 黑人操中国人逼视频| 后天国语完整版免费观看| 男女做爰动态图高潮gif福利片 | 日本黄色视频三级网站网址| 精品午夜福利视频在线观看一区| 久久精品人人爽人人爽视色| 在线观看www视频免费| 国产精品久久电影中文字幕| av在线天堂中文字幕 | a在线观看视频网站| 欧美在线黄色| 99久久99久久久精品蜜桃| 成人黄色视频免费在线看| 在线永久观看黄色视频| 亚洲va日本ⅴa欧美va伊人久久| 热99re8久久精品国产|