生态学报  2015, Vol. 35 Issue (14): 4625-4632

文章信息

赵明月, 赵文武, 刘源鑫
ZHAO Mingyue, ZHAO Wenwu, LIU Yuanxin
不同尺度下土壤粒径分布特征及其影响因子--以黄土丘陵沟壑区为例
Comparative analysis of soil particle size distribution and its influence factors in different scales: a case study in the Loess Hilly-gully area
生态学报, 2015, 35(14): 4625-4632
Acta Ecologica Sinica, 2015, 35(14): 4625-4632
http://dx.doi.org/10.5846/stxb201311272828

文章历史

收稿日期:2013-11-27
修订日期:2014-09-09
不同尺度下土壤粒径分布特征及其影响因子--以黄土丘陵沟壑区为例
赵明月1, 2, 赵文武1, 2 , 刘源鑫1, 2    
1. 地表过程与资源生态国家重点实验室, 北京师范大学, 北京 100875;
2. 北京师范大学 资源学院, 北京 100875
摘要:土地利用方式和环境因子是影响土壤粒径的重要因素,尺度不同其影响效应差异明显。研究不同尺度表层土壤粒径与土地利用、环境因子的关系对评价黄土丘陵沟壑区土壤物理性质具有积极意义。选择黄土丘陵沟壑区安塞集水区和沐浴小流域作为研究区,探讨两个尺度上表层土壤粒径含量、分布特征及其与土地利用类型和环境因子的关系。研究结果表明:(1)研究区域内表层土壤颗粒主要为砂粒和粉粒,在小流域和集水区尺度上,各粒径百分含量属于中等变异,D值为弱变异,但随着研究区由沐浴小流域变化到安塞集水区,粒径和D值的空间变异性均有所提高;(2)尺度不同,土地利用对土壤粒径的影响效应不同,在沐浴小流域不同土地利用类型的砂粒含量从高到低依次为荒草地>农用地>林地>灌木林地>园地,在集水区的变化顺序依次为荒草地>灌木林地>林地>农用地>园地,粉粒含量的次序均与砂粒相反,小流域土壤粒径分形维数D依次为灌木林地>荒草地>林地>园地>农用地,在安塞集水区为灌木林地>荒草地>农用地>林地>园地;(3)在沐浴小流域和安塞集水区,土壤颗粒分形维数与土壤砂粒、黏粒百分含量呈极显著正相关关系,与土壤粉粒百分含量呈极显著负相关关系,且砂质壤土D值高于壤土和粉壤土;(4)在沐浴小流域,植被盖度、相对海拔和坡位相对于土地利用类型、坡向和坡度对土壤粒径影响较大;而在安塞集水区植被盖度、土地利用类型和坡度对土壤粒径影响较大。
关键词土地利用    环境因子    土壤粒径    分形维数    尺度    
Comparative analysis of soil particle size distribution and its influence factors in different scales: a case study in the Loess Hilly-gully area
ZHAO Mingyue1, 2, ZHAO Wenwu1, 2 , LIU Yuanxin1, 2    
1. State Key Laboratory of Earth Surface Processes and Resource Ecology, Beijing Normal University, Beijing 100875, China;
2. College of Resources Science and Technology, Beijing Normal University, Beijing 100875, China
Abstract:The land use types and environmental factors have a great significance on soil particle size distribution (PSD). However the effects are different in different scales. Research on the relationship between topsoil particle size distribution and land use & environmental factors can provide a basis for assessing topsoil physical properties in loess hilly-gully area. The study was conducted in the Muyu small watershed and the Ansai catchment in the loess hilly area of the Loess Plateau in China, and the comparative analysis of soil particle size distribution and its influence factors in different scales were carried out. The research results showed that: (1) The topsoil in the loess hilly-gully areas of the Loess Plateau are mainly composed of sand and silt in both small watershed and catchment. The percentage of each particle size is medium spatial variability, and the spatial variability of D is far below of the particle size. But it appears clearly that the spatial variabilities are scale dependent and it increases with the increase of study area. (2)The effects of land use types on soil particle size are different at two scales. The sand content varied with land use in the order: grassland> farmland> woodland> shrub land> orchard land in the Muyu small watershed, and grassland> shrub land> woodland> orchard land> farmland in the Ansai catchment. The silt content varied with land use in the opposite order of the sand content. The D value varied with land use in the order: shrub land> grassland>woodland> orchard land> farmland in the Muyu small watershed, and shrub land>grassland>farmland>woodland>orchard land in the Ansai catchment. (3) In both Muyu small watershed and the Ansai catchment, D is very significant positive correlation with the soil sand content and clay content, significant negative correlation with soil silt content, and D value is much higher in the sandy loam than that in the loam and silt loam. (4) Vegetation cover, altitude and position have more influence on the soil particle size than other factors in the Muyu small watershed, while vegetation cover, land use, and slope influence the soil particle much more greatly in the Ansai catchment.
Key words: land use    environmental factors    soil particle size    fractal dimension    scale    

