譚 艷，吳承禎，①，洪 偉，陳建忠，肖應(yīng)忠，陳 燦

（1.福建省高校森林生態(tài)系統(tǒng)過(guò)程與經(jīng)營(yíng)重點(diǎn)實(shí)驗(yàn)室福建農(nóng)林大學(xué)桉樹(shù)研究中心，福建 福州 350002；2.福建省建陽(yáng)市林業(yè)局，福建 建陽(yáng) 354200）

近年來(lái)，作為重要的速生豐產(chǎn)林樹(shù)種之一，桉樹(shù)（Eucalyptus spp.）的生產(chǎn)和種植越來(lái)越受到各級(jí)政府和學(xué)者的重視。鄧恩桉（Eucalyptus dunnii Maiden）是一種抗寒、速生性桉樹(shù)種類(lèi)，其主要用途是作為紙漿原料［1］，目前已成為中國(guó)南方地區(qū)工業(yè)原料林基地的主要造林樹(shù)種。鄧恩桉為外來(lái)引種樹(shù)種，耐寒性強(qiáng)，目前對(duì)該種的引種繁育、栽培技術(shù)及耐寒能力等方面已有大量的研究報(bào)道［1-6］，且有學(xué)者探討了其林地土壤養(yǎng)分與土壤酶活性之間的關(guān)系［7］。鄧恩桉要求在土層深厚、水肥條件較好的立地條件下種植，否則生長(zhǎng)速率會(huì)受到一定影響，因而，土壤條件決定了鄧恩桉的生長(zhǎng)質(zhì)量，土壤養(yǎng)分狀況也成為鄧恩桉種植和培養(yǎng)過(guò)程中營(yíng)養(yǎng)調(diào)控的重要依據(jù)。土壤pH能影響土壤養(yǎng)分有效性的發(fā)揮，是土壤的重要屬性之一。因此，研究鄧恩桉林地土壤pH的空間分布特征，對(duì)于鄧恩桉的栽培管理及其種植地土壤生產(chǎn)力的提高具有重要意義。

目前，地統(tǒng)計(jì)學(xué)方法已被證明是研究土壤養(yǎng)分和其他屬性空間變異規(guī)律的最佳方法之一［8］。該方法可同時(shí)考慮到樣本值的大小和樣本的空間位置及樣本的間距，彌補(bǔ)了經(jīng)典統(tǒng)計(jì)學(xué)方法忽略空間方位的缺憾。近年來(lái)，地統(tǒng)計(jì)學(xué)方法的應(yīng)用領(lǐng)域不斷拓寬，從土壤到地質(zhì)、礦業(yè)，從景觀生態(tài)格局到生態(tài)學(xué)，從地下水水文過(guò)程到水資源等［9］，但應(yīng)用地統(tǒng)計(jì)學(xué)方法研究林地土壤養(yǎng)分與土壤pH的相關(guān)報(bào)道較少，而關(guān)于鄧恩桉林地土壤pH的空間變異分析研究尚未見(jiàn)報(bào)道。

作者以福建省建陽(yáng)市7年生鄧恩桉林地為研究對(duì)象，測(cè)定了林地不同土層（0～20、20～40和40～60 cm）土壤的 pH，運(yùn)用 ArcGIS9.2軟件的Geostatistic Analyst模塊對(duì)土壤pH的空間變異特征及其分布格局進(jìn)行了分析，以期掌握鄧恩桉林地土壤pH的空間變異規(guī)律，為鄧恩桉種植管理及林地長(zhǎng)期生產(chǎn)力的維護(hù)奠定研究基礎(chǔ)。

1 研究區(qū)概況和研究方法

1.1 研究區(qū)概況

建陽(yáng)市地處福建省北部、武夷山南麓，地理坐標(biāo)為東經(jīng)117°31′～118°38′、北緯27°06′～27°43′，總面積3378.2 km2，丘陵山地面積占總面積的90.4％；土壤類(lèi)型為紅壤和黃壤，土層深厚肥沃。該地區(qū)屬中亞熱帶季風(fēng)性氣候，光熱資源豐富；冬短夏長(zhǎng)、氣候宜人，靜風(fēng)多、溫差大，雨季集中；年平均氣溫18℃，年平均無(wú)霜期282 d，年平均降雨量1700～2400 mm，年均日照時(shí)數(shù)1802 h，適宜動(dòng)植物生長(zhǎng)繁衍。

