环境科学  2024, Vol. 45 Issue (3): 1674-1683   PDF    
黄河下游典型湿地土壤养分及其生态化学计量特征
王传盈1,2,3, 王凯月1,2,3, 王浩然1,2,3, 张梦迪1,2,3, 周云凯1,2,3     
1. 河南大学地理与环境学院, 开封 475004;
2. 河南大学黄河中下游数字地理技术教育部重点实验室, 开封 475004;
3. 河南省地球系统观测与模拟重点实验室, 开封 475004
摘要: 土壤中碳氮磷钾是植物生长发育所需的必要养分元素, 其含量和生态化学计量比能够反映土壤质量与养分限制状况.选取黄河下游河南段的花园口黄河浮桥湿地为研究对象, 运用方差分析、冗余分析和线性回归拟合等研究方法, 分析黄河下游典型湿地土壤有机碳(SOC)、全氮(TN)、全磷(TP)、全钾(TK)、碱解氮(AN)、速效磷(AP)、速效钾(AK)及其生态化学计量特征, 明确土壤养分的限制性元素, 揭示影响湿地土壤养分及其生态化学计量比的关键理化因子.结果表明, 湿地土壤ω(SOC)、ω(TN)、ω(TP)、ω(TK)、ω(AN)、ω(AP)和ω(AK)的均值分别为:5.46 g·kg-1、0.60 g·kg-1、0.28 g·kg-1、17.06 g·kg-1、13.75 mg·kg-1、6.54 mg·kg-1和158.56 mg·kg-1, 其空间分布呈现出从河岸向滩地波动增加的变化趋势, 高植被盖度区普遍高于低植被盖度区, 且土壤SOC、TN、TP和TK之间具有较强的相关性;土壤C/P、C/K、N/P和N/K呈现出和土壤养分一致的变化趋势, 而C/N则相反;土壤SOC、TN、AN、N/P和N/K变异系数超过50.00%, 其空间差异性较强;湿地土壤C/N均值为11.882, 接近中国土壤均值, 而C/P和N/P的均值分别为49.119和4.516, 两者均低于中国土壤平均水平, 并且土壤N/P远小于14, 研究区土壤受N限制;土壤黏粒占比和电导率共同解释了土壤养分变异与其生态化学计量比变异的61.4%和43.9%, 是影响花园口黄河浮桥湿地土壤养分和生态化学计量比的主要理化因子.研究结果有助于加深对黄河下游湿地土壤养分及其影响因素的认知, 为黄河下游湿地生态修复与管理提供重要科学依据.
关键词: 土壤养分      生态化学计量比      黄河下游湿地      限制性元素      土壤理化性质     
Nutrients and Ecological Stoichiometry Characteristics of Typical Wetland Soils in the Lower Yellow River
WANG Chuan-ying1,2,3 , WANG Kai-yue1,2,3 , WANG Hao-ran1,2,3 , ZHANG Meng-di1,2,3 , ZHOU Yun-kai1,2,3     
1. College of Geography and Environmental Science, Henan University, Kaifeng 475004, China;
2. Key Laboratory of Geospatial Technology for the Middle and Lower Yellow River Regions, Ministry of Education, Henan University, Kaifeng 475004, China;
3. Henan Key Laboratory of Earth System Observation and Modeling, Henan University, Kaifeng 475004, China
Abstract: Carbon, nitrogen, phosphorus, and potassium in the soil are the necessary nutrient elements for plant growth, and their contents and ecological stoichiometry can reflect the status of soil quality and nutrient limitation. The Huayuankou Yellow River Floating Bridge Wetland in the lower Yellow River was selected as the research object. The methods of ANOVA, redundancy analysis, and linear regression fitting were used to study the contents of organic carbon (SOC), total nitrogen (TN), total phosphorus (TP), total potassium (TK), alkaline nitrogen (AN), available phosphorus (AP), available potassium (AK), and their ecological stoichiometric ratios as well as the limiting elements of soil nutrients, and the key physicochemical properties that affect soil nutrients and their ecological stoichiometry in the wetland were revealed. The results showed that the mean values of ω(SOC), ω(TN), ω(TP), ω(TK), ω(AN), ω(AP), and ω(AK) in wetland soil were 5.46 g·kg-1, 0.60 g·kg-1, 0.28 g·kg-1, 17.06 g·kg-1, 13.75 mg·kg-1, 6.54 mg·kg-1, and 158.56 mg·kg-1, respectively, which showed an increasing trend from the river bank to the shoaly land and were generally higher at the high vegetation coverage areas than at the low vegetation coverage areas. There were significant correlations among SOC, TN, TP, and TK. Soil C/P, C/K, N/P, and N/K showed a consistent trend with soil nutrients, whereas C/N showed the opposite. The coefficients of variation of SOC, TN, AN, N/P, and N/K in the soil exceeded 50.00%, with significant spatial differences. The average value of C/N in wetland soil was 11.882, which was close to the average level of soils in China, whereas the average values of C/P and N/P were 49.119 and 4.516, respectively, both of which were lower than the average level of soils in China, and the N/P of soil was far less than 14, which indicated that N was limited in the soil. The proportion of clay and electrical conductivity combined to explain 61.4% and 43.9% of the variation in the soil nutrients and their ecological stoichiometry, respectively, which were the dominant soil physicochemical properties affecting the soil nutrients and their ecological stoichiometry of Huayuankou Yellow River Floating Bridge Wetland. The research results are helpful to improve our knowledge of nutrients and their influencing factors in the wetland soil of the lower Yellow River and provide an important scientific basis for the ecological restoration and management of the wetland in the lower Yellow River.
Key words: soil nutrient      ecological stoichiometry      lower Yellow River wetland      limiting element      soil physicochemical properties     

湿地是陆地生态系统与水生生态系统相互作用、相互影响的过渡区域[1], 是地球表面重要的生态系统.湿地独特的植被、水文和土壤等生态特征使得地表的各种元素在湿地汇集, 并且在复杂的生态过程驱动下, 养分元素在湿地土壤与植物之间进行循环、转化[2].碳(C)、氮(N)、磷(P)和钾(K)是湿地土壤养分的主要组成部分, 同时也是植物生长所需的必要常量元素, 其含量能够作为衡量土壤肥力水平的重要指标[3].土壤养分的供给能力决定着植物与微生物的生长发育, 进而影响生态系统的功能[4].

