垦殖对黄河三角洲盐渍土碳氮分布特征的影响*

李贤红 陈为峰 宋希亮 王曼华 胡 琴 邓 从

(山东农业大学资源与环境学院，山东泰安 271018)

随着土壤资源的日益紧张以及人口数量的剧增，盐碱土资源利用与改良成为世界性的重大课题之一［1］。因此，研究盐碱地开发利用后的土壤性质变化对指导盐碱地开发具有重要的意义。国内外关于盐渍土开发利用的研究多关注水盐运移及盐分改良［1-3］。近年来，有关盐渍土开发利用的研究集中于土地利用方式以及管理模式［4-5］，研究多关注表层土壤［6-8］，且对垦殖历史的研究较短［8-9］，对于深层次、大时间尺度的农业利用碳氮分布规律研究尚欠缺。黄河三角洲作为中国最年轻的土地，区域土壤形成发展的历史，实际就是土地不断被垦殖利用的历史。近年来，盐碱地开发利用进一步加快，基于强灌强排的主导治理措施和高潜水位的特征可能会对区域垦殖利用过程中土壤碳变化产生较大的影响。前期，主要开展了黄河三角洲滨海新生湿地碳氮循环的研究［10-12］，而对盐碱荒地垦殖前后土壤碳行为过程的研究尚未引起足够的重视。土壤无机碳作为碳循环过程的一部分，紧密联系全球碳循环，如果忽视无机碳，不利于精确和正确估算碳储存的变化［13］。土地利用方式显著影响土壤无机碳变化过程［14］和土壤碳氮等生源要素的耦合关系［15-16］。随着黄河三角洲高效经济发展以及“渤海粮仓”建设，研究大尺度的盐渍土变化特征，对后续土地利用具有重要的现实指导意义。

本文选取黄河三角洲典型不同垦殖年限土壤为研究对象，以样地空间变化代替时间动态，系统研究盐碱地垦殖前后及不同垦殖历史典型土层的土壤碳（有机碳、无机碳、全碳）、氮（全氮、碱解氮）的分布特征，阐明变化规律，有利于从整体上准确评估土地利用剧烈变化情景下的土壤碳氮响应过程，为区域土壤碳氮循环模型的建立提供思路和参考数据，还可为区域未利用盐碱地的合理开发与科学管理提供指导，具有重要的理论意义和实践价值。

1 材料与方法

1.1 研究区概况

研究区位于黄河三角洲东营市仙河镇，地理坐标118°35′57″E～118°55′49″E，37°53′20″N～37°56′22″N。选取开垦不同年限的农用地及未垦殖盐碱荒地作为研究样地。该区域年平均气温13.3℃，平均无霜期206 d，年均降水量537 mm，年平均蒸发量1 885 mm。区域土壤类型主要为盐化潮土，土壤质地类型以粉砂壤土为主，含有砂质壤土、黏壤土。研究区内农作物以棉花、玉米、小麦为主，自然植被以芦苇、碱蓬、蒿类为主。地下水矿化度在空间上分布的总体规律表现为，由内地向沿海方向地下水矿化度逐渐增加，近内陆部分矿化度较低，一般2 g L-1左右；沿海地带地下水矿化度较高，局部地区甚至大于50 g L-1，其他大部分地区地下水矿化度大于10 g L-1。研究区灌溉水源来自黄河，灌水方式为大水漫灌，每年平均洗盐灌水量约为2 250 m3hm-2。

1.2 野外调查与土壤样品采集

研究区位于1934年海岸线以下，根据成国栋［17］对黄河三角洲形成时间研究，黄河三角洲按照形成顺序分为11个叶瓣，研究区位于第9个叶瓣上，并与第10、11叶瓣有所重叠，陆地形成时间约在1934年以后。其成土年龄与成土环境接近。1963年成立军马场开始开发［18］，至今盐碱地垦殖历史已超过50 a，可视为160多年以来整个黄河三角洲盐碱地开发进程的一个典型阶段，是开展黄河口地区土地利用变化及生态效应研究的最佳天然实验室。采样区根据土地利用变化研究有关文献［19］，以及土地现状图结合实地调查，以未垦殖盐碱荒地为对照，分别选取2014年开发的耕地（垦殖小于5 a）、2005年开发的耕地（垦殖约10～15 a）、1995年开发的耕地（垦殖约20～25 a）、1985年开发的耕地（垦殖约30～35 a）和1956年左右开发的耕地（垦殖约50～60 a），确定了不同垦殖年限区域土壤进行研究，详见图1。不同垦殖年限区域基本概况见表1。