建陽(yáng)市2004年開(kāi)始先后引種鄧恩桉3000 hm2。本研究對(duì)象為建陽(yáng)市潭城街道城關(guān)村7年生鄧恩桉林地，林地土壤為紅壤，造林面積6.6 hm2，種植密度為1665株·hm-2，株距和行距分別為2和3 m；鄧恩桉的平均胸徑11.9 cm、平均樹(shù)高13.1 m。林地前茬為馬尾松毛蟲(chóng)（Dendrolimus punctatus Walker）重度危害樹(shù)木枯死1 a后進(jìn)行采伐的采伐跡地，屬陽(yáng)坡地、坡向?yàn)槿?，但立地質(zhì)量較差，土層較厚，坡度在25°以下。

1.2 土壤樣品采集與處理及pH測(cè)定方法

于2011年4月中旬在研究區(qū)內(nèi)劃出南北向長(zhǎng)660 m、東西向長(zhǎng)90 m的采樣區(qū)，按照30 m等距劃分網(wǎng)格，在網(wǎng)格的4個(gè)頂點(diǎn)各采集1個(gè)土樣，理論上應(yīng)有采樣點(diǎn)92個(gè)，但是受地形和鄧恩桉分布狀況的影響，共設(shè)采樣點(diǎn)65個(gè)（詳見(jiàn)圖1），基本涵蓋了7年生鄧恩桉的主要栽植區(qū)域。65個(gè)采樣點(diǎn)的土壤均按0～20、20～40和40～60 cm分3層采集，共采集土壤樣品195個(gè)，自然風(fēng)干。

參照文獻(xiàn)［10］采用電位法測(cè)定土壤pH，測(cè)定時(shí)土和水的質(zhì)量比為1∶2.5，每個(gè)土樣重復(fù)測(cè)定3次，結(jié)果取平均值。

1.3 數(shù)據(jù)分析與處理

分別對(duì)鄧恩桉林地各層土壤的pH數(shù)據(jù)進(jìn)行描述性統(tǒng)計(jì)分析，統(tǒng)計(jì)特征值包括平均值、中位值、標(biāo)準(zhǔn)差、變異系數(shù)、偏度、峰度等。通過(guò)ArcGIS9.2軟件的地統(tǒng)計(jì)分析模塊Geostatistic Analyst對(duì)鄧恩桉林地各層土壤pH進(jìn)行半方差函數(shù)理論模型相關(guān)參數(shù)和誤差的檢驗(yàn)分析，從而選出最優(yōu)半方差函數(shù)理論模型。在此基礎(chǔ)上，通過(guò)ArcGIS9.2地統(tǒng)計(jì)軟件的普通克里格（Kriging）法對(duì)鄧恩桉林地各層土壤pH進(jìn)行空間局部插值分析和插值標(biāo)準(zhǔn)差分析，并以插值標(biāo)準(zhǔn)差分析評(píng)價(jià)空間插值效果。

圖1 鄧恩桉林地土樣采樣點(diǎn)分布圖Fig.1 Distribution diagram of soil sam p le p lots in Eucalyptus dunnii Maiden forest

2 結(jié)果和分析

2.1 鄧恩桉林地土壤pH的描述性統(tǒng)計(jì)分析

對(duì)7年生鄧恩桉林地各層土壤pH的描述性統(tǒng)計(jì)分析結(jié)果見(jiàn)表1。由表1可見(jiàn)：0～20 cm土層pH最小值為pH3.31，最大值為pH3.93；20～40 cm土層pH最小值為pH3.30，最大值為pH3.86；40～60 cm土層pH最小值為pH3.33，最大值為pH3.90。鄧恩桉林地3層土壤pH的變幅較小，均小于0.62，其平均值自上而下逐漸增大。分別對(duì)0～20和20～40 cm土層的pH以及20～40和40～60 cm土層的pH進(jìn)行配對(duì)樣本t檢驗(yàn)，檢驗(yàn)結(jié)果表明：雙尾顯著性水平均為0.000，表明該鄧恩桉林地0～20與20～40 cm土層的pH以及20～40和40～60 cm土層的pH均存在顯著性差異。

表1 7年生鄧恩桉林地不同土層土壤pH描述性統(tǒng)計(jì)分析結(jié)果Table1 Analysis results of descriptive statistics of pH in different soil layers of seven-year-old forest of Eucalyptus dunnii Maiden