自1958年Redfield比率提出后, 生态化学计量学便逐渐成为新的研究领域并得到了快速发展[2, 5], 其为研究多种元素间相互关系提供了新的途径, 并成为分析土壤养分循环和确定土壤养分限制性元素的重要方法[6 ~ 9].从已有的研究进展来看, 生态化学计量研究区域更多集中在陆地生态系统和水生生态系统[9 ~ 13], 而对于湿地生态系统关注较少.然而湿地生态系统具有过渡性的特点, 其养分含量及生态化学计量比有别于上述两种生态系统, 且其影响因素更为复杂多样, 加强湿地生态系统养分和生态化学计量特征的研究可进一步丰富和完善营养元素生物地球化学循环理论与方法[2].

黄河作为世界含沙量最大的河流, 其所携带的泥沙在流动过程中逐渐沉积并在中下游区域形成了众多湿地和滩涂.然而由于全球气候变化以及人类活动的强烈影响, 黄河部分区段的湿地出现一定程度的退化, 为了解湿地土壤现状, 诸多学者对黄河湿地土壤养分及其化学计量特征进行了研究.屈凡柱等[5]与Qu等[14]研究了黄河三角洲不同芦苇湿地土壤的碳氮磷化学计量特征, 发现土壤C/N高于全国及全球平均水平而N/P则低于平均水平, 湿地土壤整体受到N限制.刘展航等[3]研究互花米草不同入侵年限下的湿地土壤化学计量特征发现, 黄河河口湿地土壤C/P和N/P均受土壤电导率、容重和含水量的影响.Li等[15]研究发现土壤含水量、盐分和pH是影响黄河三角洲湿地土壤碳氮磷及其化学计量比的重要因素.

目前, 围绕黄河三角洲或河口湿地土壤养分已开展了部分研究, 但对于黄河中下游湿地土壤养分缺乏较为深入的认知.因此, 本文选取黄河下游河南段的花园口黄河浮桥湿地为研究对象, 分析湿地土壤养分及其生态化学计量比的分布特征, 确定黄河下游河南段湿地土壤的限制性元素, 通过探究土壤理化性质和土壤养分及其生态化学计量比的关系, 揭示影响土壤养分及其生态化学计量比的关键理化因子, 以期为黄河湿地生态修复提供重要的科学依据.

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

黄河下游湿地河南段属于暖温带大陆性季风气候区, 夏季炎热多雨, 冬季寒冷干燥;年均温处于12~15℃之间, 年均降水量为500~900 mm, 但降水年际变化较大, 时空分布不均匀[16].本研究所选择的湿地位于河南省郑州市惠济区的花园口黄河浮桥周边区域(34°54′28.49″ ~ 34°54′30.74″N, 113°42′30.43″ ~113°42′43.75″E), 该湿地受人类活动影响相对较小, 自然条件保护较好, 植被覆盖程度较高, 优势植物主要有水稗草、芦苇和落叶阔叶林.

1.2 样品采集与测定

2022年6月在花园口黄河浮桥湿地进行样品采集.由河岸向外缘滩地方向依次设置5条平行于河岸的样带, 样带宽1.5 m, 各样带彼此间隔5~8 m, 且5条样带分别设置在不同的植被覆盖区.其中, 样带1位于河岸边缘, 毗邻黄河, 无植被覆盖;样带2位于水稗草群落, 植物分布较为稀疏;样带3位于芦苇群落, 植被盖度较高;样带4无明显优势植物, 植被盖度整体偏低;样带5位于落叶阔叶林, 距黄河最远, 林下存在草本植物, 植被盖度较高.在每条样带中分别设置3个1 m×1 m的样方, 每个样方沿对角线随机选择3个采样点, 除去表层凋落物后采集0~10 cm的土壤并进行混合, 将混合后的土壤样品装入编号自封袋中带回实验室进行化验分析.

采集的土壤样品经自然风干后去除杂物并研磨、过筛.土壤含水率采用烘干法测定(未风干土样);土壤pH利用pH仪测定;土壤电导率利用电导率仪测定;土壤颗粒占比利用Malvern激光粒度仪(Mastersizer 3000)进行检测;土壤有机碳(SOC)含量采用K2Cr2O7氧化法测定;土壤全氮(TN)含量采用凯氏定氮法测定;土壤碱解氮(AN)含量采用碱解扩散法测定;土壤全磷(TP)含量采用钼锑抗比色法测定;土壤速效磷(AP)含量采用NaHCO3浸提-钼锑抗比色法测定;土壤全钾(TK)含量采用NaOH熔融-火焰光度法测定;土壤速效钾(AK)采用CH3COONH4浸提-火焰光度法测定.

1.3 统计分析

运用Excel软件对实验数据进行处理和分析, 并计算湿地土壤的养分含量、生态化学计量比和变异系数(CV), CV用于表示土壤养分及其生态化学计量比的空间变异性, 其公式为:CV =标准偏差/均值×100%, 本研究中的化学计量比均为土壤SOC、TN、TP和TK的量比.运用SPSS 26软件对实验数据进行单因素方差分析(ANOVA)和多重比较, 对比湿地不同样带土壤养分及其生态化学计量比的差异显著性, 使用Levene's test对数据方差齐次性进行检验, 当方差齐次时, 使用最小显著差异法(LSD)进行多重比较;当方差不齐时, 使用Games-Howell法进行多重比较.运用Origin 2018软件作图并对实验数据进行线性回归拟合分析, 研究湿地土壤SOC、TN、TP和TK的相互关系.运用Canoco 5软件进行冗余分析(RDA), 研究土壤理化性质对土壤养分及其生态化学计量比的影响.

2 结果与分析 2.1 土壤养分及其生态化学计量特征 2.1.1 土壤养分特征

花园口黄河浮桥湿地土壤养分特征如表 1图 1所示.研究区土壤ω(SOC)均值为5.46 g·kg-1, 不同植被盖度的样带土壤中SOC差异显著(P < 0.05), 位于落叶阔叶林的样带5土壤SOC显著高于其他样带.湿地土壤ω(TN)均值为0.60 g·kg-1, 与SOC相似, 土壤TN最高值出现在样带5, 且不同样带的土壤TN差异显著(P < 0.05).研究区土壤SOC和TN沿远离河岸的方向呈现出波动增加.湿地土壤ω(TP)均值为0.28 g·kg-1, 不同植被盖度的样带土壤TP差异较小但仍较显著(P < 0.05), 高植被盖度的样带3和样带5土壤TP相近, 且显著高于植被盖度较低的土壤.湿地土壤ω(TK)均值为17.06 g·kg-1, 其最高值出现在植被盖度较高的样带3.湿地土壤的TP和TK在不同覆盖度的样带中呈现波动变化但总体较为平稳.研究区土壤ω(AN)和ω(AP)的均值分别为13.75 mg·kg-1和6.54 mg·kg-1, AN和AP在不同覆盖度的样带间呈波动上升的趋势, 且两者的最高值均出现在高植被盖度的样带5, 不同植被盖度的土壤AN和AP差异较显著(P < 0.05).土壤ω(AK)均值为158.56 mg·kg-1, 各样带间差异显著(P < 0.05), 高植被盖度的样带3土壤AK最高.从土壤养分的整体分布来看, 高植被盖度土壤的养分普遍高于低植被盖度土壤, 并且土壤养分基本呈现出由河岸向滩地方向波动增加的变化趋势.

