南长山岛土壤理化因子空间分布特征

徐苏源， 衣华鹏， 闫潍虹， 位 彬， 张小梅， 张安安

(鲁东大学 资源与环境工程学院, 山东 烟台 264025)

根据《联合国海洋法公约》的定义，岛屿是四面环水并在高潮时高于水面的自然形成的陆地区域。因海岛四周被海水环绕，地理环境特殊，其自然环境与陆地相比具有相对独立性、独特多样性和脆弱易变性，其生态系统也兼有陆地、湿地和海洋的特征[1]。从地理区域上来看，目前国内学者对全国各地土壤的理化性质研究较广泛。其中，对内陆湖泊湿地的研究[2-3]揭示了湿地湖区土壤含水量和有机质含量的正相关关系，对河岸带土壤进行的分析[4]表明河岸带土壤有机质含量的空间分异明显且主要积聚在地表。有的探讨了植被恢复模式[5]、沙漠化过程[6]、森林采伐方式[7]或轮作方式[8]对土壤理化因子的影响，结果体现出不同生境中上述因子和土壤理化性质有密切的相关性。同时从各地海岛土壤的研究工作[9-11]中可以看出，海岛土壤研究主要采用了正态分布检验、相关性分析、方差分析和多重比较等地统计学手段，目前土壤理化性质的研究已经有了较完整的思路和方法。文献表明，南长山岛的乔木碳储量高于山东省平均值[12]，地形、土壤和乔木对草本植物多样性有显著影响[13]，也有对南长山岛海岸潮间带[14]进行的研究，但是对南长山岛的土壤理化性质研究较少。南长山岛是庙岛群岛中离大陆最近，受人类活动影响最大的岛屿。近年来随着长岛旅游资源的开发和旅游业的发展，越来越多的游客涌入南长山岛，对南长山岛的生态产生了巨大的压力。土壤是植物生长的重要基础，其养分含量对植物的数量和质量有明显影响，其理化性质是了解海岛的土壤特征及演变规律的重要基础。通过研究南长山岛的土壤理化性质，揭示不同生境下土壤的空间异质性和相关性，可以为地表植被的保护和有序开发提供科学依据，为海岛生态保护提供参考。

1 材料与方法

1.1 研究区概况

南长山岛，行政上隶属于山东省烟台市长岛县，位于胶东、辽东半岛之间，黄渤海交汇处，南距蓬莱7 km，位于37°53.43′—37°57.39′N，120°42.92′—120°45.59′E，为庙岛群岛最南端。南长山岛岛屿面积为13.21 km2，海岸线长21.6 km，属于低山丘陵地形，整体上东高西低，最高海拔为156.1 m；气候属于东亚温带季风型气候，年均温12 ℃，最冷月均温-1.6℃，最热月均温24.5 ℃，年均降水量549.5 mm，降水集中在6—9月，≥10 ℃积温在4 000～4 600 ℃之间，年日照时数在2 500～2 800 h之间，岛上无较大的河流和湖泊，岩石类型主要为石英岩和板岩，土壤主要为棕壤和褐土，植被类型主要为常绿针叶林、落叶阔叶林、灌丛、草丛、滨海盐生沙生植被，以及少量果林和农作物分布。森林覆盖率为44.8%，乔木以黑松和刺槐为优势种。

1.2 样带设置和样品采集

2016年10月19—21日对南长山岛进行了实地考察，主要在岛东部丘陵区采集土壤样本，使用手持GPS进行定位，沿山脊往各个方向共设置了15个样带、70个样方，采集表层20 cm内土壤，每个样方内采集四角和中间土壤混合，得到有效土壤样品68份。用环刀采集了土壤，密封后带回实验室，作为容重待测样品。

1.3 试验及数据处理

各样方的土壤指标值是3个重复样品的平均值。土壤含水量和容重使用环刀法，采用靛酚蓝比色法测定铵态氮，采用四苯硼钠比浊法测定速效钾，采用浸提—钼锑抗比色法测定速效磷，采用重铬酸钾法测定有机质，pH值采用5∶1水土比的酸度计测定。土壤养分标准参照我国第二次土壤普查分级表[15]，使用Excel 2010对数据进行统计处理，使用SPSS Statistics 24进行统计分析、正态分布检验、方差分析和多重比较。