鄧恩桉林地0～20、20～40和40～60 cm土層的土壤pH變異系數(shù)分別為3.90％、3.93％和4.05％，變異幅度很小，說(shuō)明該林地3層土壤的pH均屬于弱變異［11］，并呈現(xiàn)隨土層加深自上而下逐漸增大的趨勢(shì)。對(duì)該鄧恩桉林地3層土壤的pH分別進(jìn)行非參數(shù)單樣本K-S檢驗(yàn)，結(jié)果表明：0～20、20～40和40～60 cm土層土壤pH的P值分別為0.909、0.523和0.710，各土層土壤pH均服從正態(tài)分布，滿(mǎn)足地統(tǒng)計(jì)學(xué)分析要求。

2.2 鄧恩桉林地土壤pH的半方差分析

利用ArcGIS9.2軟件的Geostatistic Analyst模塊獲得7年生鄧恩桉林地各層土壤pH的空間變異半方差函數(shù)理論模型的相關(guān)參數(shù)和誤差檢驗(yàn)值（表2），并對(duì)4種半方差函數(shù)的理論模型（圓狀模型 Circular model、球狀模型Sphericalmodel、指數(shù)模型Exponential model和高斯模型Gaussian model）在0階、1階、2階、3階趨勢(shì)效應(yīng)下預(yù)測(cè)誤差的各指標(biāo)進(jìn)行比較，篩選出最優(yōu)半方差函數(shù)模型，最終分析確定該鄧恩桉林地0～20、20～40和40～60 cm土層土壤pH空間變異的最優(yōu)模型分別為0階高斯模型、1階指數(shù)模型和1階球狀模型（表2）。

半方差函數(shù)模型具有4個(gè)重要的參數(shù)，即塊金值（C0）、變程、基臺(tái)值（C0+C）和偏基臺(tái)值。其中，塊金值與基臺(tái)值的比值稱(chēng)為塊金效應(yīng)，用以表明土壤屬性的空間變異特征；根據(jù)該比值大小可將相關(guān)土壤屬性劃分為強(qiáng)空間相關(guān)、中等空間相關(guān)和弱空間相關(guān)3個(gè)等級(jí)［12］。供試的7年生鄧恩桉林地0～20、20～40和40～60 cm土層土壤pH的塊金效應(yīng)分別為0.3896、0.4303、0.6043（表3），均位于0.25～0.75之間，說(shuō)明該7年生鄧恩桉林地土壤pH的空間相關(guān)性中等。

表2 7年生鄧恩桉林地不同土層土壤pH的最優(yōu)半方差模型預(yù)測(cè)誤差Table2 Prediction errors of the best sem ivariogram model of pH in different soil layers of seven-year-old forest of Eucalyptus dunnii Maiden

表3 7年生鄧恩桉林地不同土層土壤pH最優(yōu)半方差模型參數(shù)Table3 Parameters of the best sem ivariogram model of pH in different soil layers of seven-year-old forest of Eucalyptus dunnii Maiden

長(zhǎng)軸變程和短軸變程分別是相關(guān)性消失最慢和最快方向的變程，經(jīng)計(jì)算，該7年生鄧恩桉林地不同土層土壤pH的長(zhǎng)軸變程為145.684～167.809 m，短軸變程為55.7835～85.1553 m，也即不同土層土壤pH空間相關(guān)性消失的范圍在55～168 m，而本研究的采樣距離為30 m，因此，該采樣尺度在變程以?xún)?nèi)，符合采樣距離要求。長(zhǎng)軸變程與短軸變程的比值稱(chēng)為各向異性比（Anisotropic ratio），其值接近或等于1時(shí)，空間異質(zhì)性為各向同性，反之為各向異性。由表3可見(jiàn)：該7年生鄧恩桉林地不同土層土壤pH的各向異性比均大于1，所以呈現(xiàn)明顯的各向異性，符合空間插值的要求。

2.3 鄧恩桉林地土壤pH的普通克里格插值及其標(biāo)準(zhǔn)差分析

