栓皮栎人工林土壤有机碳空间分布特征及其影响因素

doi:10.13466/j.cnki.lyzygl.2023.04.010

摘要/Abstract

摘要：

森林土壤碳库是全球生态系统碳库的重要组成部分,对人工林土壤有机碳空间分布及其影响因素进行研究,有助于了解人工林土壤碳汇动态变化规律。以河南省国有登封林场栓皮栎人工林为研究对象,采用系统抽样法在林地内设置186个样地(10 m×10 m),应用地统计学方法分析了土壤有机碳空间分布特征,构建并使用结构方程模型识别出土壤有机碳主要影响因素。结果表明,栓皮栎人工林土壤有机碳含量平均值为18.59 g/kg,其空间分布呈斑块状,高低值区分布较为分散;而平均土壤有机碳密度为2.47 kg/m²,其空间分布特征与土壤有机碳含量表现出高度一致性,但整体分布更均匀。相关分析结果表明,土壤有机碳含量和密度与胸径、树高等林分因子,与土壤全氮、土壤速效磷等土壤因子呈显著正相关关系(P<0.05);结构方程结果表明,土壤有机碳含量和土壤有机碳密度受土壤化学性质影响均最大,其中土壤全氮的贡献度最高。研究对于揭示人工林土壤碳汇形成机制,制定土壤碳汇调控技术措施有一定的参考价值。

关键词: 栓皮栎, 土壤有机碳, 空间分布, 结构方程模型

Abstract:

Forest soils play an important role in the global carbon cycle,and research on the spatial distribution of soil organic carbon and its influencing factors in plantation forests can help to understand the dynamic changes of soil carbon sink in plantation forests.Taking Quercus variabilis plantation in Dengfeng Forest farm,Henan Province as the research object,186 sample squares(10 m×10 m)were set in the forest by systematic sampling method,and the spatial distribution characteristics of soil organic carbon were analyzed by geostatistical methods,and the main factors influencing soil organic carbon were identified by constructing and using structural equation models.The results showed that the average soil organic carbon content of Quercus variabilis plantation was 18.59 g/kg,and its spatial distribution was patchy,with a more dispersed distribution of high and low value areas;while the average soil organic carbon density was 2.47 kg/m²,and its spatial distribution characteristics showed a high consistency with the soil organic carbon content,but the overall distribution was more uniform;the results of correlation analysis showed that soil organic carbon content and density were significantly and positively correlated(P<0.05)with stand factors such as diameter at breast height and tree height,and with soil factors such as total soil nitrogen and soil available phosphorus;the results of structural equation showed that both soil organic carbon content and soil organic carbon density were most influenced by soil chemical properties,with the highest contribution of total soil nitrogen.This study can provide some references for revealing the mechanism of soil carbon sink formation in artificial forests and developing technical measures for soil carbon sink regulation.

Key words: Quercus variabilis plantation, soil organic carbon, spatial distribution, structural equation model

中图分类号:

李思颖, 郭昊天, 陈晓蔚, 周梦丽, 靳姗姗, 闫东锋. 栓皮栎人工林土壤有机碳空间分布特征及其影响因素[J]. 林业资源管理, 2023, 0(4): 80-89.

LI Siying, GUO Haotian, CHEN Xiaowei, ZHOU Mengli, JIN Shanshan, YAN Dongfeng. Spatial Distribution Characteristics of Soil Organic Carbon in Quercus variabilis Plantation and Its Influencing Factors[J]. FOREST RESOURCES WANAGEMENT, 2023, 0(4): 80-89.

图/表 8

图1

表1

表2

图2

图3

表3

图4

表4

参考文献 36

