硅是土壤和岩石的基本成分，是地壳上含量居第二位的元素，在全球生物地球化学循环和缓解全球气候变化方面发挥着重要作用[1]。硅可以提高植物光合作用的效率、改变营养元素的化学计量比、提高植物对害虫和病原体的抵抗力、增强植物对干旱和重金属的耐受能力，进而提高农作物的质量和产量[2-3]。硅循环和碳循环之间的交互作用可以调节大气中CO2的浓度，而通过河口输送至滨海水体中的可溶性硅的量影响滨海区域富营养化问题，因此硅的迁移转化过程对于维持地球化学循环过程具有重要意义[1,4]。海洋沉积物中的生物硅可以指示海洋古生产力的变化，硅藻控制的生物泵过程对大气CO2有净吸收效应进影响大气中CO2浓度[5]。植物生长过程中吸收硅元素可以形成植硅体，植硅体形成过程中封闭的那部分有机碳可以在环境中持留上千年，是重要的稳定性碳库[6]。湿地富含生物固定的无定形氧化硅，对全球硅循环具有重要影响[7]。一项对湿地植物封存植硅体碳能力的评估表明，全球约5.7×108hm2的湿地植物每年能固定大气中CO2的量高达4.39×107t[8]。滨海盐沼湿地可以有效捕集生物硅，具有较高的溶解性硅循环能力，因此相比于其他滨海水体具有较高的硅富集能力，而滨海潮间带被认为是滨海河口区域重要的硅汇和硅处理区[7]。尽管滨海盐沼湿地在硅的转化、迁移和储藏方面发挥着重要作用，但盐沼湿地硅的循环和储存还有待进一步了解。此外，了解硅元素在不同滨海湿地土壤中的分布特征及其影响因素对于深入了解滨海湿地碳汇功能具有重要意义。

土壤硅循环包含矿物风化过程、生物吸收、矿物相硅的吸附/解吸，因此土壤硅可以划分为不同的硅库[9]。土壤中的硅分为有机硅和无机硅，无机硅可以分为晶态硅(硅酸盐矿物、结晶二氧化硅)和可提取态硅(水溶态硅、交换态硅、胶体态硅和无定形硅),其中可提取态硅对植物的营养意义高于晶体硅[10]。由于不同植物对硅的需求不同，而且植物从土壤生物硅库吸收的硅量远超过从岩石风化释放吸收的硅量，无硅源补充的情况下，土壤中可提取态硅会因植被吸收和淋溶作用而逐渐损失[1,10]。在流域水平上，硅的物质量平衡受到易溶性硅的强烈影响[4],草酸提取态硅的量对土壤可利用硅具有正向影响[11]。

草酸提取态硅可以表征弱结晶硅、弱序性倍半氧化物中的硅和封闭于弱结晶铁氧化物中的硅，与弱结晶铁氧化物结合的硅在铁氧化物溶解时可以变成植物可利用性硅[4,12]。由溶解性硅酸和结合在矿物表面的硅组成的植物可利用性硅被认为是活性最高的硅组分[13]。连二硫酸钠-柠檬酸钠可以提取与结晶态次生铁氧化物结合的硅，焦磷酸钠可以提取与无定形铁铝氧化物结合的硅[9]。本研究选择非淹水芦苇湿地、潮汐芦苇湿地、淡水芦苇湿地和盐地碱蓬湿地作为研究样地，探究草酸提取态硅、连二硫酸钠-柠檬酸钠提取态硅和焦磷酸钠提取态硅在4种湿地中的含量变化及分布特征，研究结果可为丰富滨海湿地硅元素的分布信息提供有力支撑。

1、材料与方法

1.1 样品采集及测试分析

本研究的研究区位于山东省东营市的黄河三角洲自然保护区。黄河三角洲位于暖温带半湿润大陆性季风气候区，气候特征为四季分明、温差明显、雨热同期。年平均降雨量约为594 mm, 主要集中在6～9月，年平均蒸发量为2 049 mm, 年平均气温为11.7 ℃～12.6 ℃[14]。