由于該鄧恩桉林地各層土壤pH數(shù)據(jù)經(jīng)驗(yàn)證均符合正態(tài)分布，因此采用普通克里格法對(duì)7年生鄧恩桉林地不同土層土壤pH進(jìn)行了最優(yōu)插值分析，并繪制不同土層土壤pH的空間分布圖，結(jié)果見(jiàn)圖2。由圖2可見(jiàn)：該鄧恩桉林地不同土層土壤pH等值線(xiàn)整體上既體現(xiàn)相似性，又體現(xiàn)出一定的差異性。不同土層土壤pH分布均呈現(xiàn)出帶狀分布，且?guī)罘植嫉娘@著度隨土壤深度的增加而增強(qiáng)。不同土層土壤pH的空間變異格局均呈現(xiàn)由西向東遞減的趨勢(shì)；0～20、20～40和40～60 cm土層土壤pH分布格局相似，說(shuō)明各層土壤pH之間的空間相關(guān)性極強(qiáng)。0～20 cm土層土壤pH最高值集中在西南方向的斑塊區(qū)域，并向東北方向和東面逐漸減小，最小值出現(xiàn)在東北方向（圖2-a）；20～40和40～60 cm土層土壤pH最高值均出現(xiàn)在西南側(cè)和西北側(cè)，并向東面呈波浪帶狀遞減，且變化比較均勻，最低值均出現(xiàn)在東北方向（圖2-b，c）。該林地20～40與40～60 cm土層間土壤pH的分布格局相似度高于0～20與20～40 cm土層之間的相似度。

進(jìn)一步分析該鄧恩桉林地不同土層土壤pH經(jīng)普通克里格插值分析之后產(chǎn)生的插值標(biāo)準(zhǔn)差的空間分布，結(jié)果見(jiàn)圖3。由圖3可見(jiàn)：該林地0～20、20～40和40～60 cm土層土壤pH插值標(biāo)準(zhǔn)差的變化范圍分別為0.1071～0.1767、0.0957～0.1243、0.1079～0.1282，pH插值標(biāo)準(zhǔn)差較小，表明各土層土壤pH插值結(jié)果較好。由圖3還可見(jiàn)：該鄧恩桉林地不同土層土壤pH的插值標(biāo)準(zhǔn)差均呈現(xiàn)出一定的規(guī)律性，大體上以南部和中部為中心，向東北方向逐漸擴(kuò)散增大，在東北部達(dá)到最大值。0～20 cm土層土壤pH插值標(biāo)準(zhǔn)差的空間分布以低值帶為中心向四周增大（圖3-a）；20～40 cm土層土壤pH插值標(biāo)準(zhǔn)差的空間分布出現(xiàn)許多低值小斑塊，插值標(biāo)準(zhǔn)差圍繞低值小斑塊向四周增大，并將小斑塊連在一起，形成1個(gè)與采樣點(diǎn)布局相似的網(wǎng)格（圖3-b）。該林地20～40與40～60 cm土層間土壤pH插值標(biāo)準(zhǔn)差分布圖相似度較高（圖3-b，c）。

圖2 7年生鄧恩桉林地不同土層土壤pH空間分布圖Fig.2 Spatial distribution maps of pH in different soil layers of seven-year-old forest of Eucalyptus dunnii Maiden

圖3 7年生鄧恩桉林地不同土層土壤pH插值標(biāo)準(zhǔn)差空間分布圖Fig.3 Spatial distribution maps of interpolation standard deviation of pH in different soil layers of seven-year-old forest of Eucalyptus dunnii Maiden

3 結(jié)論和討論