于2016年8月，选择未垦殖盐碱荒地裸地及垦殖不同年限（小于5 a、10～15 a、20～25 a、30～35 a、50～60 a）的6个典型地块样地，土地利用方式均为小麦/玉米连作，以其1 m深度土壤为研究对象，利用直径3 cm土钻采集0～100 cm土壤样品。每个垦殖年限选择三个相距不小于1 000 m的地块，每个地块采用五点取样法，每个样点为相距不小于5 m的三钻混合，每个点面按0～10、10～20、20～30、30～40、40～60、60～80、80～100 cm的7个层次进行采样，同时分层取环刀土壤样品。

1.3 土壤理化性质分析

土样风干，磨细过0.15 mm筛，使用元素分析仪（Elementar Vario Macro，德国）测定土壤全氮、全碳［19］，使用1 mol L-1的盐酸去除土样中的碳酸盐后，使用总有机碳分析仪（Elementar Vario TOC，德国）测定土壤有机碳［20］，无机碳为全碳与有机碳差减值，并通过气量法验证。土样过1 mm筛后，参照《土壤农业化学分析方法》［21］，1︰5土水比浸提液，质量法测定土壤盐分，电位法测定pH，烘干测定土壤容重，碱解扩散法测定碱解氮，碳酸氢钠浸提—钼锑抗比色法测定有效磷，火焰光度法测定速效钾。

图1 取样点分布图Fig. 1 Distribution of sampling points

表1 不同垦殖年限区域概况Table 1 General information of the sample plots different in cultivation history

1.4 数据处理

采用Microsoft Excel 2007处理数据，用SPSS 19.0软件进行单因素方差分析和显著性分析（最小显著差异（LSD）法检验不同处理的差异性），相关图表制作使用Origin 9.0完成。

2 结 果

2.1 垦殖年限对盐碱地土壤碳的影响

土壤基本理化性状见表2。通过分析不同垦殖年限土壤有机碳、无机碳和全碳含量的变化（图2），结果发现，未垦殖荒地对照各土层之间差异不显著（p＞0.05），盐碱地垦殖显著增加1 m土体土壤有机碳、无机碳、全碳含量。通过进一步对1 m土壤碳含量进行方差分析，结果发现，不同垦殖年限土壤有机碳、无机碳、全碳含量均显著高于未垦殖盐碱荒地（p＜0.05），除垦殖10～15 a的无机碳，均随着垦殖年限的增加整体呈增加趋势，土壤碳（有机碳、无机碳、全碳）含量在垦殖30～35 a后达到相对稳定的水平。

通过分析不同垦殖年限土壤有机碳、无机碳和全碳的垂直分布特征，发现盐碱土垦殖主要影响表层（0～40 cm）碳分布，随着土层深度的增加该影响降低，具有一定的表聚性，40～100 cm碳含量无明显变化规律。不同垦殖年限40～100 cm平均土壤无机碳在一定程度上受地下水埋深（表1）的影响。随着垦殖年限的增加，直至垦殖30～35 a以上，表层（0～40 cm）土壤有机碳显著高于底层（40～100 cm）土壤有机碳（p＜0.05）。个别处理（垦殖＜5 a、20～25 a）土壤有机碳40～100 cm含量显著高于其他处理（p＜0.05）；土壤无机碳在垦殖10～15 a处理的0～40 cm土层出现降低，在＜5 a、20～25 a和30～35 a的深层土壤（40～100 cm）有一定的聚集。

表2 不同垦殖年限土壤基本理化性状（0～1m）Table 2 Basic physical and chemical propertiesof the soil in sample plot relative to cultivation history（0～1m）

图2 不同垦殖年限不同土层土壤有机碳、无机碳和全碳含量Fig. 2 Contents of soil organic carbon, inorganic carbon and total carbonin sample plot, relative to cultivation history

2.2 垦殖年限对盐碱地土壤氮的影响

不同垦殖年限全氮变化范围在0.17～1.27 g kg-1之间，平均值为0.51 g kg-1（图3A）。通过对1 m土层全氮方差分析，盐碱荒地垦殖后，全氮含量显著增加，且随垦殖年限的延长而增加，具体表现为，垦殖20～25、30～35、50～60 a显著高于＜5、10～15 a。

通过分析土壤全氮垂直变化规律，结果表明，随着土层深度的增加，不同垦殖年限处理全氮均表现降低趋势。未垦殖盐碱荒地0～10 cm全氮显著高于10～100 cm土层（p＜0.05）。表层（0～40 cm）土壤全氮含量与未垦殖盐碱地相比均显著提高，随着垦殖年限的增加表聚性特征更加明显（除20～25 a），个别年限（20～25 a）全氮含量出现底层（40～100 cm）累积现象。垦殖30 a后表层（0～40 cm）全氮含量显著高于底层（40～100 cm）。

图3 不同垦殖年限不同土层土壤全氮(A)、碱解氮(B)含量Fig. 3 Contents of total nitrogen (A) and alkaline nitrogen (B)in sample plot relative to cultivation history

