楊陽(yáng),段海俠
(遼寧工程技術(shù)大學(xué)資源與環(huán)境工程學(xué)院,123000,遼寧阜新)
我國(guó)主要地帶性土壤顆粒組成的分形維數(shù)存在逐漸遞增的趨勢(shì),成土母質(zhì)對(duì)土壤顆粒組成的分形維數(shù)有較大影響,易風(fēng)化和風(fēng)化程度高的母質(zhì)上發(fā)育的土壤,其質(zhì)地較細(xì),相應(yīng)的顆粒組成的分形維數(shù)也較大。分形維數(shù)的大小反映土壤質(zhì)地中黏粒、粉粒和砂粒質(zhì)量分?jǐn)?shù)的變化,隨黏粒質(zhì)量分?jǐn)?shù)的增多分形維數(shù)增大,隨砂粒質(zhì)量分?jǐn)?shù)的增多分形維數(shù)減小。我國(guó)學(xué)者將土壤顆粒質(zhì)量分形研究大量運(yùn)用到不同類型的土壤[1]、不同土地利用方式、不同植被類型的研究中,取得了很大的進(jìn)展。這些研究揭示出土壤顆粒分形維數(shù)與土壤有機(jī)C、全N、土壤粉黏粒質(zhì)量分?jǐn)?shù)、總孔隙度和毛管孔隙度之間存在明顯的線性正相關(guān),而與土壤砂粒質(zhì)量分?jǐn)?shù)和非毛管孔隙度具有線性負(fù)相關(guān)關(guān)系,它能較好表征植物群落變化對(duì)土壤顆粒組成[2]和孔隙結(jié)構(gòu)[3]狀況的影響,用粒徑的數(shù)量分布表征的土壤分形特征也具有這個(gè)特點(diǎn)。在有植被覆蓋的土壤上,由于根系和腐殖質(zhì)的參與,加強(qiáng)了土壤的團(tuán)聚能力,進(jìn)而降低了土壤顆粒的分形維數(shù)[4]。由于根據(jù)不同粒徑顆粒的質(zhì)量分布獲得的土壤顆粒尺寸分布(particle size distribution,簡(jiǎn)稱PSD)分形維數(shù)是在假定土壤密度均一的前提下得出的,該假定顯然與實(shí)際情況不符,所以受到一些學(xué)者的質(zhì)疑;因此,研究者們提出了PSD 體積分布的分形方法。土壤顆粒體積分形維數(shù)和質(zhì)量分形維數(shù)均可表示土壤固有屬性,但體積分形維數(shù)可以克服質(zhì)量分形維數(shù)計(jì)算模型假設(shè)不同粒級(jí)的土壤顆粒具有相同的密度這一不足,更具合理性[5],特別是激光衍射(Laser Diffraction,LD)技術(shù)的應(yīng)用極大地推動(dòng)了體積PSD 分形研究的發(fā)展[6]。
Mandelbrot 創(chuàng)立的分形論已成為描述自然界中復(fù)雜和不規(guī)則幾何形體一個(gè)有效工具,分形論在地學(xué)和土壤制圖學(xué)研究領(lǐng)域也得到了一定應(yīng)用。土壤是具有分形特征的系統(tǒng),由于土壤內(nèi)部的物理、化學(xué)、生物等過程的相互影響以及各種地質(zhì)過程和人為措施的干擾,導(dǎo)致土壤在形態(tài)、結(jié)構(gòu)、功能等方面表現(xiàn)為復(fù)雜的自然體。土壤是具有分形特征的系統(tǒng),由形狀和大小各異的顆粒組成,表觀上反映出不規(guī)則的幾何形體。研究排土場(chǎng)不同土地利用類型下的土壤分形特征、同種土地利用類型不同層次土壤分形特征及與土壤基本性質(zhì)間的相關(guān)性有很大的實(shí)用價(jià)值和指導(dǎo)意義。筆者通過對(duì)海州露天礦排土場(chǎng)現(xiàn)場(chǎng)勘查、土壤采樣分析,了解排土場(chǎng)不同土地利用類型下土壤分形特征和同種土地利用類型不同層次土壤分形特征,利用分形理論研究不同立地土壤顆粒分形特征與植物種類和土壤養(yǎng)分狀況之間的關(guān)系,探討煤炭開發(fā)對(duì)該區(qū)土壤環(huán)境的影響程度,為保護(hù)礦區(qū)土壤生態(tài)條件和恢復(fù)土壤生態(tài)功能提供科學(xué)依據(jù)。
海州露天礦位于遼寧省南部,地理坐標(biāo)為E 121°26'、N 42°02',海拔200 ~215 m。海州露天礦排土場(chǎng)位于阜新煤田的東部,東西長(zhǎng)7 km,南北寬3 km,占地16 km2,平面呈弧條帶,地勢(shì)東高西低,最高海拔320 m,排土場(chǎng)相對(duì)標(biāo)高32 ~200 m,邊坡的坡度為30°~42°,平臺(tái)較為平坦。土質(zhì)是由露天礦排出的頁(yè)巖、砂巖構(gòu)成,部分土地得到了良好的風(fēng)化。排土場(chǎng)是地帶性土壤在各種巖石風(fēng)化物殘積母質(zhì)上以黃土、紅土母質(zhì)發(fā)育的淋溶褐土,以褐土性土為代表,土層較薄。排土場(chǎng)基質(zhì)除地帶性土壤外,還夾雜著大量粉沙巖、礫巖、煤頁(yè)巖等巖石成分。
土壤樣品采集于2012 年5 月,采樣點(diǎn)設(shè)置在海州礦有代表性的3 個(gè)不同土地利用類型的樣區(qū)上,分別為海州礦排土年限為30 a 以上的人工榆樹(Ulmus pumila)林地、10 a 以上的人工草地及耕地。每個(gè)樣區(qū)按“Z”字形選取5 個(gè)樣點(diǎn),每個(gè)樣點(diǎn)10 cm 一個(gè)層次,一直到100 cm,部分樣點(diǎn)由于土石緊密,取樣不足100 cm,每個(gè)樣區(qū)的5 個(gè)樣點(diǎn)同一層次的土樣混合在一起,共計(jì)27 個(gè)土樣。
樣品采回后,剔除植物根系和石塊等雜質(zhì),自然風(fēng)干后用四分法將土樣磨碎,分別過1 mm 和0.5 mm 篩孔。取過1 mm 篩孔的樣品約5 g 放入燒杯中,依次加入適量蒸餾水和雙氧水,在電爐上加熱并不斷攪拌,直至無氣泡產(chǎn)生,以去除土樣中的有機(jī)質(zhì),然后用六偏磷酸鈉浸泡24 h 使土粒分散。處理后的土樣用LS-POP(6)激光粒度儀測(cè)定土壤粒徑,所測(cè)粒徑為0.002 ~0.5 mm,并借助儀器附帶軟件中的分級(jí)功能,得到任意2 粒徑之間的體積分?jǐn)?shù)(%)[7],并計(jì)算分形維數(shù)[8-11]。按美國(guó)土壤質(zhì)地分類系統(tǒng)的土壤粒徑分級(jí)分為<0.002、0.002 ~0.05、0.05 ~0.1、0.1 ~0.25、0.25 ~0.5、0.5 ~1 和1 ~2 mm 共7 個(gè)級(jí)別[12]。土壤pH 值用PHS-3C 型酸度計(jì)測(cè)定;速效鉀用火焰光度法;全鉀采用氫氧化鈉熔融法測(cè)定;全氮用半微量開氏法消煮[13]7-59,全自動(dòng)定氮儀測(cè)定;堿解氮用1 mol/L KCl 浸提-流動(dòng)分析儀測(cè)定[13]394-400;土壤含水量用烘箱法測(cè)定。
不同土地利用類型、同種土地利用類型下不同層次的土壤顆粒組成和化學(xué)性質(zhì)不同,導(dǎo)致土壤顆粒的分形維數(shù)也不同。按照土壤顆粒粒徑的不同,分別利用不同模型進(jìn)行擬合分析,得出不同粒徑的分形維數(shù)。
由土壤顆粒體積分形維數(shù)的物理意義可以判斷,土壤顆粒的組成狀況與土壤顆粒體積分形維數(shù)D 值有密切關(guān)系。根據(jù)美國(guó)標(biāo)準(zhǔn)[12],劃分實(shí)驗(yàn)數(shù)據(jù)中不同樣區(qū)和土層黏粒、粉粒、極細(xì)砂粒、細(xì)砂粒、中砂粒、粗砂粒及極粗砂粒的平均累積體積分?jǐn)?shù)(%)。表1 為不同土層深度范圍內(nèi)3 塊樣地土壤顆粒機(jī)械組成分析結(jié)果,可以看出,在土層深度為0 ~30 cm 內(nèi),黏粒(<0.001 mm)質(zhì)量分?jǐn)?shù)在垂直剖面上,隨土層深度的增加,草地不斷增加,耕地逐漸減小,人工林地減小。土壤的物理性黏粒(<0.05 mm)隨土層的變化則是草地遞增,耕地遞減,人工林地遞增。草地、人工林地土壤的物理性黏粒和黏粒的質(zhì)量分?jǐn)?shù)在遞增,而耕地卻在不斷減小。表明土壤覆蓋度的不同,必然導(dǎo)致土壤風(fēng)力侵蝕程度的不同,從而影響土壤質(zhì)地、結(jié)構(gòu)以及其他物理性質(zhì)的差異。
表1 不同土地利用類型不同土層的土壤粒徑分布Tab.1 Average percentage contents of soil particles in different layers of different land use types of soil
根據(jù)不同樣地不同土層的各級(jí)土壤粒徑做D值線性回歸方程,通過公式計(jì)算土地復(fù)墾區(qū)人工林、耕地和草地區(qū)土樣的體積分形維數(shù)D,結(jié)果見表2。
表2 不同土地利用類型不同土層的土壤分形維數(shù)DTab.2 Soil fractal dimension D of land use types and different soil layers
分形維數(shù)的大小反映了土壤質(zhì)地中黏粒、粉粒和砂粒體積分?jǐn)?shù)的變化,隨黏粒體積分?jǐn)?shù)的增多分形維數(shù)增大,隨砂粒體積分?jǐn)?shù)的增多分形維數(shù)減小,同時(shí)土壤顆粒分形維數(shù)與所對(duì)應(yīng)的水分特征曲線的擬合分形維數(shù)呈現(xiàn)出良好的一致性[5]。分形維數(shù)均值表現(xiàn)為耕地>草地>榆樹林地,農(nóng)業(yè)上大多數(shù)的土壤平均值為2.645。表層(0 ~20 cm)土壤分形維數(shù)的大小依次為草地>榆樹林地>耕地,深層(20 ~60c m)土壤分形維數(shù)的大小依次為耕地>草地>榆樹林地。
由表2 可以看出,3 個(gè)樣區(qū)的D 值在2.078 6 ~2.717 5 之間,平均D 值為2.424 1,低于李德成等[14]得出的東北黑面砂土土壤分形維數(shù)平均值為2.704 的結(jié)果。樣區(qū)的地面植被狀況較差,黏粒體積分?jǐn)?shù)較少,特別是人工林,大塊巖石較多,故D 值較小。這是由于海州露天礦排土場(chǎng)是大型煤礦排土場(chǎng),體表多巖石、礦體,土地復(fù)墾區(qū)土質(zhì)是由露天礦排出的頁(yè)巖、砂巖構(gòu)成,部分土地得到了風(fēng)化,地下水資源匾乏,常年干旱少雨,風(fēng)蝕作用嚴(yán)重,土壤質(zhì)地較粗,黏粒、粉粒體積分?jǐn)?shù)較低所致。耕地土壤細(xì)砂粒體積分?jǐn)?shù)越多,分形維數(shù)越小;黏粒體積分?jǐn)?shù)越多,分形維數(shù)越大。草地土壤中粗砂粒體積分?jǐn)?shù)越多,分形維數(shù)越小;黏粒體積分?jǐn)?shù)越多,其分形維數(shù)越大??梢?,土壤分形維數(shù)不僅反映了土壤顆粒的大小,而且體現(xiàn)了土壤質(zhì)地均一的程度[15]。分形維數(shù)越高,表明土壤質(zhì)地相對(duì)較好,通透性相對(duì)較差;分形維數(shù)越低,則表明土壤結(jié)構(gòu)越松散,保水保肥性能越差。
由表1、表2 和表3 可知,土壤顆粒體積分形維數(shù)與土壤顆粒組成有著密切的聯(lián)系,二者與土壤養(yǎng)分質(zhì)量分?jǐn)?shù)的相關(guān)程度基本一致,進(jìn)一步說明了土壤分形維數(shù)是土壤的基本屬性之一。對(duì)分形維數(shù)與顆粒組成進(jìn)行回歸分析,回歸方程的決定系數(shù)(R2)在0.038 5 ~0.668 2 之間,粉粒達(dá)到顯著水平(P <0.01),說明擬合效果良好。黏粒和粉粒體積分?jǐn)?shù)與D 值呈正相關(guān),回歸方程為對(duì)數(shù)函數(shù),其中粉粒與D 值回歸方程的決定系數(shù)最大為0.668 2,表明土壤粉粒體積分?jǐn)?shù)對(duì)D 值的影響最大。粉粒體積分?jǐn)?shù)小于3%時(shí),D 值隨粉粒體積分?jǐn)?shù)的增加而急速增加,然后增加趨勢(shì)變緩。極細(xì)砂粒、細(xì)砂粒、中砂粒、粗砂粒體積分?jǐn)?shù)與D 值呈負(fù)相關(guān)關(guān)系,回歸方程可以用一元二次多項(xiàng)式表示。由此得出,0.05 mm 粒徑是D 值增加或減小的臨界值。對(duì)土壤粒徑體積分?jǐn)?shù)與D 值進(jìn)行多元分析表明,決定D 值大小的是<0.002 mm 粒徑的體積分?jǐn)?shù),其他粒徑體積分?jǐn)?shù)通過影響粒徑的體積分?jǐn)?shù)間接影響D 值;所以,不同立地土壤中黏粒、砂粒體積分?jǐn)?shù)的閾值不同,它們與D 值之間的相關(guān)關(guān)系也有所不同,這與本研究的結(jié)果基本一致。本研究中黏粒、粉粒、粗砂粒體積分?jǐn)?shù)與D 值呈正相關(guān),極細(xì)砂粒、細(xì)砂粒、中砂粒、極粗砂粒體積分?jǐn)?shù)與D 值呈負(fù)相關(guān),砂粒體積分?jǐn)?shù)與D 值相關(guān)性很小。原因可能是后者所研究的土壤為工礦區(qū)排土場(chǎng)土壤,黏粒、砂粒體積分?jǐn)?shù)的閾值相對(duì)很小,對(duì)土壤體積分形維數(shù)的貢獻(xiàn)不明顯。
在經(jīng)典統(tǒng)計(jì)方法中,常用變異系數(shù)從整體上來衡量測(cè)定值的變異程度,變異系數(shù)≤0.15 時(shí)為弱變異性,變異系數(shù)在0.15 ~0.35 時(shí)為中等變異性,變異系數(shù)≥0.35 時(shí)為強(qiáng)變異性。表4 中,變異系數(shù)均<0.15,變異程度較弱。
由表4 可以看出,不同樣區(qū)平均D 值大小關(guān)系為:人工林(2.245 7)<草地(2.502 4)< 耕地(2.547 6)。由結(jié)果可知,耕地的土壤分形維數(shù)最大,原因是耕地經(jīng)常通過一系列的人工耕作來影響土壤,運(yùn)用客土、人工整地、人工施肥等方式快速改變了土壤的理化性質(zhì),這樣耕地土壤的水肥條件相對(duì)人工林和草地充足、植被覆蓋較密,使得耕地土壤黏粒和粉粒體積分?jǐn)?shù)增高,故耕地D 值較大。人工林和草地常年保持無干擾的狀態(tài)下,僅靠自然恢復(fù),所以土壤理化性質(zhì)改良比較慢。草地的土壤分型維數(shù)介于耕地和人工林之間:一方面,由于草地土壤風(fēng)化較好,土壤顆粒介于人工林與耕地之間;另一方面,草場(chǎng)上的牧草根系比較發(fā)達(dá),相對(duì)于林地蓄水保土效果較好,所以草場(chǎng)的土壤分形維數(shù)較人工林地的大。
表3 不同土地利用類型土壤顆粒體積分形維數(shù)與土壤性質(zhì)相關(guān)關(guān)系Tab.3 Soil in different land use type particle volume fractal dimension and the correlation between soil properties
通過方差分析,以土壤不同層次為影響因子,在顯著性水平為0.05 的F 分布表中查得F(9,17)=2.49,在顯著性水平為0.01 的F 分布表中查得F(9,17)=3.68;因?yàn)楫?dāng)F <顯著性水平為0.05 的F時(shí),認(rèn)為影響不顯著??芍寥啦煌瑢哟螌?duì)土壤分形維數(shù)影響不顯著。
結(jié)果表明,在排土場(chǎng)土地復(fù)墾區(qū)土壤不同層次對(duì)分形維數(shù)影響不顯著,說明在垂直方向上土壤分形維數(shù)沒有一定的相關(guān)性。黨亞愛等[16]認(rèn)為,同一地區(qū)不同土層土壤顆粒體積分形維數(shù)差異不顯著,與本研究結(jié)論相似。海州露天礦排土場(chǎng)的土質(zhì)是由采礦后排出的表土及巖層沙石等經(jīng)過壓實(shí)工藝形成的,由于時(shí)間過短,排土場(chǎng)表土不同層次之間的土壤顆粒之間沒有一定的規(guī)律,所以土壤不同土層之間的變化與土壤分形維數(shù)不相關(guān)。
土壤化學(xué)性質(zhì)的改變對(duì)于土壤改良有很重要的作用,排土場(chǎng)通過采用不同植被進(jìn)行土壤改良,導(dǎo)致土壤的化學(xué)性質(zhì)存在差異。通過測(cè)定水解N、全N、速效K 的質(zhì)量分?jǐn)?shù)和土壤pH 值,繪制各化學(xué)性質(zhì)的質(zhì)量分?jǐn)?shù)與土壤分形維數(shù)的相關(guān)曲線,分析土壤化學(xué)性質(zhì)與土壤分形維數(shù)的關(guān)系。
將土壤顆粒分形維數(shù)D 值與pH 值,全N、水解N、速效K 質(zhì)量分?jǐn)?shù)進(jìn)行相關(guān)分析,全N 質(zhì)量分?jǐn)?shù)與土壤分形維數(shù)呈多項(xiàng)式曲線關(guān)系:y=-0.018 2x2-0.135 5x+2.640 2,R2=0.699 4;水解N 質(zhì)量分?jǐn)?shù)與土壤分形維數(shù)呈多項(xiàng)式曲線關(guān)系:y=0.270 1x2-0.683 7x+2.543 5,R2=0.548 8;速效K 質(zhì)量分?jǐn)?shù)與土壤分形維數(shù)呈對(duì)數(shù)曲線關(guān)系:y=-0.379 3lnx+4.762 3,R2=0.420 4。表明海州露天礦排土場(chǎng)土壤分形維數(shù)D 值與土壤全N、水解N、速效K 質(zhì)量分?jǐn)?shù)為負(fù)相關(guān),與pH 值不相關(guān)。
各元素質(zhì)量分?jǐn)?shù)與分形維數(shù)的回歸擬合方程在幾種模型擬合中,擬合效果最佳,即R2的值最大,求得全氮|R|=0.84,水解氮|R|=0.74,速效K|R|=0.65。當(dāng)|R|<0.19 為極低相關(guān)、|R|<0.39 為低度相關(guān)、|R|<0.59 為中度相關(guān)、|R|<0.89 為高度相關(guān)、|R|<1.00 為極高相關(guān)。由此可見,全N、水解N 和速效K 質(zhì)量分?jǐn)?shù)與分形維數(shù)均呈高度負(fù)相關(guān)關(guān)系,分形維數(shù)與穩(wěn)滲率相關(guān)系數(shù)的絕對(duì)值大于分形維數(shù)與初滲率相關(guān)系數(shù)的絕對(duì)值,分形維數(shù)的大小對(duì)土壤含氮含鉀的影響程度是不同的,土壤顆粒分形維數(shù)增大,全N、水解N 和速效K 的質(zhì)量分?jǐn)?shù)隨之變小,而分形維數(shù)對(duì)全N 和水解N 質(zhì)量分?jǐn)?shù)的影響要大于對(duì)速效K 質(zhì)量分?jǐn)?shù)的影響。分形維數(shù)D 與全N、水解N 和速效K 質(zhì)量分?jǐn)?shù)均為負(fù)相關(guān),結(jié)合分形維數(shù)D 與土壤質(zhì)地的關(guān)系,隨著質(zhì)地由砂到黏,土壤顆粒的分形維數(shù)增加,全N、水解N和速效K 的質(zhì)量分?jǐn)?shù)減少。27 個(gè)土樣涉及不同的土地利用類型,有林地、草地和耕地,由于經(jīng)營(yíng)管理方式的不同,土壤的水解N、全N、速效K 質(zhì)量分?jǐn)?shù)及pH 值存在明顯的差異。人工林的全N、水解N及速效K 質(zhì)量分?jǐn)?shù)遠(yuǎn)高于耕地和草地,原因是樹木的根系較深,吸收1 m 以下土層的氮元素較多,而牧草與玉米(Zea mays)的根系較短,吸收1 m 以上土層的氮元素較多,人工林和草地的pH 值比耕地低,說明人工林與草地的保肥效果比耕地的低。
表4 不同土地利用類型區(qū)土壤顆粒體積分形維數(shù)及其擬合方程決定系數(shù)Tab.4 Volume fractal dimensions of soil particles in different areas of land use soil particles and the determination coefficients of the fitting equation
1)不同土地利用類型對(duì)土壤分形維數(shù)影響顯著,其分形維數(shù)的變化與當(dāng)?shù)氐耐寥拉h(huán)境狀況有關(guān),而土壤不同層次對(duì)土壤分形維數(shù)影響不顯著,說明在垂直方向上土壤分形維數(shù)沒有一定的相關(guān)性。不同土地利用類型區(qū)土壤分形維數(shù)表現(xiàn)出不同的變化規(guī)律,所以不同土地利用類型對(duì)土壤質(zhì)地的改良效果不同。
2)林地、耕地和草地樣區(qū)的D 值在2.078 6 ~2.717 5 之間,平均D 值為2.424 1。不同樣區(qū)平均D 值大小關(guān)系為人工林(2.245 7)<草地(2.502 4)<耕地(2.547 6),原因是耕地經(jīng)過多年的客土,人工整地,人工施肥,水肥條件相對(duì)充足,黏粒和粉粒體積分?jǐn)?shù)增高,故D 值較大。
3)水解N、全N 質(zhì)量分?jǐn)?shù)與土壤分形維數(shù)D 呈多項(xiàng)式曲線關(guān)系,速效K 質(zhì)量分?jǐn)?shù)與土壤分形維數(shù)呈對(duì)數(shù)曲線關(guān)系,土壤全N、水解N、速效K 質(zhì)量分?jǐn)?shù)與排土場(chǎng)土壤分形維數(shù)D 值為負(fù)相關(guān),與pH 值不存在相關(guān)性。
由于分形理論是一門發(fā)展中的學(xué)科,許多概念都還在不斷發(fā)展與深化之中;因此,分形幾何在土壤科學(xué)研究中的應(yīng)用還有許多問題需要解決。直接運(yùn)用分形數(shù)學(xué)模型從不同的角度對(duì)土壤結(jié)構(gòu)進(jìn)行分形描述雖已有了一定的基礎(chǔ);但還需要不斷完善,比如,新的分形方法的引入、利用分形維數(shù)與其他參數(shù)關(guān)系對(duì)土壤進(jìn)行改良。實(shí)驗(yàn)過程中由于時(shí)間有限和地形、地質(zhì)等因素的限制,取樣可能會(huì)缺少代表性,因此在對(duì)工礦區(qū)排土場(chǎng)進(jìn)行研究和探討時(shí)需要花費(fèi)大量的時(shí)間和精力。工礦區(qū)排土場(chǎng)為人為堆積,很多實(shí)驗(yàn)數(shù)據(jù)可能會(huì)不合常理,但都屬正?,F(xiàn)象。要對(duì)這些真實(shí)數(shù)據(jù)進(jìn)行整理和分析,為保持?jǐn)?shù)據(jù)的真實(shí)性加大研究力度,為工礦區(qū)排土場(chǎng)的治理做出應(yīng)有的貢獻(xiàn)。
在實(shí)驗(yàn)、資料整理和數(shù)據(jù)處理過程中,得到了呂向楠、李葉鑫和班云云等同學(xué)的大力幫助,在此表示衷心感謝。