在研究区内选择淡水芦苇湿地、潮汐芦苇湿地、非淹水芦苇湿地和盐地碱蓬湿地作为研究样地，于2018年10月采集0～20 cm深度土壤并按0～<5、5～<10和10～20 cm间隔分割土壤，每种湿地设置3个样点重复，每种湿地共计采集9个土样。将土样装入自封袋内带回实验室，剔除肉眼可见的石块、植物残体，室温条件下自然风干土样，经研钵研磨后，一部分土样过20目筛，用于测定土壤电导率和pH;另一部分过100目筛，用于测定草酸提取硅、连二硫酸钠-柠檬酸钠提取硅、焦磷酸钠提取硅、土壤有机质、总氮质量分数以及其他土壤指标。不同提取态硅元素的处理和测定方法如下：(1)草酸提取硅的测定。称取0.3 g过100目筛的土壤于三角瓶中，加入30 mL 0.2 mol/L的草酸-草酸铵混合溶液，调节pH=3以后，在暗处振荡4 h。(2)连二硫酸钠-柠檬酸钠提取硅质量分数。称取2 g研磨过100目筛的土壤于三角瓶中，加入1 g连二亚硫酸钠和50 mL 0.3 mol/L的柠檬酸钠溶液后，振荡16 h。(3)焦磷酸钠提取硅质量分数。称取0.3 g研磨过100目筛的土壤于三角瓶中，加入30 mL 0.1 mol/L的焦磷酸钠，并调节pH至10.5。以上提取过程结束后，加入0.2 mL 0.5 mol/L的MgCl2溶液，在1 000 r/min转速下离心15 min, 然后将上清液过0.45 μm滤膜，用ICP-OES测定滤液中硅的质量分数。土壤有机质质量分数、总氮质量分数、容重、电导率、pH、粉砂粒、砂粒和含水率的测定参照赵海晓等[15]的研究。

1.2 数据分析

采用单因素方差分析(one-way-ANOVA)对不同样地土壤草酸提取态硅、连二硫酸钠-柠檬酸钠提取态硅、焦磷酸钠提取态硅和土壤理化性质进行显著性差异分析，当p<0.05 时，认为具有显著性差异。在单因素方差分析中，当方差齐性时，采用Turkey检验；在方差非齐性时，采用Welch检验(校正的F检验)。皮尔森相关分析判定影响三种提取态硅的潜在控制因素。ANOVA和皮尔森相关分析采用SPSS 22.0 软件包来实现，土壤DOC含量变化图采用OriginPro 2021软件绘制。

2、研究结果

2.1 四种湿地土壤中3种易溶态硅和土壤理化性质的对比

四种湿地0～20 cm土壤中硅元素和土壤理化性质(取0～<5、5～<10和10～20 cm土壤的均值)之间的对比如表1所示。非淹水芦苇湿地土壤的草酸提取态硅质量分数的平均值为0.16 g/kg, 显著低于潮汐芦苇湿地(0.28 g/kg,p<0.05),但和淡水芦苇湿地、盐地碱蓬湿地土壤草酸提取态硅质量分数之间无显著性差异(p>0.05)。连二硫酸钠-柠檬酸钠提取态硅和焦磷酸钠提取态硅在4种湿地之间均无显著性差异(p>0.05),仅从数值上来看，两种硅质量分数的最小值均出现在非淹水芦苇湿地，分别为0.55和0.27 g/kg。潮汐芦苇湿地土壤连二硫酸钠-柠檬酸钠提取态硅质量分数的平均值在4种湿地中最高(0.76 g/kg),而焦磷酸钠提取态硅质量分数平均值的最高值出现在盐地碱蓬湿地(0.34 g/kg)。

表1 四种湿地0～20 cm土壤中硅和土壤理化性质的差异(平均值±标准差，n=36)

2.2 草酸提取态硅在各湿地不同土壤层次的分布特征

草酸提取态硅在4种湿地0～20 cm土壤中的剖面分布如图1所示。沿0～20 cm土壤剖面，草酸提取态硅在非淹水芦苇湿地呈现先上升后下降的趋势(p<0.05),最高值出现在5～<10 cm (0.20 g/kg),最低值均出现在10～20 cm土壤层(0.11 g/kg)。潮汐芦苇湿地和盐地碱蓬湿地0～<5和5～<10 cm土壤草酸提取态硅质量分数之间无显著性差异(p>0.05),但均显著高于10～20 cm土壤草酸提取态硅含量(p<0.05),这两种湿地土壤草酸提取态硅沿0～20 cm土壤剖面总体呈现下降趋势。淡水芦苇湿地0～<5、5～<10和10～20 cm土壤的草酸提取态硅之间无显著性差异(p>0.05),仅就数值来说，该湿地草酸提取态硅含量沿剖面呈逐渐下降趋势。

图1 不同湿地土壤中草酸提取态硅的分布