表 1 花园口黄河浮桥湿地土壤养分特征 Table 1 Characteristics of soil nutrients in Huayuankou Yellow River Floating Bridge Wetland

不同小写字母表示土壤养分在不同样带下差异性达到显著(P < 0.05) 图 1 花园口黄河浮桥湿地不同样带土壤养分特征 Fig. 1 Characteristics of soil nutrients in different transects in Huayuankou Yellow River Floating Bridge Wetland

从湿地整体来看, 土壤TN和AN的变异系数较高, 分别为70.28%和63.52%, 变异系数高于其他土壤养分, 空间变异性较强;土壤SOC、AK和AP的变异系数超过40.00%, 属中等程度变异;土壤TP和TK的变异系数最低, 均低于20.00%, 空间变异性较弱.从变异系数来看, 土壤TN和AN在花园口黄河浮桥湿地中分布的空间差异最大, 而TP和TK在湿地中分布较为均匀但有效性养分空间差异较大.

线性回归拟合结果表明(图 2), 花园口黄河浮桥湿地土壤的SOC、TN、TP和TK之间呈现显著的相关性(P < 0.05), 养分存在强烈的相互作用关系.湿地土壤SOC与TN含量的拟合程度优于SOC与TP和TK的拟合, 而TN与TK的拟合程度优于TN与TP的拟合.对比线性回归拟合结果的R2可以发现, 土壤SOC和TN的相关性最强且极显著(P < 0.01), 土壤TP和TK相关性较弱, 但相关性仍较为显著(P < 0.05).

图 2 湿地土壤SOC、TNTPTK的关系 Fig. 2 Relationships of SOC, TN, TP, and TK in wetland soil

2.1.2 土壤生态化学计量比特征

花园口黄河浮桥湿地土壤生态化学计量特征如表 2图 3所示.研究区土壤的C/N介于8.945~14.045, 均值为11.882, 从河岸到滩地方向湿地土壤C/N呈波动下降趋势, 最低值出现在高植被盖度的样带5.湿地土壤的C/P、C/K、N/P和N/K分别介于19.513~75.642、0.441~1.628、1.412~8.484和0.032~0.182, 其均值分别为49.119、1.006、4.516和0.093, 4个比值在不同植被盖度的样带间差异较为显著(P < 0.05), 且从河岸到滩地方向均呈现出波动上升的变化趋势, 最高值均出现在高植被盖度的样带5.研究区土壤P/K介于0.016~0.023, 均值为0.021, 除样带4显著低于其他样带外(P < 0.05), 其余样带间差异不显著.从生态化学计量比的整体分布来看, 除C/N和P/K外, 研究区土壤生态化学计量比均呈现由低植被盖度区向高植被盖度区波动增加的特征.

表 2 花园口黄河浮桥湿地土壤生态化学计量比特征 Table 2 Characteristics of soil ecological stoichiometric ratio in Huayuankou Yellow River Floating Bridge Wetland

不同小写字母表示土壤生态化学计量比在不同样带下差异性显著(P < 0.05) 图 3 花园口黄河浮桥湿地不同样带土壤生态化学计量比特征 Fig. 3 Characteristics of soil ecological stoichiometric ratio in different transects in Huayuankou Yellow River Floating Bridge Wetland

从变异系数来看, 湿地土壤的N/P和N/K变异系数分别为59.36%和61.31%, 属于中等变异程度, 相较于其他化学计量比, N/P和N/K的空间变异性最强.土壤C/P和C/K的变异系数分别为43.68%和43.26%, 空间变异性属中等程度, 在整个湿地中变化较大.土壤C/N和P/K的变异系数最低, 均低于20.00%, 其在空间内的分布较为均匀.

2.2 土壤理化性质对土壤养分及其生态化学计量比的影响 2.2.1 土壤理化性质对土壤养分的影响

以土壤养分为响应变量, 土壤理化性质为解释变量进行RDA.在分析前, 首先对土壤养分和土壤理化性质进行去趋势对应分析(DCA)检验, 检验结果表明排序轴最长梯度的长度为0.25, 数值远小于3.00, 表明土壤养分和土壤理化性质适合进行RDA.

根据RDA结果, 轴1和轴2的特征值分别为0.627 7和0.132 3, 前两轴的累计解释率为76.00%.土壤养分和土壤理化性质间的相关性在第一轴达到0.858 4, 在第二轴达到0.995 8, 相关性较高, 分析结果具有可信性.由表 3可知, 各土壤理化性质对土壤养分影响的重要性表现为:黏粒占比 > 电导率 > 粉粒占比 > 含水率 > pH > 砂粒占比, 其中黏粒占比、电导率和粉粒占比对土壤养分具有显著影响(P < 0.05), 其对土壤养分的贡献率分别为46.2%、34.1%和14.5%, 且分别解释了土壤养分35.3%、26.1%和11.1%的变异, 而其余土壤理化因子对土壤养分的影响并未达到显著水平.依据RDA排序结果可知(图 4), 土壤黏粒占比与土壤所有养分均呈极显著正相关, 土壤电导率和粉粒占比与TP、TK、AN、AP和AK呈显著正相关关系.上述分析表明, 土壤黏粒占比、电导率和粉粒占比是影响花园口黄河浮桥湿地土壤养分的主要土壤理化因子.

表 3 土壤理化性质对土壤养分重要性排序 Table 3 Importance ranking table of soil physicochemical properties to soil nutrients

红色空心箭头表示土壤理化因子, 黑色实心箭头表示土壤养分 图 4 土壤理化性质与土壤养分的RDA排序结果 Fig. 4 RDA ranking results of soil physicochemical properties and soil nutrients

2.2.2 土壤理化性质对土壤生态化学计量比的影响