尺度是指在研究某一物体或现象时所采用的空间或时间维,往往以粒度和幅度来表达。尺度效应不仅依赖于事物的本身,而且与整个环境或背景有关[1, 2, 3, 4]。土壤由不同粒级土粒组成,各粒级的含量差异很大。土壤粒径分布PSD)影响土壤的水力特性、土壤肥力状况以及土壤侵蚀等,是土壤重要的物理特性之一[5, 6, 7, 8]。由于成土母质、气候、地形、植被和人类干扰程度不同,不同尺度的土壤性质存在空间变异性。在较大尺度上,土壤主要受母质、气候、地形等因子的控制;在较小尺度上,主要受土壤特性、植被覆盖、微立地以及干扰的影响[9, 10]。其中,土地利用方式和环境因子是影响土壤粒径的重要因素[11, 12, 13]。目前,国内外土壤粒径空间变异性的研究已有很多文献[14, 15, 16, 17, 18],但是在不同空间尺度上讨论土壤粒径的空间分布还较为少见。从不同尺度上探讨土壤粒径的含量和空间变异性,分析表层土壤粒径含量及其与环境因子的关系,有利于从不同尺度系统把握土壤粒径的空间变异特征及影响因素。

黄土高原地区地形破碎,土地利用类型复杂多样,土壤侵蚀程度差异较大。对于该地区土壤物理性质的研究较多,但主要集中于特定区域不同土地利用类型下土壤的物理性质差异研究,而忽略不同尺度间的对比分析。本研究选择黄土丘陵沟壑区安塞集水区和集水区内的沐浴小流域,从两个尺度上分析土壤粒径与环境因子的关系,以期为黄土丘陵沟壑区土壤侵蚀的防治和土地利用结构调整提供理论参考。

1 材料与方法 1.1 研究区概况

安塞集水区和沐浴小流域地处黄土高原腹地(108°5′44″—109°26′18″E,36°30′45″—37°19′3″N),地形地貌复杂多样,梁峁起伏、沟壑纵横,属典型的黄土丘陵沟壑区。该区域属中温带大陆性半干旱季风气候,年平均气温8.8 ℃,年平均降水量505.3 mm,其中74%集中在6—9月。集水区面积1334.00 km2,属于强度水土流失区,也是西北典型的生态环境脆弱区。土壤以黄绵土为主,约占总面积的95%。土壤成土母质主要有黄土和洪积冲积两大类,黄土土层深厚,质地组成以粉沙为主,抗蚀能力低。沐浴小流域位于安塞县郝家坪乡,面积1.30 km2,是安塞集水区的典型小流域。该小流域地处安塞集水区中游,土地利用方式以农耕地、林地、荒草地、灌木林地为主,农耕地面积较小且分布零散,总体人为干扰较少。

1.2 土壤采样

根据土地利用、地形、土壤特征等因素,综合考虑流域水系分布与道路可达性,经实地踏勘于2012年7、8月在安塞集水区和沐浴小流域分别选取了79、35个采样点。每个样点内设置5个重复随机采样。土壤粒径采用比重计法测定。样点分布情况见图 1[19],样点的土地利用类型和坡位分布情况见表 1

图1 研究区位置示意和采样点分布图 Fig.1 Location of study area and the distribution of sampling sites
表1 采样点在土地利用类型及坡位上的分布 Table 1 The distribution of sampling sites under different land use and slop positions
土地利用类型与坡位 Land use types and slope position小流域尺度 Small watershed scale集水区尺度 Watershed scale
灌木林地Shrubland411
林地Woodland1324
园地Orchard15
荒草地Grassland1424
耕地Farmland315
上坡Upslope1022
中上坡Upper slope812
中坡Middle slope712
中下坡Middle-lower slope413
下坡Down slope620
1.3 分形维数计算与统计分析

