張 娜，張紅玲，屈忠義，張棟良

（1.內(nèi)蒙古農(nóng)業(yè)大學水利與土木建筑工程學院，內(nèi)蒙古呼和浩特010018；2.寧夏回族自治區(qū)水利科學研究院，寧夏銀川750021）

內(nèi)蒙古河套灌區(qū)不同土層有機質(zhì)空間變異的分形

張 娜1，2，張紅玲2，屈忠義1，張棟良1

（1.內(nèi)蒙古農(nóng)業(yè)大學水利與土木建筑工程學院，內(nèi)蒙古呼和浩特010018；2.寧夏回族自治區(qū)水利科學研究院，寧夏銀川750021）

以內(nèi)蒙古河套灌區(qū)解放閘灌域為研究對象，以區(qū)域內(nèi)間隔4 km等間距112個取樣點共448個樣本為基礎，通過多重分形及聯(lián)合多重分形方法，對有機質(zhì)含量在不同土層（0～20 cm、20～40 cm、40～70 cm與70～100 cm）的空間變異特性及其相互關(guān)系進行研究。結(jié)果顯示，有機質(zhì)含量0～20 cm至70～100 cm土層的多重分形譜寬度依次為 0.933、0.795、0.824、0.805，即在0～20 cm土層空間變異性最大，而20～40 cm土層空間變異性最小。不同土層多重分形譜曲線均表現(xiàn)為左鉤狀（Δf＞0），表明在灌域內(nèi)土壤有機質(zhì)的空間分布中，數(shù)值較大的數(shù)據(jù)占主導地位。同時從土壤表層到埋深 1m處 D（1）／D（0）依次為 0.983、0.980、0.973、0.965，其值逐漸減小，即空間分布逐漸由密集區(qū)域向稀疏區(qū)域轉(zhuǎn)化。聯(lián)合多重分形分析表明0～20 cm土層與20～40 cm、40～70 cm土層空間變異相關(guān)程度較高，而與70～100 cm土層空間變異相關(guān)程度較低。

有機質(zhì)；不同土層；空間異質(zhì)性；多重分形；聯(lián)合多重分形；內(nèi)蒙古河套灌區(qū)

由于受到自然因素和人為因素的共同作用，土壤特性在空間分布上表現(xiàn)出不同程度的變異性［1－2］，且研究區(qū)域和尺度不同時，土壤特性空間分布的復雜性和異質(zhì)性程度并不相同［3－4］。此外，土壤是三維自然空間實體，在水平方向和垂直方向上都存在空間變異性［5］。垂直方向上不同土層土壤特性的空間變異特征及其主要影響因素并不完全相同［6］，且垂直方向的不同土層土壤特性空間變異性之間存在一定的相互關(guān)系。這種相互關(guān)系的研究有助于呈現(xiàn)土壤特性在三維空間上的變異特征，使土壤過程的模擬更貼近真實情況，也有助于實現(xiàn)通過某一土層土壤特性的空間變異特征反映其它土層土壤特性的空間變異特征。因此，定量分析不同土層土壤特性的空間變異特征及其相互關(guān)系已成為目前水土科學的研究熱點之一。

近年來，各國學者對不同尺度土壤特性參數(shù)的空間變異進行了探討，從不同角度提出了一些值得借鑒的研究方法，如地質(zhì)統(tǒng)計學、方差分析、自相關(guān)分析、譜分析、小波分析和分形理論等［7］。分形理論作為一種能夠解釋研究對象復雜結(jié)構(gòu)并定量表征均一程度的工具，已被廣泛應用于土壤特性空間分布研究中［8－11］。多重分形方法將分形集聯(lián)合起來形成多重分形譜用以描述土壤特性的空間變異特性。該方法不需要數(shù)據(jù)符合某一種分布，且其結(jié)果不受研究對象取樣尺度的影響，是以研究對象尺度為基礎劃分多尺度特征的結(jié)果。同時聯(lián)合多重分形理論可以定量分析和確定同一幾何支撐下不同研究對象在多尺度上的相互關(guān)系［12－13］。