以土壤生态化学计量比为响应变量, 土壤理化性质为解释变量进行RDA.在分析前, 首先对土壤生态化学计量比和土壤理化性质进行DCA检验, 检验结果表明排序轴最长梯度的长度为0.39, 数值远小于3.00, 表明土壤生态化学计量比和土壤理化性质适合进行RDA.

根据RDA结果, 轴1和轴2的特征值分别为0.651 0和0.204 1, 前两轴的累计解释率为85.51%.土壤理化性质和土壤生态化学计量比的相关性在第一轴达到0.933 3, 在第二轴达到0.756 8, 相关性较高, 分析结果具有可信性.根据表 4结果, 各土壤理化性质对土壤生态化学计量比影响的重要性表现为:黏粒占比 > 粉粒占比 > pH > 电导率 > 含水率 > 砂粒占比, 其中黏粒占比和电导率对土壤生态化学计量比具有显著影响(P < 0.05), 其对土壤生态化学计量比的贡献率分别为39.5%和17.5%, 且分别解释了生态化学计量比30.4%和13.5%的变异, 而其他土壤理化因子对土壤生态化学计量特征的影响均未达到显著水平.依据RDA排序结果可以发现(图 5), 黏粒占比与C/P、C/K、N/P和N/K呈极显著正相关关系, 与C/N和P/K呈负相关关系;电导率与C/N和P/K呈正相关关系, 与C/P、C/K、N/P和N/K呈负相关关系.上述分析表明, 土壤黏粒占比和电导率是影响花园口黄河浮桥湿地土壤生态化学计量比的主要土壤理化因子.

表 4 土壤理化性质对土壤生态化学计量比重要性排序 Table 4 Importance ranking table of soil physicochemical properties to soil ecological stoichiometric ratio

红色空心箭头表示土壤理化因子, 黑色实心箭头表示土壤生态化学计量比 图 5 土壤理化性质与土壤生态化学计量比的RDA排序结果 Fig. 5 RDA ranking results of soil physicochemical properties and soil ecological stoichiometric ratio

3 讨论 3.1 湿地土壤养分及其生态化学计量特征

土壤作为植物生长必不可少的重要环境因子, 土壤养分是植物养分的重要来源[17], 植物在受土壤养分影响的同时也会影响土壤养分的含量和分布.一方面, 不同植物养分分配机制的差异会使植物选择性地吸收养分, 影响土壤中各养分的平衡[18];另一方面, 植物通过光合作用从大气中吸收太阳能和二氧化碳, 并以植物残体和凋落物的形式向土壤提供大量的碳、氮等养分[2, 18, 19], 改变土壤养分含量, 影响土壤养分循环[20].植被覆盖在很大程度上影响土壤养分含量[21], 具有高植被覆盖度的土壤, 其植物生物量较高, 根系及其凋落物归还到土壤中的数量多[22], 凋落物和腐殖质分解能够改变土壤微生物的群落结构和活动并促进土壤团聚体的形成, 而土壤团聚体能够通过多种途径保护有机质, 利于养分积累[23], 从而使得土壤养分含量相对较高.本研究中, 样带3中的芦苇生产力较高、生物量较大, 其凋落物就地分解后, 养分返还到土壤中, 土壤养分含量相对较高;落叶阔叶林的枯枝落叶较多, 且林下生存着草本植物, 这使得样带5中由植物返还到土壤的养分也较为丰富, 促进了土壤养分的积累.而其他样带的植被盖度相对较低甚至无植被覆盖, 植物生物量较小, 由植物返还到土壤中的养分远少于样带3和样带5, 不同植被覆盖的土壤养分含量存在显著差异, 这与之前的研究具有相似性[7, 24, 25].此外, 植被覆盖与土壤侵蚀密切相关, 高植被覆盖的土壤具有茂盛的地上植物群落、地表凋落物及复杂的地下植物根系, 能够一定程度上缓解水力侵蚀、风力侵蚀等引起的土壤养分迁移与流失, 有利于土壤养分的保持[26].相较于其他样带, 样带3和样带5的植被覆盖度较高, 可以一定程度上减轻因大气降水或丰水期河水淹没所引起的养分流失, 表层土壤因而保持了较高的养分.

植被覆盖度的差异同样会引起土壤生态化学计量比的变化, 植物通过养分吸收和凋落物分解输入等方式影响土壤养分分布, 植物凋落物的高输入能够为微生物提供更丰富的营养条件并提高微生物生物量, 同时, 植物的根系分泌物可以通过提高或减弱酶活性以调节微生物与土壤的相互作用[27], 进而影响微生物对土壤中有机物质的矿化过程, 引起土壤生态化学计量分布特征的变化[28].土壤C/N主要受植物凋落物数量的影响, 输入到土壤中的凋落物越多, 微生物分解的基质越充足, 有机质分解速度越快, 土壤C/N越低[29].本研究中, 从黄河河岸到滩地的植被盖度整体呈现上升趋势, 土壤中输入的凋落物数量及有机质分解速率不断增加, 因而C/N整体呈下降趋势.虽然土壤养分整体均呈现增加的趋势, 但由于土壤SOC和TN在样带间的变化程度高于TP和TK, 因而C/P、C/K、N/P和N/K呈现由河岸向滩地方向不断增加的特点.此外, 因植被覆盖差异而引起的土壤养分变化也是导致湿地土壤化学计量比变化的重要原因.

对比整个湿地土壤各养分的变异系数可以发现, 花园口黄河浮桥湿地土壤中SOC、TN和AN的空间变异性较大, 而TP和TK的空间变异性较小, 这主要是由于不同元素的来源差异而引起的.有研究表明, 土壤P和K主要受到土壤成土母质的影响, 其来源主要是岩石的风化分解[2, 30], 花园口黄河浮桥湿地土壤是受到黄河多次泛滥冲积而形成的黄河沉积土质[31], 在整个湿地中土壤的成土母质无明显差异, 因而湿地土壤TP和TK的变异性较小.土壤C和N除受到成土母质的影响外, 还受植物凋落物及残体分解的影响[32, 33], 花园口黄河浮桥湿地不同样带的植被类型及其覆盖程度差异较大, 植物对养分的吸收利用以及通过凋落物返还到土壤中的养分含量不同, 因而土壤SOC、TN和AN存在较大的空间变异性.土壤养分的空间差异会引起其生态化学计量比的空间变化, 由于土壤C和N的较大空间变异性以及P和K的相对稳定性, 花园口黄河浮桥湿地土壤C/P、C/K、N/P和N/K具有较大的空间变异性, 而P/K则变化较小.此外, 由于C和N是结构性元素, 其积累与消耗的比例较为固定, 因而土壤中的C/N也相对稳定[34].