从不同土壤层来看，非淹水芦苇湿地0～<5 cm深度土壤的草酸提取态硅质量分数(0.16 g/kg)显著低于潮汐芦苇湿地(0.30 g/kg)和盐地碱蓬湿地(0.23 g/kg,p<0.05),而淡水芦苇湿地0～5 cm土壤草酸提取态硅质量分数(0.38 g/kg)与其他三种湿地之间均无显著性差异(p>0.05)。就5～<10 cm土壤层来看，非淹水芦苇湿地土壤草酸提取态硅质量分数显著低于潮汐芦苇湿地(0.33 g/kg,p<0.05),但这两类湿地与淡水芦苇湿地(0.20 g/kg)和盐地碱蓬湿地(0.35 g/kg)土壤草酸提取态硅质量分数之间均无显著性差异(p>0.05)。10～20 cm土壤层草酸提取态硅质量分数最高值出现在潮汐芦苇湿地(0.21 g/kg),显著高于非淹水芦苇湿地、淡水芦苇湿地(0.11 g/kg)和盐地碱蓬湿地(0.15 g/kg)土壤草酸提取态硅质量分数(p>0.05)。

2.3 连二硫酸钠-柠檬酸钠提取态硅在各湿地不同土壤层次的分布特征

连二硫酸钠-柠檬酸钠提取态硅质量分数在4种湿地0～20 cm土壤中的剖面分布如图2所示。非淹水芦苇湿地、潮汐芦苇湿地和淡水芦苇湿地10～20 cm土壤层连二硫酸钠-柠檬酸钠提取态硅均显著低于5～<10 cm土壤层(p<0.05),但0～<5 cm与5～<10 cm、10～<20 cm土壤连二硫酸钠-柠檬酸钠提取态硅的质量分数之间均无显著性差异(p>0.05),这三种湿地土壤连二硫酸钠-柠檬酸钠提取态硅质量分数沿0～20 cm剖面总体呈现下降趋势。淡水芦苇湿地0～<5 cm土壤连二硫酸钠-柠檬酸钠提取态硅质量分数与5～<10 cm土壤连二硫酸钠-柠檬酸钠提取态硅质量分数之间均无显著性差异(p>0.05),但均显著高于10～20 cm土壤连二硫酸钠-柠檬酸钠提取态硅质量分数(p<0.05)。

图2 不同湿地土壤中连二硫酸钠-柠檬酸钠提取态硅的分布

显著性差异，不同大写英文字母表征同一湿地中连二硫酸钠-柠檬酸钠提取态硅在不同土层之间的显著性差异，显著性水平为p<0.05。

非淹水芦苇湿地0～<5 cm土壤连二硫酸钠-柠檬酸钠提取态硅质量分数(0.67 g/kg)显著低于潮汐芦苇湿地(0.87 g/kg,p<0.05),但这两种湿地0～5 cm土壤连二硫酸钠-柠檬酸钠提取态硅与淡水芦苇湿地(0.78 g/kg)和盐地碱蓬湿地土壤连二硫酸钠-柠檬酸钠提取态硅(0.83 g/kg)质量分数之间均无显著性差异(p>0.05)。非淹水芦苇湿地和淡水芦苇湿地5～<10 cm土壤连二硫酸钠-柠檬酸钠提取态硅的质量分数显著低于潮汐芦苇湿地和盐地碱蓬湿地(p<0.05),但非淹水芦苇湿地和淡水芦苇湿地土壤连二硫酸钠-柠檬酸钠提取态硅质量分数之间、潮汐芦苇湿地和盐地碱蓬湿地土壤连二硫酸钠-柠檬酸钠提取态硅之间均无显著性差异(p>0.05)。10～20 cm土壤连二硫酸钠-柠檬酸钠提取态硅的最高值出现在潮汐芦苇湿地(0.94 g/kg),显著高于其他三种湿地(p<0.05);最低值出现在非淹水芦苇湿地和淡水芦苇湿地，而且这两种湿地土壤连二硫酸钠-柠檬酸钠提取态硅质量分数显著低于盐地碱蓬湿地(p<0.05)。

2.4 焦磷酸钠提取态硅在各湿地不同土壤层次的分布特征