有機質(zhì)含量是反映土壤地力的重要指標，其空間變異性的研究具有重要意義。然而，目前應用分形理論進行土壤有機質(zhì)的空間分布特征研究，主要集中在田間尺度及某一土層上，而大區(qū)域尺度上不同土層空間變異性的相互關(guān)系研究較少。為此，本文以內(nèi)蒙古河套灌區(qū)解放閘灌域內(nèi)間隔4 km等間距112個取樣點共448個樣本為基礎，通過多重分形及聯(lián)合多重分形方法，研究土壤有機質(zhì)不同土層空間變異及其相互關(guān)系，對于全面理解土壤有機質(zhì)空間分布具有重要意義，并為農(nóng)業(yè)技術(shù)研究中野外采樣系統(tǒng)設計提供一定的參考。

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

1.1 研究區(qū)概況

內(nèi)蒙古河套灌區(qū)位于黃河流域上游，是我國第三大灌區(qū)且為亞洲最大的灌區(qū)，總土地面積1.12×104km2，其中灌溉地面積達到5.70×103km2。本文以河套灌區(qū)解放閘灌域作為研究對象，解放閘灌域是河套灌區(qū)第二大灌域，地處東經(jīng) 106°43′～107°27′，北緯 40°34′～41°14′，東南緊鄰黃河，西北地處狼山腳下。西北以總排干溝為界，東南則以總干渠為界，大部分位于內(nèi)蒙古自治區(qū)杭錦后旗所轄境內(nèi)。解放閘灌域總控制面積2.1568×105hm2，灌溉面積1.421×105hm2，其中耕地 1.2193×105hm2，林果地面積 5.507×103hm2，牧草地面積 1.461×104hm2。灌域南北長約87 km，東西寬約78 km，呈三角形狀，解放閘灌域?qū)僦袦貛Ц咴蜌夂?，風大雨少，氣候干燥。地形較為平緩，地勢由東南向西北微度傾斜，坡度約為0.02%。海拔在1 030～1 046m之間。

1.2 樣點布設及土樣采集

于2011年8月15日至9月10日期間，根據(jù)解放閘灌域?qū)嶋H控制面積，通過ArcGIS在數(shù)字化圖（1∶12.5萬）上按照4 km×4 km進行網(wǎng)格設計，按0～20、20～40、40～70 cm與 70～100 cm共4層進行取樣。取樣區(qū)域北部以總排干溝為界，東部以四排干溝為界。采樣點位置均用GPS定位，遍布整個灌域，基本屬于平均分布。經(jīng)統(tǒng)計，共取有效點112個。取樣點位置如圖1所示。

1.3 測試項目及方法

按照土壤學實驗規(guī)范，土樣經(jīng)過風干，研磨，過篩（2mm）后，參考《土壤、水、植物理化分析教程》［14］中的濃硫酸－重鉻酸鉀容量法測定。

1.4 研究方法

多重分形方法主要用于分析單一變量的空間變異性，描述多重分形的參量有兩種，一種為 D（q）和q，另一種為α（q）和 f（q），計算過程可參考文獻［12－13］。

式中，D（q）為廣義維數(shù)；q為質(zhì)量概率；pi（δ）為統(tǒng)計矩；δ為尺度；pi（δ）α（q）為研究對象μ的奇異指數(shù)；f（q）為研究對象奇異指數(shù)為奇異指數(shù)的維數(shù)分布函數(shù)；奇異指數(shù)α（q）和奇異指數(shù)維數(shù)分布函數(shù)f（q）的關(guān)系曲線稱為多重分形譜。

其中研究對象的多重分形譜的寬度越大，表示研究對象的空間變異性越強；若研究對象的多重分形譜左偏趨勢比較明顯，表示研究對象的高值信息對研究對象的空間分布特征影響顯著，若研究對象的多重分形譜右偏趨勢比較明顯，表示研究對象的低值信息對研究對象的空間分布特征影響顯著。

將多重分形拓展為聯(lián)合多重分形后，可利用聯(lián)合多重分形方法確定不同變量在多尺度上的相關(guān)性，計算公式可參考文獻［15－17］。