杨培玲等[20]将Katz的计算方法进行了改进,通过粒径分布与对应的重量分布相联系,提出了用粒径重量分布资料表征的土壤分维模型。本研究根据此模型,计算方法如下:

假定由大于某一粒径Ri (Ri Ri+1,i=1,2,3,…,i)的土粒构成的体积V(r > Ri)可表示为:

式中,Ri是特征尺度,Cvλv是描述颗粒形状、尺度的常数。

Ri=0时,(1)式变为计算全部土壤颗粒总体积VT,即V(r>Ri) = VT = Cv

Ri=Rmax(Rmax为最大粒径)时,λv=Rmax

因此有:

调整(2)式可得:

对上式两边同时取对数,即可得到土壤颗粒组成的分形维数计算公式:

对于某一土壤各级粒径的土粒含量进行整理,并以lg()为纵坐标,lg( )为横坐标作图,得到直线的斜率即为3-D,即可得到土壤分形维数D

采用SPSS 18.0对表层土壤容重和粒径进行描述性统计分析,采用CANOCO 4.5[20, 21]对土壤分形维数、粒径和土地利用及环境因子之间的关系进行典范对应分析分析(CCA)[22]。CCA分析需要两个矩阵,其一是物种数据矩阵,即本研究中的土壤粉粒含量、砂粒含量、黏粒含量和分形维数D,其二是环境因子数据矩阵,即本研究中的土地利用类型和环境因子。

2 结果与分析 2.1 不同尺度土壤粒径和土壤颗粒分形维数特征分析

对土壤粒径和分形维数进行描述性统计(表 2)。由表 2可知,研究区域内表层土壤颗粒以砂粒和粉粒为主,其中沐浴小流域的土壤砂粒、粉粒和黏粒的含量分别为50.27%、38.83%和10.89%;安塞集水区表层土壤砂粒、粉粒和黏粒的含量分别为46.20%、43.06%和10.74%。在小流域和集水区尺度上,各粒径百分含量空间差异性不大,属于中等变异(10%<CV<100%),但随着研究区由小流域到集水区,变异性增加。

表2 土壤粒径和分形维数的描述性统计 Table 2 Descriptive statistics of soil particle size and fractal dimension
统计值 Statistics尺度 Scales样本数 Sample number砂粒 Sand粉粒 Silt黏粒 ClayD
平均值Average沐浴3550.2738.8310.892.14
安塞7946.2043.0610.742.09
最小值Minimum沐浴3528.7611.597.841.86
安塞7922.7510.26.251.71
最大值Maximum沐浴3578.7661.5914.452.55
安塞7982.7563.615.452.55
标准差Standard deviation沐浴3510.510.451.590.15
安塞7912.6212.082.070.16
变异系数Coefficient of variation沐浴3520.8926.9114.357.01
安塞7927.3228.0519.327.66

根据粒径分级及公式(2—4),对粒径进行分形维数的计算,得到土壤PSD分形维数D值。沐浴小流域D值介于1.86—2.55之间,平均值为2.14,变异系数为7.01;安塞集水区D值介于1.71—2.55之间,平均值为2.09,变异系数为7.66。在不同研究尺度D值均属于弱变异性,但随着研究区由沐浴小流域到安塞集水区,各粒径含量和D值的空间变异系数均有增加。

2.2 不同土地利用类型的土壤粒径分布与分形特征

不同土地利用类型的土壤粒径和分形维数见表 3。由表 3可知,不同土地利用类型下的砂粒含量在沐浴小流域含量从高到低依次为:荒草地>农用地>林地>灌木林地>园地,在集水区的次序为:荒草地>灌木林地>林地>农用地>园地;在两个研究区粉粒含量在不同土地利用类型下的次序均与砂粒相反。黏粒含量变化幅度不大。在沐浴小流域,灌木林地的D值最高(2.18)、荒草地其次(2.17)、农用地最低(2.09),D值从高到低依次为灌木林地>荒草地>林地>园地>农用地。在安塞集水区,依然是灌木林地的D值最高(2.11),园地最低(2.06),从高到低依次为灌木林地>荒草地>农用地>林地>园地。