焦磷酸钠提取态硅在4种湿地0～20 cm土壤中的剖面分布如图3所示。非淹水芦苇湿地、潮汐芦苇湿地、盐地碱蓬湿地0～<5 cm和5～<10 cm土壤焦磷酸钠提取态硅质量分数之间无显著性差异(p>0.05),但均显著高于10～20 cm土壤焦磷酸钠提取态硅质量分数(p>0.05)。这3种湿地土壤焦磷酸钠提取态硅质量分数沿0～20 cm土壤剖面呈下降趋势。淡水芦苇湿地0～<5 cm土壤焦磷酸钠提取态硅质量分数与5～<10 cm、10～20 cm土壤焦磷酸钠提取态硅质量分数之间均无显著性差异(p<0.05),但5～<10 cm土壤焦磷酸钠提取态硅质量分数显著高于10～20 cm土壤(p>0.05)。

图3 不同湿地土壤中焦磷酸钠提取态硅的分布

如图3所示，4种湿地0～<5 cm土壤焦磷酸钠提取态硅质量分数之间无显著性差异(p>0.05)。非淹水芦苇湿地和淡水芦苇湿地5～<10 cm土壤焦磷酸钠提取态硅质量分数显著低于和盐地碱蓬湿地(p<0.05),这3种湿地0～<5 cm土壤焦磷酸钠提取态硅质量分数与潮汐芦苇湿地之间均无显著性差异(p>0.05)。在10～20 cm土壤层，非淹水芦苇湿地和淡水芦苇湿地焦磷酸钠提取态硅质量分数显著低于潮汐芦苇湿地和盐地碱蓬湿地(p<0.05),而非淹水芦苇湿地和淡水芦苇湿地焦磷酸钠提取态硅质量分数之间、潮汐芦苇湿地和盐地碱蓬湿地焦磷酸钠提取态硅质量分数之间均无显著性差异(p>0.05)。

2.5 环境因子对不同形态硅的影响分析

草酸提取态硅、连二硫酸钠-柠檬酸钠提取态硅和焦磷酸钠提取态硅质量分数与土壤理化性质之间的相关关系如表2所示。草酸提取态硅质量分数与土壤有机质、总氮、粉砂粒质量分数和含水率呈显著正相关关系(p<0.01),与容重、pH和砂粒呈显著负相关关系(p<0.01)。连二硫酸钠-柠檬酸钠提取态硅质量分数与土壤有机质、总氮和粉砂粒质量分数之间具有显著正相关关系(p<0.05),与容重、pH和砂粒质量分数之间具有显著负相关关系(p<0.01)。焦磷酸钠提取态硅质量分数与土壤有机质、总氮、粉砂粒质量分数和含水率之间具有显著正相关关系(p<0.05),与pH和砂粒质量分数之间具有显著负相关关系(p<0.01)。

表2 不同形态硅质量分数与土壤理化性质之间的相关性分析

3、讨论

3.1 植被类型和水文条件的影响

湿地沉积物中的硅主要来自于淹水输入的硅、底栖原生沼泽硅藻产生的硅和植物腐败后埋入土壤中的植硅体。湿地系统内硅的生物地球化学循环过程是由可以富集硅元素的植物周转促成的[16]。因此，植物类型是影响湿地硅循环的重要因素[4]。芦苇因具有较高的生产力、密集无性繁殖系生长能力和干枯植物体留存时间长的特点，被认为是潜在的营养库，而且芦苇湿地的土壤生物硅质量分数高于凤仙花湿地、柳树湿地和柳叶菜湿地[7]。此外，湿地水文条件(淹水频率和淹水持留时间等)在硅的转化过程也发挥着控制作用。淹水可以输入大量含有植硅体和草酸提取态硅的悬浮颗粒物，频繁的淹水还可以促进植硅体的溶解，而植硅体的溶解是湿地土壤中重要的硅源[16]。Li等[8]认为水位也是影响土壤硅组分的重要因素，水位管理可以通过促进有机质和植硅体的累积、提高土壤-植物之间的硅循环来提高土壤硅的生物有效性[8]。淹水造成的还原环境使得与硅结合的铁被还原从而释放出更多的硅。本研究中，就平均值来说，盐地碱蓬湿地0～20 cm土壤3种提取态硅质量分数与水文条件不同的芦苇湿地之间并无显著性差异，而且在部分土层盐地碱蓬湿地土壤3种提取态硅质量分数高于非淹水芦苇湿地，这表明除植被类型以外的因素如潮汐淹水和其他因素也会影响土壤硅质量分数。

3.2 土壤理化性质的影响