圖1 取樣點分布示意圖Fig.1 Distributionmap of sampling sites

聯(lián)合 多 重 分 形 參 數(shù) α1（q1，q2）、α2（q1，q2）、f（α1，α2）的三維圖稱為研究對象1和研究對象2的聯(lián)合多重分形譜。當確定了聯(lián)合多重分形譜參數(shù)α1（q1，q2）、α2（q1，q2）、f（α1，α2），可以通過以下兩種方法分析研究對象1和研究對象2之間的相互關(guān)系：一種方法是將聯(lián)合多重分形譜三維圖轉(zhuǎn)化為二維圖，用灰度值表示 f（α1，α2），然后將 f（α1，α2）投影到α1（q1，q2）、α2（q1，q2）的平面上，得到研究對象1和研究對象2聯(lián)合多重分形譜的灰度圖，如果研究對象1和研究對象2聯(lián)合多重分形譜的灰度圖相對集中且沿對角線方向延伸，則可定性的說明研究對象1和研究對象2具有某些相同的分布或者相關(guān)程度較高，如果研究對象1和研究對象2聯(lián)合多重分形譜的灰度圖分布比較離散，則反之；另一種方法是計算α1（q1，q2）、α2（q1，q2）的相關(guān)系數(shù)，如果α1（q1，q2）、α2（q1，q2）的相關(guān)性顯著，則可定量的說明在多尺度上研究對象1與研究對象2的相關(guān)程度高，反之，可定量的說明研究對象1與研究對象2的的相關(guān)程度低。

2 結(jié)果與分析

2.1 有機質(zhì)含量的統(tǒng)計特征值

不同土層土壤有機質(zhì)的統(tǒng)計特征結(jié)果（表1）顯示，土壤有機質(zhì)含量的均值隨著土層深度的增加而減少，即從土壤表層到埋深1米處土壤有機質(zhì)含量遞減。這一規(guī)律產(chǎn)生的原因為：對于原始土壤來說，微生物是土壤有機質(zhì)的最早來源，而隨著生物的進化和成土過程的發(fā)展，動植物殘體就成為土壤有機質(zhì)的基本來源。在通常的自然植被條件下，土壤有機質(zhì)絕大部分直接來源于土壤上生長的植物殘體和根系分泌物，故土壤表層的有機質(zhì)含量較高，且土壤表層晝夜溫差較大，微生物分解活動較弱，土壤有機質(zhì)更易于積累。

表1 土壤有機質(zhì)含量描述性統(tǒng)計特征值Table 1 Descriptive statistics of soil organicmatter concentration

變異系數(shù)的大小反映總體內(nèi)部各樣本之間的變異程度，一般可將樣本的變異程度分為3級：CV＜10%為弱變異性，10%≤CV≤100%為中等變異性，CV＞100%為強變異性［13］。由表1可知，不同土層有機質(zhì)含量的變異系數(shù)介于37.52%～60.40%，均屬中等變異性。變異系數(shù)隨著土層深度的增加而增加。

2.2 基于多重分形配分函數(shù)的分形特征判定

根據(jù)多重分形配分函數(shù)的性質(zhì)［17］，當 log（ε）～log（χq（ε））關(guān)系為線性時，表明研究區(qū)域內(nèi)變量分布具有多重分形特征［18］。擬合不同土層下土壤有機質(zhì) log（ε）與 log（χq（ε））的關(guān)系，圖 2給出了雙對數(shù)坐標下不同 q值對應的配分函數(shù)值。由圖2可知，在 －10≤q≤10時，log（ε）～log（χq（ε））之間成一定的線性關(guān)系，表明研究區(qū)域內(nèi)土壤有機質(zhì)含量具有多重分形特征。表2給出了不同土層的有機質(zhì)含量 log（ε）－log（χq（ε））擬合的斜率 K和決定性系數(shù)R2，可知 log（ε）－log（χq（ε））之間均有良好的線性關(guān)系，決定性系數(shù) R2位于0.81～1.00之間，在整個空間尺度上表現(xiàn)出良好的標度不變性，同時，當 q取不同值時，log（ε）與 log（χq（ε））之間擬合直線的斜率均不同，這表明了有機質(zhì)含量在不同土層的空間分布上具有多重分形特征。