由图3B可见，不同垦殖年限土壤碱解氮平均值为30.75 mg kg-1。垦殖初期，土壤碱解氮大幅增加，随着垦殖年限的增加，碱解氮增幅降低。碱解氮变化规律与全氮变化规律相似。通过对1 m土壤碱解氮方差分析发现，不同垦殖年限处理土壤碱解氮含量显著高于未垦殖盐碱荒地（p＜0.05），具体表现为：随着垦殖年限的增加1 m土体土壤碱解氮含量整体呈增加的趋势；垦殖年限50～60 a有机碳含量显著高于30～35 a，30～35 a显著高于10～25 a，10～15 a显著高于＜5 a。上述结果说明，盐碱地垦殖后会显著增加1 m土壤土体碱解氮含量。

不同垦殖年限处理下，1 m土体土壤碱解氮垂直分布特征表现为，随着土层深度的增加土壤碱解氮含量整体呈降低趋势；垦殖年限对碱解氮的影响主要作用于表层（0～30 cm），荒地0～10 cm碱解氮显著高于10～100 cm；随着垦殖年限的增加，至垦殖10 a之后0～30 cm土壤碱解氮显著高于30～100 cm（p＜0.05）。底层（30～100 cm）土壤碱解氮含量在垦殖初期出现一定程度的增加，具体表现为，垦殖初期5～25 a增加，30～60 a后碱解氮维持在稳定水平。

2.3 垦殖盐碱地土壤中碳氮相关性

通过对0～40 cm土层有机碳、全氮变化分析（图4A），发现随着垦殖年限的增加，土壤有机碳、全氮含量均呈增加趋势（除20～25 a），两者之间变化趋势高度一致。本研究发现，土壤垦殖主要影响0～40 cm土壤有机碳、全氮含量。

由图4B可以看出，不同垦殖年限各土层土壤有机碳与全氮呈极显著的（n=126，p＜0.01）线性正相关关系，且相关系数较高，二者之间可以用一元线性方程（y=0.1109x-0.012 8）进行拟合。这表明土壤有机碳会促进全氮的增加，土壤全氮的增加也会促进有机碳的增加，改善土壤养分含量，应注重两者之间的相关性。

2.4 垦殖年限对盐碱地土壤碳氮比（C/N）的影响

分析土壤C/N的变化有助于深入理解土壤有机碳氮的积累过程及土壤质量的变化趋势。由图5可知，不同垦殖年限盐碱地土壤碳氮比为6.56～16.87，垦殖后碳氮比平均值（9.35）较未垦殖荒地（10.89）有所降低；通过对1 m土体C/N方差分析，盐碱荒地垦殖后，土壤C/N显著降低，C/N随着垦殖年限的增加出现先降低后稳定的趋势，具体表现为：盐碱荒地开垦后至垦殖20～25 a，土体C/N降低至8.41，垦殖30～35 a C/N维持在稳定水平。根据表层（0～40 cm）以及底层（40～100 cm）碳氮比方差分析，虽然垦殖后随着垦殖历史的变化C/N出现波动，最终30～35 a后维持在9.2左右。

3 讨 论

3.1 不同垦殖年限下盐碱地土壤碳变化的原因

不同垦殖年限1 m土体有机碳含量随土层深度的增加而降低，这与夏江宝等［6］对0～40 cm盐碱地土地利用的研究结果一致。研究发现，不同垦殖年限土壤有机碳含量在40～100 cm无显著性变化趋势，该研究结果与李晓光等［9］对枸杞、柽柳能够增加1 m土体各土层有机碳的研究结果不同，其原因可能是土地利用方式不同，本研究作物为小麦/玉米连作，根系主要分布在表层（0～40 cm）。40～100 cm可能是该区域成土母质本身有机碳的含量普遍偏低。此外，垦殖年限＜5 a、20～25 a 的40～100 cm层次有机碳含量显著高于其他处理，出现深层积累，其原因可能是垦殖＜5 a、20～25 a的地下水埋深较其他处理浅（表1），深层土壤受水文频繁变化的累积影响。

图4 0～40 cm土层有机碳、全氮随着垦殖年限的变化（A）及其相关性（B）Fig.4 Total nitrogen and organic carbon contents in the topsoil layer (0～40 cm) of the sample plot relative to cultivation history and their relationship (B)

图5 不同垦殖年限土壤碳氮比Fig. 5 Soil C/N in sample plot relative to cultivation history

本研究发现，无机碳占全碳的比重更大，不同垦殖年限的土壤无机碳、全碳及有机碳具有相似的变化趋势。本研究中，表层（0～40 cm）无机碳含量显著高于底层（40～100 cm）。该研究结果与邓彩云等［22］的研究结果相反，分析其原因，一方面可能为区域以及成土母质的差异；另一方面，灌溉可能会改变土壤碳输入、输出以及储量，黄河水作为区域唯一灌溉水源，无机碳含量高达35 mg L-1以上［23］。垦殖10～15 a土壤无机碳较低，可能是因为该研究区地下水位较深（表1）。因此，研究无机碳储量不能脱离土地所属区域的地下水位以及随时间的变化情况等因素，尤其是以黄河三角洲为代表的高潜水位滨海盐渍土。