2 结果与分析

2.1 土壤理化性质的描述性统计分析

对采集的28个样本进行描述性统计分析，结果详见表1。由土壤养分等级分级标准可知，样点的铵态氮、速效磷和有机质属于四级，土壤较贫瘠。速效钾属于二级，土壤钾含量丰富。大部分土壤为棕壤，由于淋溶作用较强，加上成土母质是酸性岩风化物，导致土壤呈现微酸性，pH值大部分在5.0～6.0之间；少部分发育在黄土母质上的土壤为褐土，pH值大于6.5。土壤容重大部分在1.2～1.6之间，部分土壤中砾石含量较多。南长山岛植被群落中的乔木以刺槐和黑松为优势种，多为20世纪60—70年代人工种植，部分采样点土壤被深厚的枯枝落叶覆盖，腐殖质层较厚。

对各项指标计算变异系数，变异强度由小到大依次为pH值、容重、速效磷、砾石占比、有机质、速效钾、含水量、铵态氮。根据土壤变异系数的划分等级，Cv<0.1为弱变异性；0.11为强变异性[16]。据此，除了pH值属于弱变异性，其他项目都是中变异性；其中容重、速效磷、砾石占比、有机质变异系数较小，范围在0.3～0.5之间；速效钾、含水量、铵态氮变异系数较大，范围在0.5～0.95之间。土壤变异系数可以表示土壤各项肥力指标的离散程度，其大小受地表植被生长状况、植被种类、土壤质地等因素的影响。铵态氮，速效钾、含水量变异系数较大，说明其在空间上分布不均，推测因为受植被长势特别是乔木长势的空间差异与地形坡度差异的影响，导致土壤养分含量呈现空间分布上的不均衡。对土壤理化性质的各项指标进行正态分布检验，从峰度和偏度的统计数值可以看出，土壤的容重、速效磷、有机质和pH值服从正态分布，含水量、砾石占比、铵态氮和速效钾经过对数转换后服从对数正态分布。说明采样数据具有代表性，符合统计学特征。

注：*原始数据经过对数转换后计算得出。

2.2 土壤理化因子的相关性分析

对土壤的含水量、容重、铵态氮、速效钾、速效磷、有机质和pH值进行相关性分析，结果详见表2。计算得容重与含水量呈显著负相关(r=-0.720)；土壤有机质与含水量(r=0.419)和速效钾(r=0.252)呈显著正相关，与铵态氮(r=0.053)和速效磷(r=0.029)呈不明显正相关。容重越大土壤的孔隙度就越小，土壤的保水保肥能力越弱，土壤的含水量就越小；同时随着土壤含水量的降低，植被长势就差，地表枯枝落叶层就较薄，土壤有机质积累较少，导致土壤养分含量降低。南长山岛土壤容重和含水量是影响有机质含量的因素之一；容重和土壤含水量的变化，对速效钾影响明显，对铵态氮和速效磷影响不明显。南长山岛面积较小，无明显的河流和湖泊发育，土壤水的来源主要是大气降水，而土壤容重的变化导致了含水量的变化，由此影响了地表植被的生长发育。因此含水量较高的土壤地表植被的长势更好，有机质的积累也更多。实地调查中，土壤水含量高的地方往往覆盖着较厚的凋落物层，部分样点的凋落物层达到了10 cm，腐殖质层较厚，土壤分层明显，有机质的积累也促进了土壤养分的积累。

2.3 土壤理化因子的空间异质性分析

由于本研究区的面积较小，各样地海拔高度差异不大，由海拔变化所引起的水热条件变化不明显，因此对采集的样点土壤按坡向和坡位进行分组。研究表明，影响南长山岛土壤理化因子空间分布的因素有地形、植被和人类活动等，其中地形和植被是主要的影响因素。

表2 土壤各指标相关系数