2.3 基于廣義維數(shù)譜的多重分形分析

根據(jù)多重分形廣義維數(shù)譜算法對不同土層有機質(zhì)含量進行多重分形分析，在－10≤q≤10的變化范圍內(nèi)得到土壤有機質(zhì)分布的廣義維數(shù)譜 D（q），有機質(zhì)不同土層廣義維數(shù)譜曲線 q－D（q）見圖3。

根據(jù)廣義分形維數(shù)隨質(zhì)量概率統(tǒng)計矩階數(shù)變化的趨勢對有機質(zhì)含量在研究區(qū)域上的多重分形特征進行進一步判定。當質(zhì)量概率統(tǒng)計矩階數(shù)大于0時，如果隨質(zhì)量概率統(tǒng)計矩階數(shù)的增加，廣義分形維數(shù)依次減小，則可判定研究對象的多重分形特征比較明顯［19－20］。不同土層的廣義分形維數(shù)均隨質(zhì)量概率統(tǒng)計矩階數(shù)的增大而減小，這進一步說明有機質(zhì)在研究區(qū)域內(nèi)的空間分布格局具有多重分形特征，可以利用多重分形方法研究其空間異質(zhì)性。αmax－αmin的差異來獲取研究對象的空間分布及其不規(guī)則特征是多重分形分析的一個重要方面。Δα越大表示分布越不均勻［21］，多重分形譜α－f（α）曲線呈左鉤狀時，即Δf＞0，說明數(shù)據(jù)中較大值的數(shù)量占主導地位，反之，較小數(shù)值的數(shù)量占主導地位［21］。

表 2 log（ε）—log（χq（ε））擬合的斜率 K和決定性系數(shù) R2Table 2 The slope K and determination coefficients R2 of the fitting line between log（ε）and log（χq（ε））

圖2 不同土層土壤有機質(zhì)的log（ε）與log（χ（qε））關(guān)系Fig.2 The relationship of log（ε）and log（χ（qε））of soil organicmatter in different layers

圖3 不同土層有機質(zhì)含量的廣義分形維數(shù)Fig.3 Generalizedmultifractal dimension of soil organicmatter concentration at different soil layer

圖4給出土壤有機質(zhì)在10≥q≥－10情況下不同土層的多重分形譜α－f（α）曲線，表3給出土壤有機質(zhì)含量的多重分形譜參數(shù)。

定量描述土壤有機質(zhì)分布的離散程度可以通過信息熵維數(shù)和盒維數(shù)的比值 D（1）／D（0）反映。D（1）／D（0）越接近于 1時表明土壤特征分布主要集中在密集區(qū)，接近于0時說明分布集中于稀疏區(qū)域［1］?？芍獜耐寥辣韺拥铰裆?1 m處 D（1）／D（0）逐漸較小，認為有機質(zhì)的空間分布逐漸由密集區(qū)域向稀疏區(qū)域轉(zhuǎn)化。

2.4 土壤有機質(zhì)空間變異性的多重分形分析

根據(jù)多重分形譜形狀和多重分形譜譜寬Δα＝

圖4 不同土層有機質(zhì)含量的多重分形譜曲線Fig.4 Multifractal spectrums of soil organic matter concentration at different soil layer

表3 有機質(zhì)含量空間分布多重分形譜參數(shù)Table 3 Multifractal parameters of soil organicmatter concentration

由圖4可知土壤有機質(zhì)的多重分形譜均為“鐘形”或類似開口向下的二次拋物線，但不同土層的多重分形譜的“頂點”、開口大小存在差異，說明在不同土層其空間分布及其不規(guī)則性是有差別的。由圖4及表3可知不同土層間多重分形譜的譜寬Δα具有一定差異：0～20、20～40、40～70 cm及 70～100 cm土層的多重分形譜寬度依次為 0.933、0.795、0.824、0.805，可見有機質(zhì)含量不同土層的空間變異性由大到小依次為：0～20、40～70、70～100 cm及 20～40 cm土層。

