研究区设在江苏省盐城市境内的东台林场，其地理位置为32°48′40″N，120°49′31″E。东台林场地处江苏省中东部黄海之滨，创建于1965年，属于暖温带和亚热带的过渡区，四季分明，常年平均气温为14.6 ℃，年均无霜期为225.0 d，年均降水量为1 051.0 mm，年均日照时数为2 169.6 h。土壤类型为脱盐草甸土，质地为砂质壤土，pH值偏碱性。试验样地于2012年2月开始设立，通过人工外源氮输入模拟未来氮沉降趋势。根据经营管理措施和立地条件基本相同的原则，在林场内选择11年生中林龄杨树Populus deltoides ‘I-35’人工林分。①样地面积：大小为20 m × 90 m，3个重复。②样地内样方设置：各个重复样地包括5块10 m × 20 m样方（随机区组实验方法排列），样方之间缓冲带宽为10 m，总面积为20 m × 90 m。③实验处理：样方施氮量梯度处理为N₀（施氮0 g·m^-2·a^-1，对照），N₁（施氮5 g·m^-2·a^-1），N₂（施氮10 g·m^-2·a^-1），N₃（施氮15 g·m^-2·a^-1），N₄（施氮30 g·m^-2·a^-1），在每年的生长季节5-10月进行施氮处理，1次·月^-1，共施6次·a^-1[16-17]。

选用2015年期间的数据。在实验样地内，于4，6，8，10，12月在各个施氮样方内随机选取5个采样点，用土钻取0~10 cm层土壤。土壤样品带回实验室后，将同一个施氮样方内的5袋土壤样品充分混合成1袋，共计75袋。将新鲜土样分为2份，一份去杂后，过2 mm的钢筛后进行SMBC的测定；另一份自然风干、去杂、过2 mm筛后进行土壤总有机碳（TOC）以及pH值等其他指标的测定。土壤TOC采用岛津TOC-VCPH分析仪测定^[18]；土壤微生物生物量碳采用氯仿熏蒸-硫酸钾浸提法测定^[19]；土壤pH值采用m（土）:m（水）=1.0:2.5电位法测定^[20]；土壤硝态氮（NO₃^--N）采用双波长紫外分光光度法测定^[21]。

利用Origin 8.5，SPSS 19.0和Excel 2016等进行数据分析和图表处理。采用重复测量方差分析和单因素方差分析不同施氮水平土壤TOC和SMBC的差异显著性，并对不同月份、不同量氮输入下的土壤TOC质量分数、SMBC质量分数做多重比较（显著性水平为0.05）。

HAUSDORFF在1919年提出了连续空间的概念，也就是空间维数是可以连续变化的，它可以是整数也可以是分数，称为豪斯道夫维数。设一个整体S划分为N个大小和形态完全相同的小图形，每一个小图形的线度是原图形的r倍, 则豪斯道夫维数D=lim[logN（r）/log（1/r）]。计算的基本原理为分形集都遵循一定的标度律，即测度随测量尺度按照一种幂指数规律而变化，即在双对数坐标中作图，并运用最小二乘法拟合一条直线，其斜率k与分形维数D之间有如下关系：D=f（k）。采用不同的测度，对应的函数D=f（k）也不相同，如利用变异函数法（semivariogram，SV）和根据功率谱密度法计算分维（power spectrum density，PSD）^[22]等。本研究采用PSD法分析土壤TOC和SMBC随月份变化的分形关系以及土壤TOC与SMBC的分形关系。分形PSD曲线具有下列表达幂函数关系：

式（1）中：f是频率，在本研究中代表月份；S（f）是PSD，在本研究中代表土壤TOC，SMBC；w是PSD曲线的线性回归所得回归直线的斜率。斜率w与分形维数D的关系为：

式（2）中：D定量表征了土壤TOC和SMBC随月份变化的复杂程度以及土壤TOC质量分数随SMBC质量分数变化的复杂程度^[23]。

氮输入对土壤TOC质量分数影响显著（P＜0.05），氮输入在不同月份对土壤TOC质量分数的影响不同（图 1）。经多因素方差分析，在4，6，8，10，12月氮添加对土壤TOC质量分数影响极显著（P＜0.01）。氮输入处理下土壤TOC质量分数在生长季6-8月显著降低，在非生长季10-12月、12-4月显著增加。进一步对不同量氮输入水平下土壤TOC质量分数做多重比较，发现不同量氮输入水平下土壤TOC质量分数变化差异极显著（P＜0.01）。

图  1  氮输入下杨树人工林土壤TOC的时间动态

Figure 1.  Dynamics of total organic carbon in poplar plantations under nitrogen treatment