注：**表示在0.01级别相关性显著；*表示在0.05级别相关性显著。

2.3.1 不同坡位土壤与植被的异质性分析 根据表3可知，经过方差分析，除了铵态氮，其他各项指标在不同坡位表现出显著差异性。从多重比较结果看，差异主要存在于坡顶和坡中，坡顶的土壤有机质和速效钾含量更高，速效磷和pH值略小，而坡底的指标值介于两者之间。根据表4可知，草本覆盖度、乔木高度、灌木高度、灌木丛径有显著性差异，坡顶和坡中差异最大，如乔木高度和灌木丛径。坡顶的乔木长势略差于坡中和坡底，但灌木长势相对较好。这与土壤呈现的规律相一致。

2.3.2 不同坡向土壤与植被的异质性分析 根据表5，不同坡向铵态氮和速效磷表现出显著差异。多重比较表明：铵态氮南坡偏低，速效磷北坡偏高，其他指标无显著差异。

根据表6可知，有显著差异的有草本覆盖度、乔木胸径和乔木郁闭度。其中西坡草本覆盖度较高，东坡乔木胸径和乔木高度较高，这说明西坡草本长势更好，东坡乔木长势更好。

表3 不同坡位土壤统计分析

注：数值为平均值±标准差，相同字母表示无显著性差异(S-N-K多重比较p<0.05)。下同。

表4 不同坡位植被统计分析

表5 不同坡向土壤统计分析

表6 不同坡向植被统计分析

3 讨 论

3.1 不同坡位土壤与植被的异质性讨论

物种在群落中的分布格局是物种与生境相互作用的结果，不仅与生境中的地形因子及土壤因子有密切关系，也受到物种的生态学特性及其与周边物种间关系的影响[17]。乔木和灌木生长状况与草本生长状况有明显的负相关性。乔木和灌木的高度较高，郁闭度更大，其枝干和叶片遮挡住阳光，使地表有效日照时间变短，有效日照面积缩小，影响了草本植物的生长。因此，坡顶灌木长势好而草本植物长势差。土壤的发育过程受到母质、气候、地形和植被等自然因素的影响，同时与土地利用、农业生产等人类活动也有着密切关系。这些成土因素在空间上的分布都不是均一的，在不同的空间尺度上表现出普遍的空间变异性。植被对土壤养分空间分布有明显的影响，坡顶的乔木郁闭度、灌木高度和郁闭度均比坡中高，速效钾和有机质也比坡中高，土壤养分和植被长势呈现出相关性。由于坡顶乔木和灌木长势更好，更多的凋落物进入土壤中，经过微生物分解并长年累积，有机质含量也随之升高。因此乔木长势好的地方，土壤中有机质含量也相对较高。人类活动对土壤养分和植被长势也有明显影响。由于坡顶和坡中位置属于长岛国家级自然保护区和长岛旅游景区的范围，而坡底靠近公墓、工厂、住宅、公路等生产生活区域，坡底人口密度和人类活动强度远远高于坡顶和坡中。因此坡底的植被和土壤受人类活动影响更大，生物圈循环受到影响，导致坡底一部分指标不呈现任何规律，土壤理化因子的空间差异性主要体现在自然状态下的坡顶和坡中。

3.2 不同坡向土壤与植被的异质性讨论

南坡属于阳坡和迎风坡，北坡属于阴坡和背风坡。受日照和降水差异影响，南坡的土壤受淋溶作用更强，水分蒸发强烈，水土流失比北坡严重，土壤与养分不如北坡，导致铵态氮南坡偏低，速效磷北坡偏高。坡位的变化使不同种类的植被产生空间差异。西坡在冬季属于西北风的迎风坡，面向大海，受强烈的冬季风影响，乔木长势不如东坡好，但草本植物长势相对较好；而东坡属于冬季风背风坡，受冬季风影响小，乔木长势较好，同时草本长势不如西坡，也印证了在坡位上乔木和草本植物长势的相互影响关系。由于本次研究采样点较少且样地在全岛范围内未均匀分布，本文所得结论还需进一步采样试验验证。