表3 不同尺度土地利用类型的土壤粒径分布和分形维数 Table 3 Soil particle size distribution and fractal dimension of different land use types at different scales
土地利用类型 Land use types尺度 Scales 砂粒含量/% Sand content 粉粒含量/% Silt content 黏粒含量/% Clay contentD
农用地Farmland 安塞 47.0042.8810.132.09
沐浴46.7543.669.582.10
荒草地Grass land安塞47.4141.8710.722.10
沐浴52.5436.5110.952.17
园地Orchard安塞41.1846.4612.362.07
沐浴42.8343.9413.232.17
灌木林地Shrubland 安塞44.2144.1411.652.11
沐浴51.5037.0511.452.18
林地Woodland 安塞46.4543.1710.382.06
沐浴48.8340.3810.792.11

土壤理化性质的变异大部分由土地利用类型变化引起,许明祥等[13]在对黄土丘陵沟壑区土壤理化性质空间变异性研究中,发现土地利用类型可以解释97%的土壤质量指标的变异性,即土地利用类型是黄土丘陵区土壤理化性质变异的主导因子。在本研究中,不同土地利用类型结合植被盖度,在不同尺度上对表层土壤分形维数的影响作用不同。随着研究区由小流域到集水区,不同土地利用类型的分形维数D是降低的。总体趋势是:植被盖度大,分形维数越高。研究区的样本土壤中黏粒含量均<15%,属于砂质土壤。这类土壤砂粒含量较高,透水性好,保蓄性差。灌木林地和荒草地砂粒含量高于灌木林地和林地,这是由于林地和灌木林地通过冠层截留降低降雨侵蚀力影响水土流失,加上枯枝落叶层对表层土壤的保护,表层土壤中细颗粒含量较多,砂粒含量较低。

研究区农用地主要种植玉米,由于人为干扰剧烈,虽然有施肥投入,但由于频繁土壤耕作和收获等因素影响,土壤表层疏松、水土流失严重,进而导致砂粒与粉粒的含量比较高,而黏粒含量比较低。灌木林地、林地和荒草地的土壤机械组成表现出较大的尺度差异,在小流域,灌木林地的砂粒含量占51.50%,远高于粉粒含量37.05%,而在集水区尺度上,砂粒含量和粉粒含量分别为44.21%和44.14%;在小流域,林地的砂粒含量为48.83%,高于粉粒含量40.38%,在集水区尺度,砂粒和粉粒含量分别为46.45%和43.17%;在小流域,荒草地的砂粒和粉粒含量分别为52.54%和36.51%,集水区分别为47.41%和41.87%。

黏粒含量的变异性较低,采样点中园地的黏粒含量最高,灌木林地和荒草地的黏粒含量高于林地和农用地。集水区林地、灌木林地以及荒草地发育较好,且在集水区西南和东北部有大量封育林场,荒草地植被覆盖度也较高。一方面由于山坡上部冲刷而来的细土粒在林地沉积,另一方面林地、灌木林地和荒草地表层植被覆盖度大,枯枝落叶和有机质含量丰富,腐殖化作用明显,具有明显的保持水土和改良土壤结构功能,所以黏粒含量较高。园地土壤机械组成的尺度差异不明显,这可能与小流域园地采样不足有关。

2.3 不同尺度表层土壤颗粒的分形维数与土壤粒径的关系

从分形维数的计算过程可知,D值与土壤颗粒粒径由大到小的累积含量有关,为了进一步探寻土壤颗粒体积分形维数与土壤颗粒组成之间的关系,对土壤颗粒体积百分含量与分形维数进行Pearson相关性分析。分析表明,在沐浴小流域和安塞集水区,土壤颗粒体积分形维数与土壤砂粒、黏粒百分含量呈极显著正相关关系(P<0.01),与土壤粉粒百分含量呈极显著负相关关系(P<0.01)。

表4 土壤分形维数D与粒径含量的相关性 Table 4 Correlation analysis between fractal dimension and soil particle size
尺度 ScalesD值与砂粒含量Pearson相关性 Pearson correlation between D-value and sand content D值与粉粒含量Pearson相关性 Pearson correlation between D-value and silt content D值与黏粒含量Pearson相关性 Pearson correlation between D-value and clay content
沐浴 Muyu0.836-0.9060.432
安塞Ansai0.737-0.8300.372