土壤养分彼此之间密切联系和相互影响[2, 19, 33].线性回归拟合结果表明, 花园口黄河浮桥湿地土壤SOC、TN、TP和TK之间具有显著相关性, 其中土壤SOC与TN的相关性最强.湿地土壤C和N的强相关性可从以下两方面解释:一方面, SOC和TN是土壤有机质的重要组成部分, C和N的输入与元素比例相对稳定的植物体、微生物及两者进入土壤的残留物有关[35];另一方面, 土壤TN能够影响土壤微生物的活性, 促进或抑制微生物对有机质的分解, 进而影响土壤中SOC所占的质量分数, 因此土壤TN与SOC的质量分数具有较为密切的联系[36].此外湿地土壤的C/N空间变化较小, 且土壤SOC与TN在样带间的变化趋势几乎完全一致, 这进一步证明了花园口黄河浮桥湿地土壤C和N具有较强的相互作用, 这与黄河三角洲湿地的研究相一致[37].

3.2 黄河下游河南段湿地土壤限制性养分元素

土壤生态化学计量比作为表征土壤养分固定、维持与分解释放的重要指标, 其在土壤养分循环与限制作用的研究中被广泛应用[38], 通过对生态化学计量比进行分析能够深入理解土壤的限制性元素.C/N与土壤有机质分解速率呈反比关系, 能够反映有机质的分解速率, 同时也能够指示土壤N的淋溶风险, 有研究表明, 当C/N低于30时, 土壤中硝酸盐的淋溶风险高, N易流失[39].花园口黄河浮桥湿地土壤C/N为11.882, 接近中国土壤平均C/N(11.900)[40], 且其数值低于30, 表明研究区土壤有机质分解速度较快, 硝酸盐淋溶风险较高, N流失程度较强.C/P是分析土壤P有效性的重要指标, 能够反映土壤P的矿化能力, 土壤中C/P越低, 越有利于有机质分解过程中P的释放, 土壤P活性越强[2].花园口黄河浮桥湿地土壤的C/P(49.119)低于中国土壤平均C/P(61.000)[40], 表明研究区湿地土壤具有高P活性且有机磷净矿化.土壤N/P通常被作为氮饱和的重要指标, 可以用于分析土壤养分的限制元素, 当土壤N/P低于14时, 认为土壤受到N的限制;当土壤N/P高于16时, 认为土壤受到P的限制[2, 41].在本研究中, 花园口黄河浮桥湿地土壤N/P(4.516)低于14, 并且与中国土壤N/P均值(5.200)[40]相比偏低, 土壤更容易受到N限制.总体来看, 相对稳定的C/N与较低的N/P和C/P表明花园口黄河浮桥湿地土壤相对富含C和P, 而N相对匮乏, 湿地土壤整体受到N的限制, 这与已有研究表明黄河三角洲湿地生态系统主要受到N限制的结论具有一致性[5, 14, 42, 43].

除与硝酸盐淋溶风险较高有关外, 湿地土壤受N限制还可能与土壤pH有关.研究表明, 在pH较高时, 土壤中异化硝酸盐还原为铵(dissimilatory nitrate reduction to ammonium, DNRA)过程和反硝化过程会同时发生, 且土壤中氮气的通量会随着pH的升高而增加[44].而土壤中的硝态氮会通过DNRA过程产生的部分铵态氮参与到氨挥发过程, 同时随pH增大, 铵态氮向氨转化速率加快[45], 导致土壤氨挥发程度较强.经实验分析, 花园口黄河浮桥湿地土壤pH均值为9.60, 这种环境下, 土壤可能会发生反硝化过程及DNRA过程而引起土壤N的损失, 在一定程度上导致土壤N的相对匮乏.

3.3 湿地土壤理化性质对土壤养分及其生态化学计量的影响

土壤养分及其生态化学计量比与土壤理化性质之间具有十分密切的联系.土壤质地能够影响土壤中养分的迁移与积累, 土壤颗粒物中黏粒和粉粒的直径较小, 具有较高的比表面积, 其持水能力较强, 对养分具有较高的吸附能力, 能够影响土壤养分的循环过程[46], 因此土壤细颗粒物含量会对土壤的养分产生影响.本研究中RDA结果显示, 湿地土壤养分均与黏粒和粉粒占比具有显著相关关系, 且两者对土壤养分的贡献率达到了60.5%, 这表明土壤黏粒和粉粒占比是影响花园口黄河浮桥湿地土壤养分的重要因素.土壤细颗粒物会通过影响土壤养分进而引起土壤养分生态化学计量比的变化[47].在花园口黄河浮桥湿地中, 土壤的黏粒占比与生态化学计量比的相关性达到极显著水平, 其贡献率达到了39.5%, 表明土壤黏粒占比是影响土壤生态化学计量比的重要土壤理化因子.

湿地土壤盐分与养分及其生态化学计量比息息相关[48], 土壤盐分通过影响微生物活性、植物生长等方式引起土壤微生物数量、活性及植物生物量的变化, 对植物凋落物的输入量和土壤中养分的积累与矿化分解过程产生影响, 进而导致土壤养分发生变化并影响养分生态化学计量比[49].已有研究表明, 盐分与土壤养分及其生态化学计量比间存在显著相关性, 刘文龙等[49]研究胶州湾芦苇潮滩湿地土壤养分生态计量特征的影响因素时发现, 土壤盐度与C/P和N/P呈显著负相关, 且相对于土壤pH与含水量, 土壤盐度对养分生态化学计量特征的影响程度更高.Yan等[50]通过对中国东南部红壤碳氮磷化学计量特征的研究中发现, 土壤电导率与SOC、TN和TP呈显著正相关, 与C/P和N/P呈负相关关系.在本研究中, 土壤电导率与养分达到极显著相关关系, 除SOC外, 其他土壤养分均与电导率呈正相关, 各土壤养分生态化学计量比与电导率也具有显著相关性, 这表明土壤盐分显著影响了花园口黄河浮桥湿地土壤的养分与其生态化学计量特征.

4 结论