由表3可知，不同土層的Δf均大于0，其在多重分形譜曲線上均表現(xiàn)為左鉤狀，表明在內(nèi)蒙古河套灌區(qū)解放閘灌域土壤有機質(zhì)的空間分布中，數(shù)值較大的數(shù)據(jù)占主導地位，其概率分布較大。

2.5 不同土層的聯(lián)合多重分形

為了探討有機質(zhì)含量在不同土層空間變異性的相互關(guān)系，同時便于用淺層易于取樣測定的數(shù)據(jù)，以減少深層取樣測試工作量。繪制了0～20 cm土層與20～40、40～70 cm及70～100 cm土層的聯(lián)合多重分 形 譜，分 別 用 αOM0～20cm、αOM20～40cm及 f（αOM0～20cm、αOM20～40cm）表示α1（q1，q2）、α2（q1，q2）、f（α1，α2），其它土層依次類推。為了便于分析將其聯(lián)合多重分形譜投影到平面上得到聯(lián)合多重分形譜的灰度圖（圖6）。從圖6可知，0～20 cm土層與20～40、40～70 cm的聯(lián)合分形譜灰度圖相對集中且沿對角線方向延伸，表明有機質(zhì)含量在0～20 cm土層的空間分布與20～40 cm及40～70 cm土層空間變異的相似性較高。因此，0～20 cm土層有機質(zhì)含量的空間變異性可較好反映其在20～40 cm及40～70 cm土層的空間變異性狀況；而有機質(zhì)含量在0～20 cm土層與70～100 cm土層的聯(lián)合多重分形灰度圖分布較為分散，說明二者的相關(guān)性很小，因此其空間變異不能相互表達。

為進一步定量化分析土壤有機質(zhì)含量0～20 cm土層與20～40、40～70 cm及70～100 cm土層空間變異性的相互關(guān)系，求解了其對應奇異指數(shù)之間的相關(guān)性。0～20 cm土層與20～40、40～70 cm及70～100 cm土層對應奇異指數(shù)的相關(guān)系數(shù)分別為0.654、0.362、0.032，且 0～20 cm土層與 20～40、40～70 cm土層在0.01水平上顯著，而與70～100 cm土層不顯著。進一步說明有機質(zhì)含量空間變異性在0～20 cm土層與20～40 cm及40～70 cm土層相關(guān)程度較高。

通過一元函數(shù)建立0～20 cm土層分別與20～40、40～70 cm及70～100 cm土層的轉(zhuǎn)換模型，見表4。

圖5 0～20 cm土層有機質(zhì)含量與其它土層的聯(lián)合多重分形譜Fig.5 Jointmultifractal spectrums between organicmatter concentration in 0～20 cm soil layer and other layers inmid-scale

可知：0～20 cm僅與20～40 cm土層的預測效果較好，而與40～70 cm及70～100 cm土層的回歸效果較差。

圖6 0～20 cm土層有機質(zhì)含量與其它土層的聯(lián)合多重分形譜的灰度圖Fig.6 Gray scale image of jointmultifractal spectrums between organicmatter concentration in 0～20 cm soil layer and other layers inmid-scale

表4 0～20 cm土層有機質(zhì)與其它土層轉(zhuǎn)化模型Table 4 The transformationmodel of organicmatter in 0～20 cm and other soil layers

這一結(jié)果與聯(lián)合多重分形分析的結(jié)果相差較大，首先聯(lián)合多重分形研究兩個具有多重分形特征的變量在幾何支撐上的關(guān)系，其涵蓋了研究對象的非均勻性和各向異性等信息，體現(xiàn)的是變量空間分布的信息，而一元函數(shù)僅僅表達了一個變量隨另一個變量的線性變化趨勢，并不能很好地體現(xiàn)變量空間分布特征。可見，一元函數(shù)雖可簡便地建立不同土層變量之間的轉(zhuǎn)換函數(shù)，但并不能很好地體現(xiàn)變量空間變化的信息。基于聯(lián)合多重分形結(jié)果建立兩個變量考慮空間變異情況的轉(zhuǎn)換函數(shù)有待于進一步深入研究。