[1]	Li Xiaojun, Li Yunfei, Xie Ting, et al. Recovery of soil carbon and nitrogen stocks following afforestation with xerophytic shrubs in the Tengger Desert,North China[J]. Catena, 2022, 214:106277. doi: 10.1016/j.catena.2022.106277
[2]	洪滔, 何晨阳, 黄贝佳, 等. 不同林龄千年桐人工林的碳含量和碳储量及碳库分配格局[J]. 植物资源与环境学报, 2021, 30(1):9-16.
[3]	陈家新, 杨红强. 全球森林及林产品碳科学研究进展与前瞻[J]. 南京林业大学学报:自然科学版, 2018, 42(4):1-8.
[4]	刘晓曼, 王超, 高吉喜, 等. 服务双碳目标的中国人工林生态系统碳增汇途径[J]. 生态学报, 2023(14):1-12.
[5]	Thom D, Keeton W S. Stand structure drives disparities in carbon storage in northern hardwood-conifer forests[J]. Forest Ecology and Management, 2019, 442:10-20. doi: 10.1016/j.foreco.2019.03.053
[6]	丁亚鹏, 张俊华, 刘玉寒, 等. 基于GWR模型的伊河流域土壤有机碳空间分布特征及影响因素分析[J]. 生态学报, 2021, 41(12):4876-4885.
[7]	连玉珍, 曹丽花, 刘合满, 等. 色季拉山西坡表层土壤有机碳的小尺度空间分布特征[J]. 北京林业大学学报, 2020, 42(9):70-79.
[8]	李龙, 姜丽娜, 白建华. 半干旱区土壤有机碳空间变异及其影响因素的多尺度相关分析[J]. 中国水土保持科学, 2018, 16(5):40-48.
[9]	Rosskopf N, Fell H, Zeitz J. Organic soils in Germany,their distribution and carbon stocks[J]. Catena, 2015, 133:157-170. doi: 10.1016/j.catena.2015.05.004
[10]	Zhou Yin, Chen Songchao, Zhu A-Xing, et al. Revealing the scale- and location-specific controlling factors of soil organic carbon in Tibet[J]. Geoderma, 2021, 382:114713. doi: 10.1016/j.geoderma.2020.114713
[11]	Sun Mengjiao, Hou Enqing, Wu Jiasen, et al. Spatial patterns and drivers of soil chemical properties in typical hickory plantations[J]. Forests, 2022, 13:457. doi: 10.3390/f13030457
[12]	秦艳培, 徐少君, 田耀武. 黄河流域河南段植被和土壤及其碳密度空间分异研究[J]. 生态环境学报, 2022, 31(9):1745-1753. doi: 10.16258/j.cnki.1674-5906.2022.09.004
[13]	McGrath D, Zhang Chaosheng. Spatial distribution of soil organic carbon concentrations in grassland of Ireland[J]. Applied Geochemistry, 2003, 18(10):1629-1639. doi: 10.1016/S0883-2927(03)00045-3
[14]	闫东锋, 张妍妍, 吕康婷, 等. 太行山南麓不同海拔梯度天然林优势树种生态位特征[J]. 生态环境学报, 2021, 30(8):1571-1580. doi: 10.16258/j.cnki.1674-5906.2021.08.003
[15]	王利彦, 周国娜, 朱新玉, 等. 凋落物对土壤有机碳与微生物功能多样性的影响[J]. 生态学报, 2021, 41(7):2709-2718.
[16]	王肖楠, 耿玉清, 余新晓. 栓皮栎林与油松林土壤有机碳及其组分的研究[J]. 土壤通报, 2012, 43(3):604-609.
[17]	和丽萍, 孟广涛, 李贵祥, 等. 金沙江头塘小流域人工林有机碳及其剖面分布特征[J]. 长江流域资源与环境, 2016, 25(3):476-485.
[18]	Yu Haiyan, Zha Tonggang, Zhang Xiaoxia, et al. Spatial distribution of soil organic carbon may be predominantly regulated by topography in a small revegetated watershed[J]. Catena, 2020, 188:104459. doi: 10.1016/j.catena.2020.104459
[19]	曹琪琪, 肖辉杰, 刘涛, 等. 乌兰布和沙漠东北部不同耕作年限农田土壤有机碳密度及其影响因素[J]. 应用生态学报, 2022, 33(10):2628-2634. doi: 10.13287/j.1001-9332.202210.017
[20]	Matheron G. Principles of geostatistics[J]. Economic Geology, 1963, 58(8):1246-1266. doi: 10.2113/gsecongeo.58.8.1246
[21]	Rossi R E, Mulla D J, Journel A G, et al. Geostatistical tools for modeling and interpreting ecological spatial dependence[J]. Ecological Monographs, 1992, 62(2):277-314. doi: 10.2307/2937096
[22]	管孝艳, 王少丽, 高占义, 等. 盐渍化灌区土壤盐分的时空变异特征及其与地下水埋深的关系[J]. 生态学报, 2012, 32(4):198-206.
[23]	苏永中, 赵哈林. 土壤有机碳储量、影响因素及其环境效应的研究进展[J]. 中国沙漠, 2002(3):19-27.
[24]	Ruiz-C M, Bienes R, Eldridge D J, et al. Vegetation cover reduces erosion and enhances soil organic carbon in a vineyard in the central Spain[J]. Catena, 2013, 104:153-160. doi: 10.1016/j.catena.2012.11.007
[25]	李龙, 秦富仓, 姜丽娜, 等. 赤峰市敖汉旗土壤有机碳含量的垂直分布及其影响因素[J]. 生态学报, 2019, 39(1):345-354.
[26]	王经波, 郑利林, 郭宇菲, 等. 鄱阳湖湿地土壤有机碳空间分布及其影响因素[J]. 长江流域资源与环境, 2022, 31(4):915-926.
[27]	祖元刚, 李冉, 王文杰, 等. 我国东北土壤有机碳、无机碳含量与土壤理化性质的相关性[J]. 生态学报, 2011, 31(18):5207-5216.
[28]	邱巡巡, 曹广超, 曹生奎, 等. 祁连山南坡农田土壤碳氮含量垂直分布特征及其影响因素[J]. 水土保持通报, 2022, 42(3):366-372.
[29]	Janssens I A, Dieleman W, Luyssaert S, et al. Reduction of forest soil respiration in response to nitrogen deposition[J]. Nature Geoscience, 2010, 3(5):315-322. doi: 10.1038/ngeo844
[30]	Treseder K K. Nitrogen additions and microbial biomass:A meta-analysis of ecosystem studies[J]. Ecology Letters, 2008, 11(10):1111-1120. doi: 10.1111/ele.2008.11.issue-10
[31]	祁金虎. 辽东山区天然次生栎林土壤有机碳含量及其与理化性质的关系[J]. 水土保持学报, 2017, 31(4):135-140.
[32]	向婷, 艾宁, 刘广全, 等. 陕北黄土区典型草地土壤水分特征及其对季节的响应[J]. 林业资源管理, 2022(2):149-156.
[33]	王越, 栾亚宁, 王丹, 等. 油松林土壤有机碳储量变化及其影响因素[J]. 浙江农林大学学报, 2021, 38(5):1023-1032.
[34]	张维俊, 李双异, 徐英德, 等. 土壤孔隙结构与土壤微环境和有机碳周转关系的研究进展[J]. 水土保持学报, 2019, 33(4):1-9.
[35]	Tian Yuhong, Liu Fenghua, Wang Tiantian. Spatial distribution of surface soil organic carbon density and related factors along an urbanization gradient in Beijing[J]. Journal of Resources and Ecology, 2020, 11(5):508-515. doi: 10.5814/j.issn.1674-764x.2020.05.008
[36]	柴华, 何念鹏. 中国土壤容重特征及其对区域碳贮量估算的意义[J]. 生态学报, 2016, 36(13):3903-3910.