研究区不同施氮水平下土壤TOC质量分数随时间变化的分形维数（D）进行分析（表 1）。由表 1看出：N₀施氮水平下土壤TOC质量分数与时间变化呈线性正相关，但未达到显著水平；N₁，N₃和N₄施氮水平下土壤TOC质量分数与时间变化呈线性显著正相关；N₂施氮水平下土壤TOC质量分数与时间变化呈线性极显著正相关，这与彭赛^[24]和肖晗冉^[25]的研究一致。研究区不同施氮水平下土壤TOC质量分数随时间变化的D变化范围是1.805~1.949，土壤TOC质量分数随时间变化的分维均很接近，且接近于2，同布朗粒子运动轨迹的维数2很接近，表明不同施氮水平下土壤TOC质量分数随时间变化具有较大的随机性和复杂性^[26]。不同施氮水平下的D从大到小依次为N₃，N₂，N₄，N₁，N₀，在N₃和N₂施氮水平下，D数值较大，在这2个施氮水平下土壤TOC质量分数随时间变化的复杂程度更高，说明中等施氮量处理对土壤TOC质量分数随时间变化的影响更显著；N₁施氮水平下的D值小于N₂，N₃和N₄施氮水平下的D值，说明N₁施氮水平下土壤TOC质量分数随时间变化的复杂程度相对较小，此施氮水平处理对土壤TOC质量分数随时间变化的影响较微弱；N₄施氮水平下D值于N₂和N₃施氮水平下的D值，说明N₄施氮水平下土壤TOC质量分数随时间变化的复杂程度相对较小，此施氮处理对土壤TOC质量分数随时间变化的影响较微弱。

氮输入对SMBC质量分数影响极显著（P＜0.01），氮输入在不同月份对SMBC质量分数的影响不同（图 2）。经多因素方差分析，6和10月氮输入对SMBC无显著影响（P＞0.05）。4，8和12月氮添加对SMBC质量分数影响极显著（P＜0.01）。氮添加处理下SMBC质量分数在生长季6-8月显著增加，在8-10月生长季末其质量分数又显著降低，在12月达到峰值。进一步对不同氮添加水平下SMBC质量分数做多重比较，发4个处理下SMBC质量分数变化差异极显著（P＜0.01）。

图  2  氮输入下杨树人工林SMBC的时间动态

Figure 2.  Dynamics of soil microbial biomass carbon in poplar plantations under nitrogen treatment

对研究区不同氮水平下SMBC质量分数随时间变化的D进行分析（表 2）。由表 2看出：N₀，N₁，N₃，N₄施氮水平下SMBC质量分数与时间变化均呈线性正相关，且达到显著水平，不同施氮水平之间的显著性水平差异不大。研究区不同施氮水平下SMBC质量分数随时间变化的D变化范围是1.728~1.963，且数值接近于2，说明SMBC质量分数随时间变化具有较大随机性和复杂性^[27]。不同施氮水平下的D从大到小依次为N₂，N₃，N₁，N₄，N₀，N₂和N₃施氮水平下的D的值较大，在这2个施氮水平下SMBC质量分数随时间变化的复杂程度更高，说明中等施氮水平处理对SMBC质量分数随时间变化的影响更显著；N₄施氮水平下SMBC质量分数随时间变化的D值小于N₂和N₃施氮水平下的D值，说明此施氮处理对SMBC质量分数随时间变化的影响较微弱；在N₁施氮水平下SMBC质量分数随时间变化的D值小于N₂，N₃和N₄施氮水平下的D值，说明此施氮水平下SMBC质量分数随时间变化的复杂程度较小，对SMBC质量分数随时间变化的影响较微弱。

从表 3可见：在时间的动态变化作用下，研究区不同施氮水平下土壤TOC质量分数随SMBC质量分数变化的D变化范围是2.207~2.342，D值普遍较大，表明不同施氮水平下土壤TOC质量分数随SMBC质量分数变化具有较大的随机性和复杂性^[28]。不同施氮水平下的D从大到小依次为N₃，N₂，N₄，N₁，N₀。在N₃和N₂施氮水平下，D值较大，在这2个施氮水平下土壤TOC质量分数随SMBC质量分数变化的复杂程度更高，说明中等施氮水平处理对土壤TOC质量分数随SMBC质量分数变化的影响更显著；N₁施氮水平下的D值小于N₂，N₃和N₄施氮水平下的D值，说明N₁施氮水平下土壤TOC质量分数随SMBC质量分数变化的复杂程度相对较小，此施氮处理对土壤TOC质量分数随SMBC质量分数变化的影响较微弱；N₄施氮水平下D值小于N₂和N₃施氮水平下的D值，说明N₄施氮水平下土壤TOC质量分数随SMBC质量分数变化的复杂程度相对较小，此施氮处理对土壤TOC质量分数随SMBC质量分数变化的影响较弱。

由表 4可知：在不同施氮水平处理下，不同月份下土壤TOC质量分数随SMBC质量分数变化的D值从大到小依次为6，10，4，8，12月。研究区不同月份下土壤TOC质量分数随SMBC质量分数变化的D值变化范围是1.650~6.149，D值差异较大，说明土壤TOC质量分数随SMBC质量分数变化的随机性和复杂性和月份有很大的相关性^[23]。