4 结 论

(1) 南长山岛大部分土壤呈酸性和微酸性，研究区域土壤较贫瘠，存在土壤结构差和养分不协调等问题。pH值变异系数较小。南长山岛土壤含水量是影响有机质含量的重要因素之一。

(2) 地形因素主要通过坡位和坡向对土壤理化因子空间分布产生影响。土壤差异主要体现在不同坡位上，植被差异主要体现在不同坡向上。

(3) 植被因素的影响体现在两个方面。乔木生长对土壤环境产生的作用，乔木长势好的地段土壤养分含量高。乔木和灌木长势影响林下草本植物长势，乔木和灌木长势好的地段草本植物长势较差，两者呈现负相关关系。

(4) 人类活动总体上为次要影响因素，但在自然环境整体格局的基础上对局部产生影响；坡底受人类活动影响较大，坡底的土壤和植被没有体现出明显的规律性。

致谢：除本文的作者外，韩春梅、李苹苹也参加了野外考察工作，潘蕊蕊、刘伟杨、张梦杰、韩亚婷、亢琼琼也参加了本研究的试验工作，在此表示感谢。

[参考文献]

[1] 石洪华,郑伟,丁德文,等.典型海岛生态系统服务及价值评估[J].海洋环境科学,2009,28(6):743-748.

[2] 秦璐,傅德平,杨军,等.艾比湖湿地自然保护区典型群落土壤理化性质分析[J].干旱区资源与环境,2011,25(8):144-151.

[3] 闫培锋,周华荣,刘宏霞.白杨河—艾里克湖湿地土壤理化性质的空间分布特征[J].干旱区研究,2008，25(3):406-412.

[4] 李冬林,韩丽,阮宏华,等.秦淮河河岸带土壤理化性质分析[J].南京林业大学学报：自然科学版,2008，32(4):17-22.

[5] 王昭艳,左长清,曹文洪,等.红壤丘陵区不同植被恢复模式土壤理化性质相关分析[J].土壤学报,2011,48(4):715-724.

[6] 杨梅焕,曹明明,朱志梅,等.毛乌素沙地东南缘沙漠化过程中土壤理化性质分析[J].水土保持通报,2010,30(2):169-172，176.

[7] 刘美爽,董希斌,郭辉,等.小兴安岭低质林采伐改造后土壤理化性质变化分析[J].东北林业大学学报,2010,38(10):36-40.

[8] 万福绪,陈平,王严星.苏北林粮间作地土壤理化性质分析[J].南京林业大学学报：自然科学版,2003，27(6):27-30.

[9] 程芳,程金平,于金莲,等.大金山岛土壤有机质与全氮的空间分布及其相关性[J].环境化学,2013,32(9):1634-1638.

[10] 孙玉冰,邓守彦,李德志,等.崇明县土壤主要理化指标的空间分布与变异特征[J].生态与农村环境学报,2010,26(4):306-312.

[11] 刘骄.福建海坛岛主要防护林树种的生物生产力研究和土壤养分特征分析[D].福建 福州：福建师范大学,2011.

[12] 王晓丽,王嫒,石洪华,等.山东省长岛县南长山岛黑松和刺槐人工林的碳储量[J].应用生态学报,2013,24(5):1263-1268.

[13] 池源,郭振,石洪华,等.南长山岛草本植物多样性及影响因子[J].华中师范大学学报：自然科学版,2015,49(6):967-978.

[14] 庄树宏,陈礼学,孙力.南长山岛岩岸潮间带底栖藻类群落结构的季节变化格局[J].海洋科学进展,2003，21(2):194-202.

[15] 全国土壤普查办公室.中国土壤[M].北京:中国农业出版社,1998:356.

[16] 雷志栋,杨诗秀,许志荣,等.土壤特性空间变异性初步研究[J].水利学报,1985(9):10-21.

[17] 叶诺楠,沈娜娉,商天其,等.浙江瑞安公益林群落结构及其与环境的相关性[J].植物学报,2017,52(4):496-510.