为深入揭示D值与各粒径的关系,将3个粒径含量与D值进行多元逐步回归分析。小流域D值与粉粒(x1)、砂粒(x2)、黏粒(x3)回归方程为D=5.918-0.5x1-0.37x2,模型检验达到极显著水平,R2=0.987。偏回归系数显著性检验结果表明,D与常数、砂粒、粉粒的偏相关关系显著,粉粒含量对D值的影响大于砂粒的影响。集水区D值的回归方程为D=2.161-0.12x2+0.41x3,模型检验达到极显著水平,R2=0.980。偏回归系数显著性检验表明,D与常数、粉粒、黏粒的偏相关关系显著,粉粒含量对D值的影响大于黏粒含量。

进一步分析D值与土壤粒径含量的关系,按照美国农业部土壤质地分类三角图标法,将土壤样品进行质地划分,采集样点的土壤类型为粉壤土、壤土和砂质壤土。对不同质地土壤D值进行统计分析。结果表明,在沐浴小流域,粉壤土、壤土和砂质壤土的D值平均数分别为1.872、2.100和2.255,安塞集水区分别为1.934、2.089和2.227。采集样品中,土壤颗粒黏粒含量平均值为10.74%,粒径含量以粉粒和砂粒为主。D值与砂粒含量正相关,砂质壤土D值高于壤土和粉壤土。这种变化在沐浴小流域和安塞集水区的表现是一致的。这与在陕北黄土丘陵沟壑区土壤粒径的相关研究结论是基本一致[12, 23, 24]的。但也有研究[14, 17]认为黏粒含量与D值相关性最强。在本研究中,土壤黏粒含量较低,且变异性小,因此与D值的相关性也较弱。

2.4 不同尺度下土壤粒径与环境因子的关系

利用单因素方差分析不同环境因子对粒径含量和分形维数D的相关性(表 5)。进一步进行CCA分析的结果见图 2

表5 土壤粒径、分形维数与环境因子的相关性分析 Table 5 Correlation analysis of soil particle size,fractal dimension and environmental factors
粒径与分形纬数 Soil particle size and fractal dimension尺度 Scales相对海拔 Relative altitude坡度 Slope 坡向 Aspect 植被盖度 Vegetation coverage 坡位 Slope position 土地利用方式 Land use types
D沐浴0.173-0.156-0.0670.1160.127-0.226
安塞-0.017-0.042-0.0360.2720.023-0.095
砂粒Sand 沐浴 0.0370.075-0.1190.0340.098-0.079
安塞-0.0180.011-0.0120.1250.048-0.123
粉粒Silt沐浴 -0.066-0.0130.114-0.047-0.1090.068
安塞0.023-0.0060.012-0.146-0.0400.099
黏粒clay沐浴 0.193-0.3270.0340.0900.0710.080
安塞-0.023-0.0310.0030.096-0.0590.169
图2 沐浴小流域 (上) 和安塞集水区 (下) 土壤颗粒和环境因子关系的二维排序图 Fig.2 The two-dimensional CCA ordination diagram of soil particle size,fractal dimension and environmental factors in the Muyu small watershed area (up) and the Ansai catchment (down)

CCA分析中,环境因子用箭头表示,箭头所处的象限表示环境因子与排序轴间的正负相关性,箭头连线的长度代表着某个环境因子与土壤粒径分布之间相关程度的大小,连线越长,说明相关性越大。反之越小。箭头连线和排序轴的夹角代表着某个环境因子与排序轴的相关性大小,夹角越小,相关性越高;反之越低。

图 2可知,D值与砂粒、黏粒含量正相关,且与砂粒相关性高于黏粒;D值与粉粒含量负相关。由环境因子之间的关系可知,在小流域坡度和土地利用类型对粒径和D值的影响表现出微弱负相关,即随着坡度增加,土地利用类型对土壤粒径的影响增大,但随着研究区由小流域变化到集水区,这种相关性消失。在小流域坡度和坡向对粒径和D值的影响表现出较强的负相关,即随着坡度增加,坡向对小流域的影响作用减弱;在集水区也表现出了相同的趋势。