土壤理化性质是影响土壤硅组分的重要因素[7,18]。土壤有机质不仅可以为单硅酸提供吸附位点，其分解过程可以产生有机酸和还原性环境从而促进结合于铁锰氧化物中硅的释放[18]。因此，本研究中草酸提取硅、连二硫酸钠-柠檬酸钠提取硅和焦磷酸钠提取硅均和土壤有机碳呈显著正相关关系，表明土壤有机碳质量分数对土壤3种易溶态硅具有正向的促进作用。邱思婷等[17]对闽江口湿地的研究发现，有效硅质量分数与土壤容重呈显著负相关关系。本研究中草酸提取态硅和连二硫酸钠-柠檬酸钠提取态硅均与土壤容重呈显著负相关关系。这可能是因为较高的土壤容重会导致土壤硬度增加从而增加根系生长阻力、降低地下生物量进而导致硅的累积量减少[17]。土壤pH是影响土壤硅的重要因素，一般认为在碱性土壤中，pH与土壤硅质量分数呈显著负相关关系[19]。本研究中3种提取态硅均与pH呈显著负相关关系，可能是因为4种湿地pH介于8.14和9.73之间。草酸提取态硅和焦磷酸钠提取态硅均与土壤含水率呈显著正相关关系，这可能是因为较高的土壤含水率有利于水溶性硅在土壤中的持留，而可提取态硅组分之间可以相互转化[20]。因此，相比于潮汐芦苇湿地，非淹水芦苇湿地土壤中较低的土壤有机质、总氮、含水率和粉砂粒以及较高的容重和pH也是造成该湿地草酸提取态硅低于潮汐芦苇湿地的原因。

3.3 剖面分布趋势的影响因素

与本研究结果相似，Struyf等[7]发现芦苇湿地土壤在不同深度的剖面上，其土壤生物硅质量分数呈梯度下降的趋势，0～<10 cm土壤生物硅质量分数高于10～20 cm。这可能与植硅体主要聚集在湿地土壤表层并随土壤深度增加而逐渐下降有关[13]。此外，风化作用、枯落物返还、地表径流和潮汐作用的影响对湿地下层土壤的影响弱于表层土壤，这也是决定3种提取态硅剖面分布趋势的重要原因[17]。

4、结论

研究结果表明，盐地碱蓬湿地和3种芦苇湿地土壤草酸提取态硅、连二硫酸钠-柠檬酸钠提取态硅和焦磷酸钠提取态硅之间无差异，说明植被类型对3种提取态硅的影响较小。非淹水芦苇湿地土壤3种提取态硅显著低于潮汐芦苇湿地，但与淡水芦苇湿地和盐地碱蓬湿地之间无差异，说明就水文条件差异是影响3种提取态硅含量和分布的因素之一。除植被类型和水文条件下，有机质、容重、pH、钠离子、氯离子、粉砂粒、砂粒和含水率是影响3种提取态硅的重要环境因子。

参考文献:

[1]刘丽君,黄张婷,孟赐福,等.中国不同生态系统土壤硅的研究进展[J].土壤学报,2021,58(1):31-41.

[5]杨功旭,姚政权,冯晗,等.北极东西伯利亚陆架表层沉积生物硅分布特征及环境意义[J].矿物岩石地球化学通报,2023,42(4):717-727.

[10]赵送来,宋照亮,姜培坤,等.西天目集约经营雷竹林土壤硅存在形态与植物有效性研究[J].土壤学报,2012,49(2):331-338.

[14]杜书栋,白军红,贾佳,等.黄河三角洲芦苇湿地土壤有机碳储量沿盐分梯度的变化特征[J].环境科学学报,2022,42(1):80-87.

[15]赵海晓,高永超,赵庆庆,等.不同水文条件下黄河三角洲湿地土壤溶解性有机碳的分布特征[J].北京师范大学学报(自然科学版),2021,57(1):51-58.

[17]邱思婷,米慧珊,高会,等.闽江河口湿地不同植被带土壤全硅的含量及分布特征[J].生态学报,2020,40(22):8306-8314.

[19]李家书,谢振翅,胡定金,等.湖北省土壤有效硅含量分布[J].热带亚热带土壤科学,1997,6(3):176-181.

[20]何冬梅,陈逸飞,苏仪,等.郭岩山不同海拔天然栲树林土壤硅形态特征[J].林业科学研究,2023,36(2):153-160.

基金资助:山东省重点研发计划(重大科技创新工程)(2021CXGC011201); 国家重点研发计划(2023YFC3709004-03); 山东省自然科学基金(ZR2022MD116); 济南市“高校20条”资助项目(2020GXRC059);

文章来源:季晓晖,赵庆庆,展海银,等.三种提取态硅在典型滨海湿地中的分布特征及影响因素分析[J].山东科学,2024,37(05):95-102.