3 結(jié) 論

通過多重分形及聯(lián)合多重分形方法分析了內(nèi)蒙古河套灌區(qū)解放閘灌域有機質(zhì)含量在不同土層的空間變異特性及其相關(guān)性，得到如下結(jié)論：

1）有機質(zhì)含量在不同土層均為中等變異性，且變異系數(shù)隨著土層深度的增加而增加。

2）內(nèi)蒙古河套灌區(qū)解放閘灌域有機質(zhì)含量在空間分布上具有多重分形特征，其在不同土層的空間變異性由大到小依次為：0～20、40～70、70～100 cm及20～40 cm土層。同時在有機質(zhì)的空間分布中，數(shù)值較大的數(shù)據(jù)占主導地位。

3）聯(lián)合多重分形表明0～20 cm土層與20～40、40～70 cm土層空間變異相關(guān)程度較高，而與70～100 cm土層空間變異的相關(guān)程度較低。

4）通過一元函數(shù)建立0～20 cm土層有機質(zhì)含量分別與20～40、40～70 cm及70～100 cm土層的轉(zhuǎn)換模型，其回歸效果0～20 cm土層僅與20～40 cm土層較好，而與40～70 cm及70～100 cm土層較差。

5）基于聯(lián)合多重分形結(jié)果建立兩個變量考慮空間變異情況的轉(zhuǎn)換函數(shù)有待于進一步深入研究。

［1］ 管孝艷，楊培嶺，任樹梅，等.基于多重分形理論的壤土粒徑分布非均勻性分析［J］.應用基礎與工程科學學報，2009，17（2）：196-205.

［2］ Jacinthe PA，LaiR.Spatial variability of soil propertiesand trace gas fluxes in reclaimed mine land of southeastern Ohio［J］.Geoderma，2006，136（3－4）：598-608.

［3］ Eghball B，Schepers JS，Negahban M，et al.Spatial and temporal variability of soil nitrate and corn yield：multifractalanalysis［J］.A-gronomy Journal，2003，95（2）：339-346.

［4］ Caniego F J，Espejo R，Martin M A，etal.Multifractal scaling of soil spatial variability［J］.Ecological Modelling，2005，182：291-303.

［5］ 邵明安，王全九，黃明斌.土壤物理學［M］.北京：高等教育出版社，2006.

［6］ 劉繼龍，馬孝義，付 強，等.不同土層土壤特性空間變異性關(guān)系的聯(lián)合多重分形研究［J］.農(nóng)業(yè)機械學報，2012，43（5）：37-42.

［7］ 王 軍，邱 揚.土地質(zhì)量的空間變異與尺度效應研究進展［J］.地理科學進展，2005，24（4）：28-35.

［8］ Kravchenko A N，Bullock DG，BoastCW.Jointmultifractal analysis of soil yield and terrain slope［J］.Agronomy Journal，2000，92：1279-1290.

［9］ Lovejoy S，Schertzer D.Multifractal，cloud radiances and rain［J］.Journal of Hydrology，2006，322：59-88.

［10］ Cheng Q.Multifractal interpolation［C］／／Proceedings of the Fifth Annual Conference of the International Association for Mathematical Geology.Trondheim，Norway：［s.n.］，1999：245-250.

［11］ Li Y，LiM，Robert H.Single and jointmultifractal analysis of soil particle size distributions［J］.Pedosphere，2011，21（1）：75-83.

［12］ Zeleke T B，Si B C.Characterizing scale－dependent spatial relationships between soil properties using multifractal techniques［J］.Geoderma，2006，134（3－4）：440-452.

［13］ Zeleke TB，SiBC.Scaling relationshipsbetween saturated hydraulic conductivity and soil physical properties［J］.Soil Sci Soc Am J，2005，69（6）：1691-1702.