經(jīng)統(tǒng)計(jì)分析，供試7年生鄧恩桉林地0～20、20～40和40～60 cm土層土壤pH均值分別為pH3.57、pH3.62和pH3.67，變化幅度分別為pH3.31～pH3.93、pH3.30～pH3.86和pH3.33～pH3.90，變異系數(shù)分別為3.90％、3.93％和4.05％，說(shuō)明隨土層深度的增加土壤的酸性減弱，pH的空間變異增強(qiáng)。經(jīng)K-S檢驗(yàn)，研究區(qū)域內(nèi)的土壤pH服從正態(tài)分布。陳朝陽(yáng)［12］的研究結(jié)果顯示：建陽(yáng)市植煙區(qū)土壤pH的變幅為pH3.96～pH7.0，且土壤有酸化的趨勢(shì)；而鄧恩桉林地的土壤pH小于植煙區(qū)的土壤pH，說(shuō)明鄧恩桉林區(qū)土壤酸化程度比植煙區(qū)更嚴(yán)重。導(dǎo)致鄧恩桉林地土壤酸化嚴(yán)重的因素可能有4個(gè)：一為酸雨的影響，建陽(yáng)市為福建省內(nèi)受酸雨影響最為嚴(yán)重的城市之一，酸雨的酸度小于pH4.5［13］；二為人工施肥的影響，研究區(qū)土壤在造林前與造林過(guò)程中經(jīng)過(guò)多次施肥；三為種植鄧恩桉的影響，有研究表明桉樹(shù)的種植會(huì)使土壤pH降低［14］；四為造林地前茬的影響，該鄧恩桉林地造林前茬為馬尾松（Pinus massoniana Lamb.）林，馬尾松生長(zhǎng)過(guò)程中從土壤中吸收大量營(yíng)養(yǎng)物質(zhì)，導(dǎo)致土壤酸化，若馬尾松死后枝體歸還土壤，將使土壤酸化狀況得到一定的緩解，但若將馬尾松從林地中大量移走，相當(dāng)于將土壤中的鹽基離子移走，馬尾松被移走的越多，土壤酸化效應(yīng)越明顯［15］。供試的鄧恩桉林地為馬尾松受蟲(chóng)害后1 a被人工采伐的林地，因馬尾松被采伐使大量鹽基離子無(wú)法歸還土壤，因此無(wú)法使土壤酸化狀況得到緩解。由于鄧恩桉林地的土壤pH較低，因此在輪作種植抗酸性弱或?qū)ν寥缐A性要求較高的植物時(shí)，應(yīng)注意對(duì)林地土壤pH的改良。

地統(tǒng)計(jì)學(xué)分析結(jié)果表明：鄧恩桉林地0～20、20～40和40～60 cm土層的土壤pH塊金效應(yīng)值分別為0.3896、0.4303和0.6043，均具有中等的空間相關(guān)性，這可能與林地土壤類(lèi)型和經(jīng)營(yíng)措施等因素有關(guān)。通?？臻g相關(guān)性越大，表明土壤屬性主要由結(jié)構(gòu)因素決定，如成土母質(zhì)、地形、土壤類(lèi)型等；反之則表明土壤屬性主要由隨機(jī)因素決定，如施肥等人為干擾活動(dòng)。鄧恩桉林地土壤pH的空間相關(guān)性處于中等水平，說(shuō)明該林地的土壤屬性是由結(jié)構(gòu)因素和隨機(jī)因素共同決定的。鄧恩桉林地土壤塊金效應(yīng)值自上而下逐漸變大，土壤pH的空間相關(guān)性越來(lái)越強(qiáng)，說(shuō)明隨土層深度的增加，土壤pH受隨機(jī)因素的影響不斷減小，而受結(jié)構(gòu)因素的影響不斷增大。具體而言，在該鄧恩桉林地范圍內(nèi)，結(jié)構(gòu)因素主要是地形等，隨機(jī)因素主要是施肥管理，即鄧恩桉林地各層土壤pH受土壤類(lèi)型、成土母質(zhì)、地形和人工施肥等因素的綜合影響，且自上而下受人為干擾的影響越來(lái)越小。

經(jīng)過(guò)對(duì)7年生鄧恩桉林地3層土壤pH的0、1、2、3階的4個(gè)半方差函數(shù)模型的預(yù)測(cè)誤差進(jìn)行比較，獲得的不同層次最優(yōu)半方差函數(shù)模型分別為：0～20 cm土層土壤pH適用0階高斯模型；20～40 cm土層土壤pH適用1階指數(shù)模型；40～60 cm土層土壤pH適用1階球狀模型。對(duì)7年生鄧恩桉林地土壤pH進(jìn)行普通克里格插值分析并繪制空間分布圖，結(jié)果表明：0～20、20～40和40～60 cm土層的土壤pH均由西向東遞減，最高值出現(xiàn)在西南部；而林地土壤pH的插值標(biāo)準(zhǔn)差呈現(xiàn)由南部和中部向周?chē)皷|北部逐漸增大的分布規(guī)律。趙汝?yáng)|等［16］認(rèn)為：離采樣點(diǎn)越近，插值標(biāo)準(zhǔn)差越小，反之則插值標(biāo)準(zhǔn)差越大越不可信。因此，在采樣點(diǎn)周?chē)牟逯禈?biāo)準(zhǔn)差較小，離采樣點(diǎn)較遠(yuǎn)的位置插值標(biāo)準(zhǔn)差較大。本研究中，鄧恩桉林地各層土壤pH插值標(biāo)準(zhǔn)差的空間分布圖的表征規(guī)律與趙汝?yáng)|等［16］的研究結(jié)果相一致。另外，插值標(biāo)準(zhǔn)差大小還與采樣點(diǎn)的數(shù)量及采樣點(diǎn)的設(shè)置有關(guān)，受地形影響在北部設(shè)置的采樣點(diǎn)較少，越往南部采樣點(diǎn)越多，因此中南部位置的插值標(biāo)準(zhǔn)差較小，精度較高。有關(guān)采樣點(diǎn)數(shù)量與采樣點(diǎn)設(shè)置對(duì)鄧恩桉林地土壤pH插值的影響及其機(jī)制還有待于進(jìn)一步研究。