根据土壤粒径与环境因子的关系(图 2)可知,在小流域,土地利用和地形因子对土壤粒径的总体影响差异较小,植被盖度、相对海拔和坡位相对于土地利用类型、坡向和坡度对土壤粒径影响较大;安塞集水区植被盖度、土地利用类型和坡度对土壤粒径影响较大。在两个研究区砂粒含量和坡度、坡位均表现为正相关,与坡向、土地利用类型负相关;粉粒则相反。且砂粒和粉粒含量均与相对海拔和植被盖度的相关性较弱。但随着研究区由小流域到集水区,各含量与各环境因子的相关性大小不同。小流域植被盖度对粒径的影响较大,这是由于植物枯枝落叶的分解能够增加土壤有机质,改善土壤结构,植物根系也能较好的保持土壤,改善其结构。但由于小流域面积小,土地利用类型差别不大,因此对土壤粒径影响较小。随着研究区由小流域到集水区,植被盖度和土地利用类型对土壤粒径的影响作用明显高于地形因子。这说明在集水区尺度,土地利用类型和植被盖度是影响土壤粒径的主要因素,这与很多研究的结果是一致的[11, 25]。地形因子中,坡位在小流域中对土壤粒径含量影响较大,坡度在集水区影响较大。小流域的坡位与土壤颗粒性质关系密切,研究区易发生土壤流失,从坡上冲刷下来的细土粒易沉积在坡下[25];以往研究[26]也表明,坡下部沉积带比坡上部侵蚀带的粉粒含量高,砂粒和黏粒含量低。在较大的集水区尺度,相对海拔对集水区的影响很小,这是由于研究区为丘陵沟壑区、地形破碎,相对海拔和坡度等地形因子共同影响局地水热平衡及土壤侵蚀过程,坡度因子对粒径的影响占主导作用。

3 结论

(1)沐浴小流域的土壤砂粒、粉粒和黏粒的含量分别为50.27%,38.83%和10.89%;安塞集水区表层土壤砂粒、粉粒和黏粒的含量分别为46.20%,43.06%和10.74%。沐浴集水区D值分布在1.86—2.55之间,平均值为2.14,安塞集水区D值分布在1.71—2.55之间,平均值为2.09。研究区由沐浴小流域到安塞集水区,各粒径百分含量和D值的空间变异性增强。在小流域尺度,土壤质地的空间变异性较小,而在集水区尺度,土壤属性的变异性有所增加。

(2)沐浴小流域不同土地利用类型的砂粒含量从高到低依次为:荒草地>农用地>林地>灌木林地>园地,集水区的次序为:荒草地>灌木林地>林地>农用地>园地。粉粒含量的次序均与砂粒相反。不同土地利用类型土壤粒径分形维数D从高到低排序,在沐浴小流域表现为灌木林地>荒草地>林地>园地>农用地,在安塞集水区依次为灌木林地>荒草地>农用地>林地>园地。

(3)小流域D值与粉粒(x1)、砂粒(x2)、黏粒(x3)回归方程为D=5.918-0.5x1-0.37x2,集水区D值的回归方程为D=2.161-0.12x2+0.41x3D值与粉粒极显著负相关,与黏粒正相关,且随着研究区尺度的增加,正相关性增强。按质地进行分析,在沐浴小流域,粉壤土、壤土和砂质壤土的D值平均数分别为1.872、2.100和2.255,安塞集水区分别为1.934、2.089和2.227。采集样品中,土壤颗粒黏粒含量平均值为10.74%,粒径含量以粉粒和砂粒为主。D值与砂粒含量正相关,砂质壤土D值高于壤土和粉壤土。这种变化在沐浴小流域和安塞集水区的表现是一致的。

(4)土壤性质的变异性是不同尺度地形、土壤母质、土地利用及人为活动等多种因子综合作用的结果。黄土丘陵沟壑区地貌类型多样、地形破碎,不同坡位、坡向、海拔影响水热条件和土壤物质的再分配。其中,在小流域,植被盖度、相对海拔和坡位对土壤粒径影响较大;在安塞集水区植被盖度、土地利用类型和坡度对土壤粒径影响较大。