3.2 不同垦殖年限下盐碱地土壤氮变化的原因

在垂直方向上，随着土体深度的增加，各取样点的土壤全氮含量逐渐降低，这与夏志坚等［24］研究结果相似。随着垦殖年限增加，表层（0～40 cm）土壤全氮含量均表现为提高，可能是由以下原因共同造成的，一是土壤垦殖后施肥增加土壤氮含量［25］，增加了全氮的来源；二是由于植物残体和凋落物输入在土壤表层发生累积［26］；三是氮沉降因素［27］。此外，全氮含量在40～100 cm变化不显著，该结果可能与该区域成土母质中氮含量普遍偏低有关。研究中垦殖年限20～25 a的土壤，有机碳、全氮含量相对其他垦殖年限有所降低，一方面，可能由于该区域近两年出现一定程度的撂荒，种植管理粗放，秸秆还田量减少；另一方面，可能与有机肥的矿化分解有关。

碱解氮是作物易吸收的有效部分［21］，因而与作物在生产过程中从土壤中吸收的氮量相关性显著。土壤碱解氮的测定在农业生产上可直接用来评价土壤的供氮强度，制定氮肥施用方案。研究表明，土壤垦殖主要对0～30 cm碱解氮有显著影响，且随垦殖年限的增加而增加。相比于其他养分指标，更趋表聚化，其原因可能是垦殖施肥、秸秆还田等田间管理措施会增加碱解氮含量［28］。

3.3 不同垦殖年限下碳氮相关性及碳氮比变化的原因

土壤碳氮之间的耦合关系能更全面地阐述土壤碳氮变异，土壤C/N空间变异特征对区域环境保护与耕地可持续利用意义重大［27］。不同垦殖年限1 m土壤有机碳与全氮显著线性正相关，该结果与前人［29］研究盐碱地表层土（0～40 cm）碳氮的结果类似，研究区域各样地各层的C/N平均为9.39，略低于全国平均水平10.2［30］。土壤C/N随垦殖年限增加而逐渐降低，垦殖30～35 a后维持在较低水平，一方面，土壤有机碳含量随垦殖年限有所增加，但在长期采用强排强灌等灌水措施的作用下，可能增加有机碳淋失；另一方面可能与氮肥利用率较低有关。

土壤经过不同年限的垦殖，土壤碳（有机碳、无机碳、全碳）、氮（全氮、碱解氮）、C/N均在垦殖30～35 a后达到相对稳定的水平，一方面，可能是由于20世纪80年代开始大量使用化肥，前期以有机肥为主，养分含量低，另一方面，可能是由于施肥种植等垦殖方式可改善土壤粒径分布结构，进而增加土壤的保肥力，在垦殖30 a时均达到相对稳定的水平［31］。

本研究主要围绕黄河三角洲区域盐碱地的垦殖年限等状况对碳氮含量分布的影响展开。在研究表层土碳氮分布的同时，也对深层次的土壤碳氮的垂直分布进行了探索，土壤垦殖有利于增加土壤碳（有机碳、无机碳、全碳）、氮（全氮、碱解氮）含量，说明垦殖对黄河三角洲盐渍土土壤养分状况有一定的改良作用。本研究发现，研究区小麦/玉米连作垦殖30～35 a以上，土壤养分基本维持稳定水平，但最终表现为C/N降低，因此，针对区域盐渍土不同垦殖年限的C/N状况，应当通过改善施肥方式、均衡施肥等措施，提高养分利用率，使C/N维持在利于农业高效利用的水平，这对于未来黄河三角洲盐渍土垦殖利用、渤海粮仓的建设具有重要的指导意义。

4 结 论

荒地经垦殖后，土壤有机碳、无机碳、全碳、全氮、碱解氮含量显著提高，且随垦殖利用年限的增加而增加，并在垦殖30～35 a后达到相对稳定的水平；1 m深土体碳氮垂直分布规律具有一定的表聚现象，随土层深度的增加而降低；不同垦殖年限主要影响耕层（0～40 cm）土壤有机碳、无机碳、全碳、全氮，而在40～100 cm层次的含量变化不显著；垦殖年限主要影响0～30 cm土层碱解氮；地下水埋深对养分深层累积有一定的影响。盐渍土不同垦殖年限土壤有机碳与全氮之间显著正相关，但垦殖后C/N呈降低趋势，约30 a后基本稳定，C/N低于全国平均水平，说明长期垦殖土壤亟需进一步提升有机碳含量水平。