［1］燕麗萍，夏 陽(yáng)，毛秀紅，等.鄧恩桉的組織培養(yǎng)［J］.林業(yè)科學(xué)，2011，47（5）：157-161.

［2］宗亦臣，鄭勇奇，常金財(cái).鄧恩桉國(guó)內(nèi)引種和繁育技術(shù)研究進(jìn)展［J］.中國(guó)農(nóng)學(xué)通報(bào)，2010，26（5）：82-86.

［3］范海蘭，洪 偉，陳建忠，等.鄧恩桉引種區(qū)極端氣溫的起伏型時(shí)間序列分析［J］.河南農(nóng)業(yè)大學(xué)學(xué)報(bào)，2008，42（4）：391-393.

［4］宋建英.鄧恩桉的組織培養(yǎng)和植株再生［J］.林業(yè)科學(xué)，2010，46（6）：138-145.

［5］郭祥泉，洪 偉，吳承禎，等.“t”檢驗(yàn)法則在閩北抗寒速生鄧恩桉優(yōu)株篩選的應(yīng)用［J］.應(yīng)用與環(huán)境生物學(xué)報(bào)，2009，15（3）：385-389.

［6］張寧南，吳 煒，徐大平，等.P肥對(duì)鄧恩桉生長(zhǎng)與葉片養(yǎng)分狀況的影響［J］.中南林業(yè)科技大學(xué)學(xué)報(bào)，2010，30（6）：32-37.

［7］陳少杰.鄧恩桉林地土壤酶活性與土壤養(yǎng)分關(guān)系研究［J］.安徽農(nóng)學(xué)通報(bào)，2009，15（15）：171-183.

［8］GOOVAERTS P.Geostatistics in soil science：state- of -the-art and perspectives［J］.Geoderma，1999，89（1/2）：1-45.

［9］郭懷成，周 豐，刀 谞.地統(tǒng)計(jì)方法學(xué)研究進(jìn)展［J］.地理研究，2008，27（5）：1191-1202.

［10］鮑士旦.土壤農(nóng)化分析［M］.3版.北京：中國(guó)農(nóng)業(yè)出版社，2000.

［11］王紹強(qiáng)，朱松麗，周成虎.中國(guó)土壤土層厚度的空間變異性特征［J］.地理研究，2001，20（2）：161-169.

［12］陳朝陽(yáng).南平市植煙土壤pH狀況及其與土壤有效養(yǎng)分的關(guān)系［J］.中國(guó)農(nóng)學(xué)通報(bào)，2011，27（5）：149-153.

［13］趙衛(wèi)紅.福建省城市酸性降水特征及變化趨勢(shì)［J］.環(huán)境科學(xué)與技術(shù)，2006（9）：41-43.

［14］理永霞，羅 微，茶正早，等.桉樹(shù)幼苗對(duì)土壤養(yǎng)分及pH的影響［J］.中國(guó)土壤與肥料，2008（4）：77-78.

［15］許中堅(jiān)，劉廣深，劉維屏.人為因素誘導(dǎo)下的紅壤酸化機(jī)制及其防治［J］.農(nóng)業(yè)環(huán)境保護(hù)，2002，21（2）：175-178.

［16］趙汝?yáng)|，孫焱鑫，王殿武，等.北京地區(qū)耕地土壤有機(jī)質(zhì)空間變異分析［J］.土壤通報(bào)，2010，41（3）：552-557.