(1) 花园口黄河浮桥湿地土壤养分ω(SOC)、ω(TN)、ω(TP)、ω(TK)、ω(AN)、ω(AP)和ω(AK)的均值分别为:5.46 g·kg-1、0.60 g·kg-1、0.28 g·kg-1、17.06 g·kg-1、13.75 mg·kg-1、6.54 mg·kg-1和158.56 mg·kg-1.研究区土壤养分由河岸向滩地方向呈现出波动增加的趋势, 普遍呈现出高植被盖度区高于低植被盖度区.土壤中SOC、TN和AN具有较大的空间变异性, 而TP和TK分布较为均匀, 且土壤SOC、TN、TP和TK之间具有较强的相关性.

(2) 与养分特征相似, 研究区土壤C/P、C/K、N/P和N/K也呈现出高植被盖度区高于低植被盖度区的特征, 且这些比值均呈现出由河岸向滩地方向波动增大的变化特点, 而土壤C/N则相反.研究区中土壤N/P和N/K具有较强的空间变异性, 其分布变化较大, 而C/N和P/K分布较为均匀.湿地土壤的C/N、C/P和N/P分别为11.882、49.119和4.516, 其中C/N接近中国土壤的平均水平, 而C/P和N/P均低于中国土壤平均水平, 表明了花园口黄河浮桥湿地土壤相对富含C和P, 而N匮乏, 研究区土壤整体受到N限制.

(3) 研究区湿地土壤的黏粒占比与电导率解释了土壤养分变异的61.4%和生态化学计量比变化的43.9%, 其总贡献率分别达到80.3%和57.0%, 两者是影响花园口黄河浮桥湿地土壤养分及其生态化学计量比的主要土壤理化因子.