变量	均值±标准差	最小值	最大值	变异系数/%
海拔/m	623±47.00	518	755	7.54
坡度/(°)	9±10.39	1	46	113.04
郁闭度	0.5±0.11	0.2	0.7	19.81
林龄/a	32±6.46	15	50	20.25
林分密度/(株/hm²)	1 430±275.95	600	2 600	19.30
平均胸径/cm	11.7±2.51	6.4	18.62	21.53
平均树高/m	8.1±2.35	4	15	29.02

因子	平均值±标准差	最小值	最大值	变异系数/%	偏度	峰度
SOC/(g/kg)	18.59±8.19	6.89	78.86	44.05	2.84	16.19
SOCD/(kg/m²)	2.47±1.17	1.09	12.15	47.48	24.96	3.65
STN/(g/kg)	2.15±0.96	0.82	9.34	44.69	24.30	17.52
SAN/(mg/kg)	56.22±51.44	11.78	460.83	90.97	4.46	27.57
SAP/(mg//kg)	1.11±0.11	0.92	1.64	9.70	1.26	3.07
SAK/(mg/kg)	123.33±59.42	50.17	466.57	48.18	2.20	6.78
SWC/%	11.3±2.90	4.90	17.40	25.64	0.01	-0.61
BD/(g/cm³)	1.34±0.17	0.83	1.81	13.06	0.35	0.18
Sa/a	32±6.46	15.00	50.00	20.25	-0.10	-0.23
D/cm	11.66±2.51	6.40	18.62	21.53	0.34	-0.21
H/m	8.11±2.35	4.00	15.00	29.02	0.39	-0.35
SD/(株/hm²)	1 430±275.95	600.00	2 600.00	19.30	0.46	2.98
TDR/%	4.21±1.73	1.14	14.82	41.14	1.74	7.68
TDF/%	4.46±1.53	1.41	15.35	34.21	2.40	13.89
TT/%	8.67±3.11	3.02	30.18	35.85	2.25	12.29

指标	SOC	SOCD	STN	SAN	SAP	SAK	SWC	BD	Sa	D	H	SD	TDR	TDF
SOCD	0.932^**
STN	0.936^**	0.842^**
SAN	0.556^**	0.485^**	0.616^**
SAP	0.226^**	0.144^*	0.283^**	0.199^**
SAK	0.358^**	0.320^**	0.351^**	0.291^**	0.088
SWC	0.216^**	0.048	0.286^**	0.335^**	0.133	0.129
BD	-0.164^*	0.154^*	-0.239^**	-0.265^**	-0.178^*	-0.131	-0.545^**
Sa	0.307^**	0.231^**	0.312^**	0.340^**	0.125	0.111	0.269^**	-0.226^**
D	0.325^**	0.269^**	0.308^**	0.327^**	0.098	0.241^**	0.198^**	-0.158^*	0.695^**
H	0.339^**	0.264^**	0.333^**	0.332^**	0.125	0.254^**	0.219^**	-0.202^**	0.766^**	0.698^**
SD	-0.016	-0.001	0.022	0.020	-0.011	-0.039	0.079	0.054	0.123	0.099	0.060
TDR	-0.055	-0.07	-0.064	-0.066	-0.006	-0.237^**	0.044	-0.108	-0.097	-0.091	-0.151^*	-0.079
TDF	-0.121	-0.115	-0.133	-0.092	-0.01	-0.261^**	-0.052	-0.017	-0.14	-0.092	-0.169^*	-0.113	0.732^**
TT	-0.085	-0.093	-0.095	-0.068	-0.011	-0.251^**	0.012	-0.087	-0.133	-0.104	-0.174^*	-0.102	0.945^**	0.905^**

影响因子	SOC	SOCD
STN	0.97^**	0.94^**
SAK	0.37^**	0.37^**
SAP	0.21^*	0.19^*
SWC	-0.07^**
BD	0.05^**	0.34^**