参考文献
[1] 傅伯杰, 陈利顶, 马克明, 王仰麟. 景观生态学原理及应用. 北京: 科学出版社, 2001: 16-41, 202-236.
[2] 肖笃宁, 李秀珍, 高峻, 常禹, 李团胜. 景观生态学. 北京: 科学出版社, 2003: 5-6, 130-131.
[3] Wu J G. Hierarchy and scaling: extrapolating information along a scaling ladder. Canadian Journal of Remote Sensing, 1999, 25(4): 367-380.
[4] 傅伯杰,赵文武,张秋菊,刘宇.黄土高原景观格局变化与土壤侵蚀.北京:科学出版社, 2014:2-4
[5] Giménez D, Perfect E, Rawls W J, Pachepsky Y. Fractal models for predicting soil hydraulic properties: a review. Engineering Geology, 1997, 48(3/4): 161-183.
[6] Huang G H, Zhang R D. Evaluation of soil water retention curve with the pore-solid fractal mode. Geoderma, 2005, 127(1-2): 52-61.
[7] Montero E. Rényi dimensions analysis of soil particle-size distributions. Ecological Modelling, 2005, 182(3/4): 305-315.
[8] 王德, 傅伯杰, 陈利顶, 赵文武, 汪亚峰. 不同土地利用类型下土壤粒径分形分析-以黄土丘陵沟壑区为例. 生态学报, 2007, 27(7): 3081-3089.
[9] Gui D W, Lei J Q, Zeng F J, Mu G J, Zhua J T, Wang H, Zhang Q. Characterizing variations in soil particle size distribution in oasis farmlands-A case study of the Cele Oasis. Mathematical and Computer Modelling, 2010, 51(11-12): 1306-1311.
[10] Xu G C, Li Z B, Li P. Fractal features of soil particle-size distribution and total soil nitrogen distribution in a typical watershed in the source area of the middle Dan River, China. Catena, 2013, 101: 17-23.
[11] 胡云锋, 刘纪远, 庄大方, 曹红霞, 闫慧敏. 不同土地利用/土地覆盖下土壤粒径分布的分维特征. 土壤学报, 2005, 42(2): 336-339.
[12] 董莉丽, 郑粉莉. 陕北黄土丘陵沟壑区土壤粒径分布分形特征. 土壤, 2010, 42(2): 302-308.
[13] 许明祥, 刘国彬, 赵允格. 黄土丘陵区土地利用及环境因子对土壤质量指标变异性的影响. 应用生态学报, 2011, 22(2): 409-417.
[14] 郭中领, 符素华, 王向亮, 张学会. 北京地区地表土壤分形特征研究. 水土保持通报, 2010, 30(2): 154-158.
[15] 曾宪勤, 刘和平, 路炳军, 王秀颍, 杨威. 北京山区土壤粒径分布分形维数特征. 山地学报, 2008, 26(1): 65-70.
[16] 缪驰远, 汪亚峰, 魏欣, 徐霞, 石文. 黑土表层土壤颗粒的分形特征. 应用生态学报, 2007, 18(9): 1987-1993.
[17] 张秦岭, 李占斌, 徐国策, 张铁钢, 黄萍萍, 张洋. 丹江鹦鹉沟小流域不同土地利用类型的粒径特征及土壤颗粒分形维数. 水土保持学报, 2013, 27(2): 244-249.
[18] 王富, 贾志军, 董智, 李红丽, 甄宝艳, 孙明高, 张卫平. 不同生态修复措施下水库水源涵养区土壤粒径分布的分形特征. 水土保持学报, 2009, 23(5): 113-117.
[19] 赵明月, 赵文武, 钟莉娜. 土地利用和环境因子对表层土壤有机碳影响的尺度效应分析: 以陕北黄土丘陵沟壑区为例. 生态学报, 2014, 34(5): 1105-1113.
[20] 杨培岭, 罗远培, 石元春. 用粒径的重量分布表征的土壤分形特征. 科学通报, 1993, 38(20): 1896-1896.
[21] Leps J, Smilauer P. Multivariate Analysis of Ecological Data Using CANOCO. Cambridge: Cambridge University Press, 2003.
[22] 张金屯. 数量生态学. 北京: 科学出版社, 2003.
[23] 邓良基, 林正雨, 高雪松, 张世熔. 成都平原土壤颗粒分形特征及应用. 土壤通报, 2008, 39(1): 38-42.
[24] 苏永中, 赵哈林. 科尔沁沙地农田沙漠化演变中土壤颗粒分形特征. 生态学报, 2004, 24(1): 71-74.
[25] 邱扬, 傅伯杰, 王军, 陈利顶. 黄土丘陵小流域土壤物理性质的空间变异. 地理学报, 2002, 57(5): 587-594.
[26] Martined-Mena M, Williams A G, Ternan J L, Fitzjohn C. Role of antecedent soil water content on aggregates stability in a semi-arid environment. Soil & Tillage Research, 1998, 48(1-2): 71-80.