参考文献
[1] 吕宪国. 湿地科学研究进展及研究方向[J]. 中国科学院院刊, 2002, 17(3): 170-172.
Lv X G. A review and prospect for wetland science[J]. Bulletin of Chinese Academy of Sciences, 2002, 17(3): 170-172. DOI:10.3969/j.issn.1000-3045.2002.03.004
[2] Wu Y H, Wu Z J, Jiang S M, et al. Elemental stoichiometry (C, N, P) of soil in the wetland critical zone of Dongting Lake, China: Understanding soil C, N and P status at greater depth[J]. Sustainability, 2022, 14(14). DOI:10.3390/su14148337
[3] 刘展航, 张树岩, 侯玉平, 等. 互花米草入侵对黄河口湿地土壤碳氮磷及其生态化学计量特征的影响[J]. 生态环境学报, 2022, 31(7): 1360-1369.
Liu Z H, Zhang S Y, Hou Y P, et al. Effects of Spartina alterniflora invasion on soil carbon, nitrogen, phosphorus and their ecostoichiometric characteristics in the Yellow River Estuary wetlands[J]. Ecology and Environmental Sciences, 2022, 31(7): 1360-1369.
[4] Chen Y, Li Y Q, Duan Y L, et al. Patterns and driving factors of soil ecological stoichiometry in typical ecologically fragile areas of China[J]. CATENA, 2022, 219. DOI:10.1016/j.catena.2022.106628
[5] 屈凡柱, 孟灵, 付战勇, 等. 不同生境条件下滨海芦苇湿地C, N, P化学计量特征[J]. 生态学报, 2018, 38(5): 1731-1738.
Qu F Z, Meng L, Fu Z Y, et al. The stoichiometry characterization of carbon, nitrogen and phosphorus in different reed-dominated coastal wetland habitats[J]. Acta Ecologica Sinica, 2018, 38(5): 1731-1738.
[6] Elser J J, Sterner R W, Gorokhova E, et al. Biological stoichiometry from genes to ecosystems[J]. Ecology Letters, 2000, 3(6): 540-550. DOI:10.1111/j.1461-0248.2000.00185.x
[7] 张萌, 沈雅飞, 陈天, 等. 宜阳县不同森林类型土壤化学计量特征[J]. 陆地生态系统与保护学报, 2022, 2(2): 1-8.
Zhang M, Shen Y F, Chen T, et al. Soil stoichiometric characteristics of different forest types in Yiyang County[J]. Terrestrial Ecosystem and Conservation, 2022, 2(2): 1-8.
[8] 李从娟, 雷加强, 徐新文, 等. 塔克拉玛干沙漠腹地人工植被及土壤CNP的化学计量特征[J]. 生态学报, 2013, 33(18): 5760-5767.
Li C J, Lei J Q, Xu X W, et al. The stoichiometric characteristics of C, N, P for artificial plants and soil in the hinterland of Taklimakan Desert[J]. Acta Ecologica Sinica, 2013, 33(18): 5760-5767.
[9] 高郯, 权红, 卢杰, 等. 藏东南高山松林表层土壤养分含量及其化学计量比特征[J]. 西北农林科技大学学报(自然科学版), 2021, 49(4): 62-70, 80.
Gao T, Quan H, Lu J, et al. Characteristics of topsoil nutrients and their stoichiometry in Pinus densata forest in Southeast Tibet[J]. Journal of Northwest A & F University (Natural Science Edition), 2021, 49(4): 62-70, 80.
[10] 叶玉适, 陈进, 许继军, 等. 长江源和怒江源区水体氮磷化学计量特征初探[J]. 长江科学院院报, 2019, 36(9): 6-11, 22.
Ye Y S, Chen J, Xu J J, et al. Stoichiometric characteristics of nitrogen and phosphorus in water body in the source regions of the Yangtze River and Nujiang River[J]. Journal of Yangtze River Scientific Research Institute, 2019, 36(9): 6-11, 22.
[11] Zhou J, Han X X, Brookes J D, et al. High probability of nitrogen and phosphorus co-limitation occurring in eutrophic lakes[J]. Environmental Pollution, 2022, 292. DOI:10.1016/j.envpol.2021.118276
[12] 赵雯, 黄来明. 高寒山区不同土地利用类型土壤养分化学计量特征及影响因素[J]. 生态学报, 2022, 42(11): 4415-4427.
Zhao W, Huang L M. Stoichiometric characteristics and influencing factors of soil nutrients under different land use types in an alpine mountain region[J]. Acta Ecologica Sinica, 2022, 42(11): 4415-4427.
[13] 南富森, 李宗省, 张小平, 等. 黄河北岸兰州段丘陵区土壤生态化学计量与空间变异[J]. 中国沙漠, 2022, 42(5): 167-176.
Nan F S, Li Z X, Zhang X P, et al. Spatial variation of ecological stoichiometry characteristics of hilly soil in Lanzhou section of Yellow River north shore[J]. Journal of Desert Research, 2022, 42(5): 167-176.
[14] Qu F Z, Yu J B, Du S Y, et al. Influences of anthropogenic cultivation on C, N and P stoichiometry of reed-dominated coastal wetlands in the Yellow River Delta[J]. Geoderma, 2014, 235-236: 227-232. DOI:10.1016/j.geoderma.2014.07.009
[15] Li D X, Li Y N, Xie Y L, et al. Effects of ecological restoration on soil biogenic elements and their ecological stoichiometry in the Yellow River Delta, China[J]. Frontiers in Marine Science, 2022, 9. DOI:10.3389/fmars.2022.993202
[16] 张会萍, 肖曼, 曾曼芸, 等. 河南段黄河湿地水文因子对植物多样性的影响[J]. 河南农业大学学报, 2022, 56(5): 822-829.
Zhang H P, Xiao M, Zeng M Y, et al. Effects of hydrological environment on plant diversity in the Yellow River Wetland of Henan Section[J]. Journal of Henan Agricultural University, 2022, 56(5): 822-829.
[17] 张雨鉴, 王克勤, 宋娅丽, 等. 滇中亚高山森林植物叶-凋落叶-土壤生态化学计量特征[J]. 生态学报, 2020, 40(21): 7648-7658.
Zhang Y J, Wang K Q, Song Y L, et al. C, N and P stoichiometric characteristics of leaf, litter and soil for subalpine forests in central Yunnan, China[J]. Acta Ecologica Sinica, 2020, 40(21): 7648-7658.
[18] Zhang Z S, Song X L, Lu X G, et al. Ecological stoichiometry of carbon, nitrogen, and phosphorus in estuarine wetland soils: influences of vegetation coverage, plant communities, geomorphology, and seawalls[J]. Journal of Soils and Sediments, 2013, 13(6): 1043-1051. DOI:10.1007/s11368-013-0693-3
[19] Wang L L, Zhang G H, Zhu P Z, et al. Soil C, N and P contents and their stoichiometry as affected by typical plant communities on steep gully slopes of the Loess Plateau, China[J]. CATENA, 2022, 208. DOI:10.1016/j.catena.2021.105740
[20] 张慧慧, 白云玉, 张英洁, 等. 长白山苔原带凋落物生态化学计量特征及其对模拟氮沉降的响应[J]. 生态学报, 2022, 42(21): 8795-8808.
Zhang H H, Bai Y Y, Zhang Y J, et al. Response of chemical composition and ecological stoichiometric characteristics of three types of litter to simulated nitrogen deposition in the Changbai Mountain tundra[J]. Acta Ecologica Sinica, 2022, 42(21): 8795-8808.
[21] Guo N, Degen A A, Deng B, et al. Changes in vegetation parameters and soil nutrients along degradation and recovery successions on alpine grasslands of the Tibetan plateau[J]. Agriculture, 2019, 284. DOI:10.1016/j.agee.2019.106593
[22] 安克俭, 魏霞, 贺燕, 等. 祁连山区4种高寒植被类型下土壤养分及含水率分布[J]. 水土保持通报, 2022, 42(3): 16-22.
An K J, Wei X, He Y, et al. Soil moisture content and nutrient distribution under four alpine vegetation types in Qilian Mountains[J]. Bulletin of Soil and Water Conservation, 2022, 42(3): 16-22.
[23] Li J W, Wang K B, Shangguan Z P, et al. Coupling and decoupling of soil carbon, nitrogen and phosphorus stocks following grazing exclusion in temperate grasslands[J]. CATENA, 2022, 211. DOI:10.1016/j.catena.2021.106003
[24] 汪宗飞, 郑粉莉. 黄土高原子午岭地区人工油松林碳氮磷生态化学计量特征[J]. 生态学报, 2018, 38(19): 6870-6880.
Wang Z F, Zheng F L. C, N, and P stoichiometric characteristics of Pinus tabulaeformis plantation in the Ziwuling Region of the Loess Plateau[J]. Acta Ecologica Sinica, 2018, 38(19): 6870-6880.
[25] 温晨, 杨智姣, 杨磊, 等. 半干旱黄土小流域不同植被类型植物与土壤生态化学计量特征[J]. 生态学报, 2021, 41(5): 1824-1834.
Wen C, Yang Z J, Yang L, et al. Ecological stoichiometry characteristics of plants and soil under different vegetation types in the semi-arid loess small watershed[J]. Acta Ecologica Sinica, 2021, 41(5): 1824-1834.
[26] Chang E H, Li P, Li Z B, et al. The impact of vegetation successional status on slope runoff erosion in the loess plateau of China[J]. Water, 2019, 11(12). DOI:10.3390/w11122614
[27] Bhattacharyya S S, Furtak K. Soil-plant-microbe interactions determine soil biological fertility by altering rhizospheric nutrient cycling and biocrust formation[J]. Sustainability, 2023, 15(1). DOI:10.3390/su15010625
[28] Bing H J, Wu Y H, Zhou J, et al. Stoichiometric variation of carbon, nitrogen, and phosphorus in soils and its implication for nutrient limitation in alpine ecosystem of Eastern Tibetan Plateau[J]. Journal of Soils and Sediments, 2016, 16(2): 405-416. DOI:10.1007/s11368-015-1200-9
[29] 张子琦, 焦菊英, 陈同德, 等. 拉萨河流域洪积扇不同植被类型土壤化学计量特征[J]. 生态学报, 2022, 42(16): 6801-6815.
Zhang Z Q, Jiao J Y, Chen T D, et al. Soil stoichiometry characteristics of different vegetation types in alluvial fans of the Lhasa River basin[J]. Acta Ecologica Sinica, 2022, 42(16): 6801-6815.
[30] Augusto L, Achat D L, Jonard M, et al. Soil parent material-A major driver of plant nutrient limitations in terrestrial ecosystems[J]. Global Change Biology, 2017, 23(9): 3808-3824. DOI:10.1111/gcb.13691
[31] 聂永博. 黄河下游滩涂湿地生态系统演变特征及服务价值评估[D]. 郑州: 华北水利水电大学, 2022.
Nie Y B. Evolution characteristics and service value evaluation of beach wetland ecosystem in the lower Yellow River[D]. Zhengzhou: North China University of Water Resources and Electric Power, 2022.
[32] 朱秋莲, 邢肖毅, 张宏, 等. 黄土丘陵沟壑区不同植被区土壤生态化学计量特征[J]. 生态学报, 2013, 33(15): 4674-4682.
Zhu Q L, Xing X Y, Zhang H, et al. Soil ecological stoichiometry under different vegetation area on loess hilly-gully region[J]. Acta Ecologica Sinica, 2013, 33(15): 4674-4682.
[33] 陶冶, 张元明, 周晓兵. 伊犁野果林浅层土壤养分生态化学计量特征及其影响因素[J]. 应用生态学报, 2016, 27(7): 2239-2248.
Tao Y, Zhang Y M, Zhou X B. Ecological stoichiometry of surface soil nutrient and its influencing factors in the wild fruit forest in Yili region, Xinjiang, China[J]. Chinese Journal of Applied Ecology, 2016, 27(7): 2239-2248.
[34] Qi D D, Feng F J, Lu C, et al. C: N: P stoichiometry of different soil components after the transition of temperate primary coniferous and broad-leaved mixed forests to secondary forests[J]. Soil and Tillage Research, 2022, 216. DOI:10.1016/j.still.2021.105260
[35] Wang Y, Ren Z, Ma P P, et al. Effects of grassland degradation on ecological stoichiometry of soil ecosystems on the Qinghai-Tibet Plateau[J]. Science of the Total Environment, 2020, 722. DOI:10.1016/j.scitotenv.2020.137910
[36] 曹小玉, 李际平, 闫文德. 不同龄组杉木林土壤有机碳与氮磷钾分布特征及耦合关系[J]. 土壤通报, 2014, 45(5): 1137-1143.
Cao X Y, Li J P, Yan W D. Variation of contents of soil organic Carbon, N, P and K and their coupling relationships in different age groups of Chinese Fir Plantations[J]. Chinese Journal of Soil Science, 2014, 45(5): 1137-1143.
[37] 孙德斌, 栗云召, 于君宝, 等. 黄河三角洲湿地不同植被类型下土壤营养元素空间分布及其生态化学计量学特征[J]. 环境科学, 2022, 43(6): 3241-3252.
Sun D B, Li Y Z, Yu J B, et al. Spatial distribution and eco-stoichiometric characteristics of soil nutrient elements under different vegetation types in the Yellow River Delta wetland[J]. Environmental Science, 2022, 43(6): 3241-3252.
[38] 王绍强, 于贵瑞. 生态系统碳氮磷元素的生态化学计量学特征[J]. 生态学报, 2008, 28(8): 3937-3947.
Wang S Q, Yu G R. Ecological stoichiometry characteristics of ecosystem carbon, nitrogen and phosphorus elements[J]. Acta Ecologica Sinica, 2008, 28(8): 3937-3947. DOI:10.3321/j.issn:1000-0933.2008.08.054
[39] Gundersen P, Callesen I, De Vries W. Nitrate leaching in forest ecosystems is related to forest floor C/N ratios[J]. Environmental Pollution, 1998, 102(1): 403-407. DOI:10.1016/S0269-7491(98)80060-2
[40] Tian H Q, Chen G S, Zhang C, et al. Pattern and variation of C∶N∶P ratios in China's soils: a synthesis of observational data[J]. Biogeochemistry, 2010, 98(1-3): 139-151. DOI:10.1007/s10533-009-9382-0
[41] Koerselman W, Meuleman A F M. The vegetation N∶P ratio: a new tool to detect the nature of nutrient limitation[J]. Journal of Applied Ecology, 1996, 33(6): 1441-1450. DOI:10.2307/2404783
[42] 卞福花, 吴秋堂, 吴梦迪, 等. 不同水盐生境下芦苇湿地植被及土壤碳氮磷生态化学计量特征[J]. 应用生态学报, 2022, 33(2): 385-396.
Bian F H, Wu Q T, Wu M D, et al. C: N: P stoichiometry in plants and soils of Phragmites australis wetland under different water-salt habitats[J]. Chinese Journal of Applied Ecology, 2022, 33(2): 385-396.
[43] Meng L, Qu F Z, Bi X L, et al. Elemental stoichiometry (C, N, P) of soil in the Yellow River Delta nature reserve: Understanding N and P status of soil in the coastal estuary[J]. Science of the Total Environment, 2021, 751. DOI:10.1016/j.scitotenv.2020.141737
[44] Stevens R J, Laughlin R J, Malone J P. Soil pH affects the processes reducing nitrate to nitrous oxide and di-nitrogen[J]. Soil Biology and Biochemistry, 1998, 30(8-9): 1119-1126. DOI:10.1016/S0038-0717(97)00227-7
[45] 孙志高, 刘景双, 于君宝, 等. 湿地土壤NH3挥发、N2O释放过程及影响因素[J]. 湿地科学, 2008, 6(3): 429-439.
Sun Z G, Liu J S, Yu J B, et al. Ammonia volatilization, nitrous oxide emission processes in wetland soil and affecting factors[J]. Wetland Science, 2008, 6(3): 429-439.
[46] Lu J N, Feng S, Wang S K, et al. Patterns and driving mechanism of soil organic carbon, nitrogen, and phosphorus stoichiometry across northern China's desert-grassland transition zone[J]. CATENA, 2023, 220. DOI:10.1016/j.catena.2022.106695
[47] Wang X G, Lü X T, Zhang H Y, et al. Changes in soil C: N: P stoichiometry along an aridity gradient in drylands of northern China[J]. Geoderma, 2020, 361. DOI:10.1016/j.geoderma.2019.114087
[48] 王维奇, 王纯, 刘白贵. 盐度对湿地枯落物分解过程中碳氮磷化学计量比的影响[J]. 中国环境科学, 2012, 32(9): 1683-1687.
Wang W Q, Wang C, Liu B G. Effect of salinity on carbon, nitrogen and phosphorus stoichiometry during the decomposition of wetland litter[J]. China Environmental Science, 2012, 32(9): 1683-1687.
[49] 刘文龙, 谢文霞, 赵全升, 等. 胶州湾芦苇潮滩土壤碳、氮和磷分布及生态化学计量学特征[J]. 湿地科学, 2014, 12(3): 362-368.
Liu W L, Xie W X, Zhao Q S, et al. Spatial distribution and ecological stoichiometry characteristics of carbon, nitrogen and phosphorus in soil in Phragmites australis tidal flat of Jiaozhou Bay[J]. Wetland Science, 2014, 12(3): 362-368.
[50] Yan X J, Yang W H, Muneer M A, et al. Land-use change affects stoichiometric patterns of soil organic carbon, nitrogen, and phosphorus in the red soil of Southeast China[J]. Journal of Soils and Sediments, 2021, 21(7): 2639-2649.