土壤的形态和演化过程都非常复杂，想要精确恰当地描述和阐释土壤属性的时空特征，运用一般的变异函数达不到对土壤时空特征定量化描述的水平^[29]。分形理论运用不同于传统技术的空间分析理论，其核心分析理念自相似理论的运用更为研究土壤各种属性的时空变异提供了一种新颖的方法^[30]。将要分析和量化的土壤理化指标利用分形维数来分析其时空特征具有简单可行的特点。

本研究结果显示：N₁氮输入水平下土壤TOC质量分数随SMBC质量分数变化的D值小于N₂，N₃和N₄氮输入水平下的D值，此氮输入水平对土壤TOC质量分数随SMBC质量分数变化的影响并不显著，而N₂和N₃氮输入水平对土壤TOC质量分数随SMBC质量分数变化的影响极显著。N₂，N₃中等施氮水平处理下，极显著影响了土壤SMBC及其代谢强度，这可能是因为中等施氮水平处理促进了林地植被、地表草本植物和灌木的生长，林木凋落物和根系产物输入上升，增加了有机质的输入，为土壤微生物提供了更多的能源使微生物群落的生物量增加并且提高了其代谢强度。这与许多研究结果相一致。门中华^[31]在研究不同的硝态氮供应水平下冬小麦Triticum aestivum植株对氮素的利用水平时发现：中等氮输入下小麦植株具有最高的根系活力及氮素同化能力，这主要是由于中等水平的氮输入能够提高土壤微生物代谢强度。马慧君等^[32]在模拟氮沉降对杨树人工林土壤微生物优势种群结构影响的研究中，也得出了中等水平的氮输入能够提高土壤微生物代谢强度的结果，说明氮施加的水平在土壤肥力的增加和植物生长具有重要意义。N₄施氮水平下土壤TOC质量分数随SMBC质量分数变化的D值小于N₂和N₃氮输入水平下的D值，并且N₄氮输入处理下的土壤TOC质量分数随时间变化的D值和SMBC质量分数随时间变化的D值均小于中等水平N₂和N₃，这说明过量氮输入可能会降低林场植物的生长量、土壤微生物生物量及其活性。王晓荣等^[33]在中亚热带栎属Quercus不同树种幼苗的生长和生物量分配对短期氮沉降的响应的研究中，发现高水平氮输入对生物量积累产生了一定的抑制作用，导致这种现象出现的原因是高氮处理植株由培养介质中吸收的氮量、植株吸氮量、根系活力、营养液pH值变化均介于中氮与低氮处理之间。门中华^[31]在不同的硝态氮供应水平下冬小麦植株对氮素的利用水平的研究中也得出了这一结论。

本研究中，不同月份土壤TOC质量分数随SMBC质量分数变化的D从大到小依次为6，10，4，8，12月，在6月D值最大，此时土壤TOC质量分数随SMBC质量分数变化最具随机性和复杂性。周义贵等^[34]在川西亚高山地区米亚罗林区研究发现云杉Picea asperata低效林土壤TOC和SMBC质量分数变化均在夏季达到显著水平，导致这种现象的原因是研究区6月气温较高，水分充足，林场植被根系生长旺盛，各种生命活动增加，林木凋落物及根系产物输入上升，使得6月土壤TOC最高，此时土壤微生物新陈代谢和各种生命活动比较旺盛，对土壤TOC的作用强烈。10月植被处于生长期的末期，植被生长利用了大量的养分，使得林场土壤TOC质量分数出现了降低的趋势。KALBITZ等^[35]在研究生物降解诱导土壤溶性有机物性质变化的实验中发现：秋季的土壤TOC质量分数出现了降低的趋势，这主要是由于此时研究区湿热多雨的气候特点，土壤微生物的作用依旧比较强烈，一直在减少的土壤TOC质量分数在土壤微生物的强烈作用下，两者质量分数的变化呈现出复杂和随机的特性。4月土壤微生物各项生命活动开始加强，但损耗较高，此时土壤TOC输入量也比较小，所以D值较小；8月虽然此时土壤微生物种群数量比较大，但大量土壤TOC输入下，土壤微生物的对土壤TOC作用强度比6，10和4月微弱；12月D最小，土壤TOC质量分数随SMBC质量分数变化受到外界影响较小，随机性和复杂性降低。周莉^[36]在岩溶环境下土壤活性有机碳和土壤呼吸动态变化的研究中发现：在冬季时土壤TOC和SMBC质量分数受外界影响较少并出现最低值，这主要是由于冬季气温全年最低，微生物活性降低，故SMBC和TOC质量分数均处于全年较低水平。

本研究应用分形理论研究了苏北沿海地区杨树人工林不同水平氮输入下土壤TOC和SMBC之间质量分数变化的关系以及土壤TOC质量分数、SMBC质量分数随时间变化的情况。运用分形理论分析阐释不同施氮水平下土壤不同形式碳之间的相互关系以及随月份变化特征简单有效，对林业生产上制定合理的施氮策略具有很大的潜力和应用前景。从现有研究成果来看，分形理论的确提出了量化土壤属性空间分布特征的新思路，可以成为土壤时空变化研究的重要理论基础^[37]。