［14］ 張 韞.土壤、水、植物理化分析教程［M］.北京：中國林業(yè)出版社，2011.

［15］ Perfect E，Centry RW，Sukop M C，et al.Multifractal Sierpinski carpets：theory and application to upscaling effective saturated hydraulic conductivity［J］.Geoderma，2006，134：240-252.

［16］ Scheuring I，RiediRH.Application ofmultifractal to the analysisof vegetation pattern［J］.Journal of Vegetable Science，1994，5（4）：489-496.

［17］ Wang D，F(xiàn)u B J，ZhaoW W，et al.Multifractal characteristics of soil particle size distribution under different land－use types on the Loess Plateau，China［J］.Catena，2008，72（1）：29-36.

［18］ Lee C K.Multifractal characteristics in air pollutant concentration time series［J］.The Air Soil Pollut，2002，135（1－4）：389-409.

［19］ Eghball B，Schepers JS，Negahban M，et al.Spatial andtemporal variability of soil nitrate and corn yield：multifractal analysis［J］.A-gron J，2003，95（2）：339-346.

［20］ Peitgen H O，Jurgens H，Saupe D.Chaos and Fractals：New Frontier of Science［M］.NewYork：Springer-Verlag，1992：737.

［21］ Kravchenko A N，CharlesW B，Donald GB.Multifractal analysisof soilspatial variability［J］.Agronomy Journal，1999，91（6）：1033-1041.

［22］ Zeleke T B，Si B C.Scaling properties of to-pographic indices and crop yield：multifractal and jointmultifractal approaches［J］.Agronomy Journal，2004，96（4）：1082-1090.

Fractal study on organicmatter spatial heterogeneity of different soil layers in Inner Mongolia Hetao Irrigation District

ZHANG Na1，2，ZHANG Hong-ling2，QU Zhong-yi1，ZHANG Dong-liang1
（1.College ofWɑter Resourcesɑnd Civil Engineering，Inner MongoliɑAgriculturɑl University，Hohhot，Inner Mongoliɑ010018，Chinɑ；2.Scientific Reseɑrch Institute ofWɑter Conservɑncy of Ningxiɑ，Yinchuɑn，Ningxiɑ750021，Chinɑ）

The Jiefangzha irrigation district in the Inner Mongolia Hetao irrigation area was taken as research subject.Based on 448 samples of112 sampling points at4 km of interval distance in the region，and usingmulti-fractal and jointmulti-fractalmethod，we studied the spatial variability of organicmatter concentration in different soil layers（0～20 cm，20～40 cm，40～70 cm，70～100 cm）and their characteristics in themultiscale correlation.Results showed that the organicmatter concentration in 20～40 cm soil layer had the smallest spatial variability，and in the surface was the strongest（Δαof 0～20 cm to 70～100 cm soil layerwas 0.933、0.795、0.824 and 0.805）.Themulti-fractal spectrum curve of different layers showed a left hook（Δf＞0），indicating that the larger data dominated.Also，the spatial distribution gradually transformed from the dense region to the sparse region from the surface to the 1m depth（D（1）／D（0）from the surface to the 1m depth was0.983、0.980、0.973 and 0.965）.The jointmulti-fractal contourmap and the corresponding singular index of different soil layers showed that the correlation between 0～20 cm and 20～40 cm、40～70 cm were relatively high，butwith 70～100 cm soil layer was relatively low.

organicmatter；different soil layers；spatial variation；multifractal；jointmultifractal；Inner Mongolia hetao irrigation district

S154.4

A

1000-7601（2017）05-0157-07

10.7606／j.issn.1000-7601.2017.05.23

2016-08-04

2016-10-27

國家自然科學基金項目（51469020）

張 娜（1988—），內(nèi)蒙古呼和浩特人，博士，研究方向為非充分灌溉理論。E-mail：710840895＠qq.com。

屈忠義（1969—），內(nèi)蒙古巴彥淖爾人，博士，教授，主要從事節(jié)水灌溉理論與區(qū)域水鹽監(jiān)測方面的科學研究。E-mail：