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

覃星銘 何丙輝 吳華英

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

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

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

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

巖溶峰叢山區(qū)； 土壤水； 地球化學性質； 土壤環(huán)境因子； 灰度關聯(lián)； 典范對應分析

引言

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

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

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

1 材料與方法

1.1 研究區(qū)概況

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

1.2 樣品采集與分析

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

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

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

2 結果與分析

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

表1為巖溶峰叢山區(qū)土壤性質的空間描述性統(tǒng)計分析，分析結果表明：土壤pH值的變異系數(shù)小于10.5%，遠低于土壤的其他地球化學指標，說明其空間變化較小；相對而言，土壤速效磷含量和有效鋅含

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

2.4.4土壤水地球化學特征對土壤機械與腐殖質組成因子的響應分析

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

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

3 討論

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

4 結論

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

1 蔣忠誠，羅為群，童立強，等．21世紀西南巖溶石漠化演變特點及影響因素[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 岳躍民，王克林，張偉，等．基于典范對應分析的喀斯特峰叢洼地土壤-環(huán)境關系研究[J]．環(huán)境科學，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è)機械學報，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è)機械學報，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)計學煙區(qū)土壤養(yǎng)分及pH值空間異質性分析[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]．生態(tài)學報，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)學報，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 張川，陳洪松， 張偉，等．喀斯特坡面表層土壤含水量、容重和飽和導水率的空間變異特征[J]．應用生態(tài)學報，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 鄧艷，曹建華，蔣忠誠，等．西南巖溶石漠化綜合治理水-土-植被關鍵技術進展與建議[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]．地球學報，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)分與水化學特征研究[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)境因子的典范對應分析[J]．環(huán)境科學，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 張瑩，劉暢，宋昂，等．基于典范對應分析的會仙巖溶濕地土壤理化性質與土壤酶活性關系研究[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 魯如坤．土壤化學分析[M]．北京：中國農(nóng)業(yè)科技出版社，1999．

19 中國地質科學院巖溶地質研究所．滲流土壤水收集器：中國，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 張金屯．數(shù)量生態(tài)學[M]．北京：科學出版社，2004:157-164．

22 張偉，劉淑娟，葉瑩瑩，等．典型喀斯特林地土壤養(yǎng)分空間變異的影響因素[J]．農(nóng)業(yè)工程學報，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 含量和空間異質性[J]．生態(tài)學報，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]．應用生態(tài)學報，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)境領域的創(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 袁道先，章程．巖溶動力學的理論探索與實踐[J]．地球學報，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]．地球科學進展，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

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

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

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