-
横坡垄作是在坡耕地上沿等高线耕作，沟垄相间，因能拦蓄雨水、增加就地入渗、减少地表径流和土壤流失[1]，成为东北黑土区、北方土石山区、西南土石山区、南方红壤丘陵区等区域广泛分布的耕作措施[2]。但横坡起垄时，垄沟方向会偏离等高线，致使垄沟内易形成低洼区域，从而诱发垄沟内雨水的汇聚和泥沙沉积。积水引发侧向二维入渗，致使土壤氮发生垂直和侧向运移。同时，垄沟内低洼处泥沙沉积导致土壤氮的“富集”，但在径流作用下会被再次侵蚀与输移，形成“富集”与“迁移” 交替局面。横坡垄作坡面土壤氮迁移的二维性、富集 / 迁移的交替性，可造成垄面、垄沟等特征点位土壤氮分布的差异，进而使得土壤氮空间分布特征明显有别于传统耕作坡面。此外，土壤氮作为表征耕地质量和土壤肥力的关键指标之一，易受到坡位、坡度等地形因素的影响[3-4]，而呈现出高度的空间变异特性。这使得横坡垄作坡面土壤氮的空间异质性更为复杂。明确横坡垄作坡面土壤氮空间分布特征是实现该措施下精准农业的基础[5]，同时可为土壤质量的改良和恢复提供重要理论依据。
-
目前，国内学者主要是利用地统计学、地理信息系统（GIS）等手段解析土壤氮在区域尺度上的空间分布异质性，如：黄土高原区、东北黑土区、紫色丘陵区、沙漠地带、喀斯特地区等[6-11]，并探析了田间管理（施肥、耕作方式）、土地利用方式、地形等因素对土壤氮空间分布的影响[12-17]，但忽略了小尺度微地形下的土壤氮的空间差异性。针对坡面尺度上土壤氮的研究主要集中于土壤氮随径流、泥沙的迁移规律与流失形态[18-21]，关于横坡垄作措施下坡面尺度上土壤氮的空间异质性尚不明确。鉴于中国坡耕地面临的土壤质量下降的严峻形势及田间精准农业的全面推广，急需针对坡耕地小尺度上土壤氮的空间变异进行深入研究，尤其是被广泛应用的横坡垄作措施下土壤氮的空间异质性。
-
北方土石山区地表土石混杂、石多土少，地面极易砂砾化或石化，渗透性强、土壤涵蓄水能力差，土壤养分贫瘠，山丘区耕地土层一半以上不到 50 cm[22]，是中国水土流失治理的重点区域[23]。此外，为获取农业高产，该区农民大量添加以 N、 P 为主的肥料，化肥年平均施用水平达 482 kg/hm2，远超发达国家设置的 225 kg/hm2 安全上限。因此，本研究以北方土石山区横坡垄作坡面为研究对象，样带法采样，采用经典统计学、地统计学相结合的方法，探究垄面、垄沟等特征点位不同坡位处全氮 （TN）、硝态氮（NO3--N）、铵态氮（NH4 +-N）水平和垂直方向上的分布规律，揭示横坡垄作坡面土壤氮的空间分布特征，以便为北方土石山区农业生产的精细化提供理论依据和科学指导。
-
1 材料与方法
-
1.1 研究区概况
-
北方土石山区主要包括海河流域和淮河流域，总面积约为 58 万 km2，山地丘陵面积达 48%。其中，淮河流域的坡耕地达 84 万 hm2，区内 15°以上的坡耕地平均每年流失 0.5～1.0 cm 表土[24]。沂蒙山区是北方土石山区的典型代表，位于山东省鲁中南低山丘陵区，丘陵山地占 84%[25]，隶属于淮河流域。研究区位于山东省蒙阴县双河峪小流域（118°7′30″E，35°38′26″N），地处沂蒙山区重要水源地-云蒙湖库区。该小流域总面积 2.50 hm2，属暖温带大陆季风气候，多年平均气温 12.8℃，平均年降水量 820 mm，降雨主要集中于 6—9 月。地貌以低山丘陵为主，土壤以花岗岩发育下的棕壤为主。耕地、园地、林地、荒坡地是该流域主要的土地利用方式，其中耕地占流域总面积的 27.5%。流域内 80% 以上的坡耕地采取横坡垄作耕作措施，该措施下的垄宽约为 80 cm、垄间距 20 cm。流域内广泛采用的横坡垄作措施为沂蒙山区乃至北方土石山区花生、地瓜的典型种植模式。
-
1.2 土壤样品采集及处理
-
以双河峪小流域内典型横坡垄作花生种植坡面为采样对象，该坡面南北长 35 m、东西宽 46 m。花生收获后，按照条带法在 10、22.5、36 m 处设置 3 条采样带，分别在每条采样带的坡上、坡中、坡下和坡脚 4 个坡位相邻 50 cm 的垄沟和垄面处采样。每个样点分 7 层深度采样：0～5、5～10、10～15、15～20、20～30、30～40 和 40～50 cm，并利用罗盘仪测定采样点的坡度。具体采样布设如图1 所示。
-
采集样品去除地表枯落物、石块等杂质且自然风干后，分别研磨过 0.1 和 0.25 mm 筛。过筛后，样品中的 TN 采用凯氏定氮法测定，NO3--N 和 NH4 +-N经 1 mol/L KCl 提取后，在紫外分光光度计上分别采用双波长系数法和氨氮-钠氏比色法进行测定。
-
图1 横坡垄作坡面采样布设
-
1.3 数据处理
-
基于 SPSS 22.0 对研究区 TN、NO3--N 和 NH4 +-N 含量进行经典统计学分析，并利用该软件中的方差分析（ANOVA）对坡位间土壤氮含量进行显著性检验，基于最小显著差数法（LSD）开展多重比较，并检验在 0.05 水平上的显著性。采用 Origin 8.6 中的单样本非参数检验（Shapiro-Wilk）对采样点土壤氮含量数据进行正态性检验，并利用此软件绘制柱状图和散点图。利用 Surfer12.0（美国 Golden Software）中的克里格方法（Kriging）进行插值，生成土壤氮空间分布等值线图，并绘制土壤氮 3D 曲面图。
-
2 结果与分析
-
2.1 横坡垄作坡面土壤氮的描述性统计特征
-
基于采样点 0～20 cm 土层深度内 TN、NO3--N 和 NH4 +-N 的含量进行经典统计学分析，获取横坡垄作坡面、垄面、垄沟等特征点位土壤氮的描述性统计特征。分析发现，垄面土壤氮含量的变化明显有别于垄沟。垄面 TN、NO3--N 和 NH4 +-N 的含量分别为 0.39～1.03 g/kg、2.88～8.63 mg/kg 和 7.03～13.31 mg/kg（表1），而垄沟相应含量分别为 0.48～0.86 g/kg、3.12～8.07 mg/kg 和 5.46～13.54 mg/kg。垄面 TN、NO3--N 和 NH4 +-N 的平均含量分别为 0.68g/kg、5.06 mg/kg 和 9.92 mg/kg，较垄沟分别增加 3.67%、5.19% 和 7.23%。参照全国土壤养分含量分级标准（TN 中水平 0.75～1.00 g/kg，低水平 0.5～0.75 g/kg），北方土石山区横坡垄作坡面 TN 含量处于中低水平。因此，北方土石山区棕壤分布地带需增施氮肥。
-
变异系数可直观反映样本变异大小。垄面 TN、 NO3--N 和 NH4 +-N 含量的变异系数为 17.64%～33.60%，略高于垄沟（18.10%～28.29%），两个特征点位土壤氮的变异均属于中等变异水平（10% ≤ CV ≤ 50% 或相对极差 >1.00）。进一步分析发现，垄面 NO3--N 的变异程度最大、TN 次之、NH4 +-N 最小；垄沟土壤氮的变异为 NH4 +-N>NO3--N>TN。对于整个坡面，土壤氮的变异隶属于中等变异，且变异程度表现为 NO3--N>TN>NH4 +-N。以上分析说明，横坡垄作坡面 TN、NO3--N 和 NH4 +-N 均具有一定的空间变异性，且变异程度主要取决于垄面这一特征点位。
-
2.2 横坡垄作坡面土壤氮的空间分布
-
垄面、垄沟内及整个坡面表层 0～20 cm 土层土壤 TN、NO3--N 和 NH4 +-N 含量的 Shapiro-Wilk 检验（表1）表明，横坡垄作坡面土壤氮含量呈正态分布。因此，可基于有效样点垄面和垄沟 TN、 NO3--N 和 NH4 +-N 含量，绘制土壤氮的空间分布等值线图和 3D 曲面图（图2），以此探析横坡垄作坡面土壤氮的空间变异特征。横坡垄作坡面 TN、 NO3--N 和 NH4 +-N 的平均含量分别为 0.65 g/kg、 4.90 mg/kg 和 9.47 mg/kg。虽然沂蒙山区棕壤土质疏松多孔，易使土壤中的部分铵态离子氧化成硝态离子，但 NO3--N 易随径流迁移且可被入渗水大量淋失，从而导致 NH4 +-N 含量高于 NO3--N。空间分布上，南北方向（垂直距离）TN、NO3--N 和 NH4 +-N 含量沿坡面呈不同的变化趋势，且 TN 的变异最大 （18.92%），其次为 NH4 +-N（10.51%） 和 NO3--N （10.09%）。东西向（水平距离），TN、NO3--N 和 NH4 +-N 含量表现为自东向西逐渐减少的趋势，变异系数分别为 8.35%～11.85%、4.20%～26.75% 和 3.82%～21.02%，且在坡脚部位土壤氮东西向的变异最大。可见，南北向、东西向土壤氮的空间分布均存在异质性，且坡面西北部应重点加强氮肥的管理。
-
2.3 横坡垄作措施下坡面土壤氮的空间分布影响因素
-
横坡垄作系统沟垄相间，垄沟内积水和泥沙沉积的发生，可导致垄面和垄沟间氮含量的差异，且土壤氮随垄面侵蚀泥沙、垄沟内沉积泥沙的迁移会造成其沿坡位分布的变异。为进一步探析土壤氮的空间分布特征，分析了坡位和垄面、垄沟等特征点位对表层 0～20 cm 土层土壤 TN、NO3--N 和 NH4 +-N 空间分布异质性的影响。
-
2.3.1 坡位对土壤氮空间分布的影响
-
垄面特征点位处，TN 含量表现为坡下 >坡脚 >坡中 >坡上（图3），但坡下和坡脚间无显著性差异。坡下、坡脚和坡中的 TN 含量较坡上分别显著增加 111.81%、75.06% 和 55.15%。垄沟、坡脚 TN 含量最大，但与坡下无显著性差异，且坡中与坡上之间差异不显著；坡脚 TN 含量较坡上显著增加 37.70%。可见，垄面 TN 含量沿坡位的变化趋势与垄沟特征点位相同，但垄面特征点位处 TN 含量的变异更为明显。
-
垄面特征点位处，4 个坡位间 NO3--N 和 NH4 +-N 含量均无显著性差异。垄沟处，坡上、坡中、坡下间的 NO3--N 含量差异不显著，但坡脚 NO3--N 含量较其他 3 个坡位显著增加 68.60%～83.85%； NH4 +-N 含量坡脚与坡上、坡下与坡中之间均差异不显著，而坡脚 NH4 +-N 含量较坡中显著增加 66.30%。可见，NO3--N 和 NH4 +-N 含量沿坡位的变异在垄沟特征点位体现的更为明显。
-
2.3.2 垄沟系统特征点位对土壤氮空间分布的影响
-
垄面与垄沟间 TN 含量的差异在坡上、坡下体现的最明显，但在坡中、坡脚位置二者之间差异不显著。垄面坡下 TN 含量较垄沟相应位置显著增加 21.56%（图3），但垄面坡上较垄沟相应位置显著减少 21.34%。垄面与垄沟间 NO3--N 含量在坡上、坡下和坡脚均存在显著差异，且在坡下体现的最为明显。垄面坡下、坡上和坡脚 NO3 +-N 含量较垄沟相应位置分别显著增加 38.16%、20.40% 和-37.47%。此外，垄面与垄沟间 NH4 +-N 含量在坡中、坡下、坡脚存在显著性差异，且在坡下差异最大。垄面坡下、坡中和坡脚 NH4 +-N 含量比垄沟相应位置分别显著增加 36.81%、26.27% 和-17.42%。可见，垄面与垄沟间土壤氮含量存在明显差异，且主要体现在坡下和坡上。
-
图2 表层土壤全氮、硝态氮、铵态氮的空间分布
-
图3 垄面、垄沟不同坡位处全氮、硝态氮和铵态氮含量
-
注：图中不同小写字母表示同一特征点位处坡位间氮含量在 P<0.05 水平差异显著；不同大写字母表示垄面和垄沟间氮含量在 P<0.05 水平差异显著。
-
2.4 横坡垄作坡面土壤氮的剖面分布特征
-
2.4.1 不同坡位处垄面、垄沟土壤氮的剖面分布
-
TN 含量的剖面分布特征与坡位密切相关，而垄面与垄沟间 TN 剖面分布的不同在坡上、坡中体现的最为明显。垄面坡上 TN 含量随土层深度增加呈减小-增加交替变化，而垄沟坡上则呈逐渐减小趋势（图4）；垄面坡中 TN 含量随土层深度增加波动减小，但垄沟坡中则为先逐渐增加而后减少趋势。垄面、垄沟坡下 TN 含量呈减小-增加交替变化，坡脚则为减小-增加-减小趋势。进一步分析发现，垄面 TN 含量在 10～15 cm 土层深度达到峰值，而垄沟则为 15～20 cm。垄面坡位间 TN 含量在 20～30 cm 处差异最大，坡中、坡下和坡脚较坡上分别增加 82.01%、65.11% 和 78.51%；垄沟处则为 10～15 cm，坡中、坡下和坡脚较坡上分别增加 104.02%、156.81% 和 114.34%。说明，横坡垄作坡面 TN 剖面分布受控于坡位、特征点位。
-
图4 垄面、垄沟不同土层深度内全氮、硝态氮和铵态氮的含量
-
同一坡位处，垄面 NO3--N 含量的剖面分布与垄沟明显不同。垄面 NO3--N 剖面分布呈逐渐减小-急剧增加-逐渐减小趋势，垄沟则呈逐渐减小趋势；垄面 NO3--N 含量在 15 cm 土层深度达到极值，垄沟则在 0～5 cm。垄面坡位间 NO3--N 含量在 20～30 cm 土层深度差异最大，此土层深度内坡中、坡下和坡脚的 NO3--N 含量较坡上分别增加 99.18%、5.71% 和 146.94%。垄沟坡位间的差异则主要体现在 10～15 cm，此处坡中、坡下和坡脚的 NO3--N 含量较坡上分别增加 11.75%、109.46% 和 0.60%。可见，NO3--N 剖面分布特征主要受控于特征点位。
-
垄面与垄沟 NH4 +-N 剖面分布的不同主要体现在坡下、坡脚。垄面坡下、坡脚 NH4 +-N 剖面分布均呈减小-增加-减小趋势；垄沟坡下为减小-增加，坡脚为逐渐减小趋势。但是坡上和坡中位置的垄面、垄沟内 NH4 +-N 含量均为减小-增加波动变化。同时，垄面、垄沟 NH4 +-N 含量在 10～15 cm 最大。垄面坡位间在 5～10 cm 处差异最明显，此处坡中、坡下和坡脚的 NH4 +-N 含量较坡上分别增加-22.26%、-19.22% 和 44.28%。垄沟坡位间在 10～15 cm 差异最大，该深度范围内坡中、坡下和坡脚的 NH4 +-N 含量较坡上分别增加 35.64%、25.12% 和 22.67%。说明，特征点位是影响 NH4 +-N 剖面分布的主要因素。
-
2.4.2 土壤氮剖面分布拟合
-
基于线性、指数、幂函数等对 TN、NO3--N、 NH4 +-N 含量的剖面分布进行拟合，发现垄面和垄沟土壤氮随土层深度增加的变化均符合指数函数关系（表2）。指数方程中垄面 TN、NO3--N 土层深度的指数系数均大于垄沟相应系数，但 NH4 +-N 呈相反趋势。表明垄沟 TN、NO3--N 含量随土层深度的递减速度较垄面更快，而垄面 NH4 +-N 含量随土层深度变化的降低幅度更大。
-
注：* 表示在 0.05 水平上显著，** 表示在 0.01 水平上极显著。
-
3 讨论
-
3.1 横坡垄作坡面土壤氮的空间变异性
-
土壤氮是土壤养分的主要组成成分之一，其空间异质性的研究是进行田间养分有效管理、精准农业合理实施的基础。坡面微地形很大程度上决定着土壤侵蚀的强弱、土壤水热分配，因而可直接影响土壤氮的空间分布[3，26]。横坡垄作作为被广泛采用的保护性耕作措施之一，坡面沟垄相间，使地表呈现有规则的高低起伏、凹凸不平的微地貌。同时，垄沟内雨水的汇聚和泥沙沉积，使得土壤氮在垄面、垄沟特征点位发生重新分配，进而改变土壤氮空间分布格局。
-
横坡垄作坡面 TN、NO3--N 和 NH4 +-N 含量自东向西逐渐降低（图2），呈现中等变异（变异系数 17.64%～33.60%）。这与在黄土高原地区农田、英国牧草地等小尺度上开展的坡面土壤氮空间变异的研究结果一致[27-28]，即小尺度上土壤氮空间变异性属于弱到中等变异。横坡垄作坡面土壤氮沿东西向出现一定空间异质性，主要与此坡面的地形和侵蚀形态密切相关。坡面东部的坡度（5.5°～16.0°） 大于西部（4.0°～7.2°），垄沟长达 46 m，导致垄沟内汇集的雨水不断流入西部。同时，野外调查发现，坡面西部部分垄沟内出现细沟，可使西部土壤氮随径流或泥沙输移进入坡面边缘的沟道，植物可直接吸收利用 NO3--N 和 NH4 +-N。可见，垄沟长度的增加会提高细沟发生几率[29]。因此，横坡垄作措施下应合理设置垄沟长度。同时，坡面出现不规则地形情况下，尤其当坡面水平方向起伏较大时，可采用斜坡垄作[30]、垄沟植草[31]等保护性措施，从而减小垄向坡度，削减垄沟内径流侵蚀量[32]。
-
以往研究认为，北方土石山区果园地总氮含量呈现坡上 >坡中 >坡下，撂荒地则为坡下 >坡上 >坡中[6]；黄土高原地区果园地土壤氮含量沿坡面呈逐渐降低趋势[7]；毛乌素地区人工林地土壤氮含量从大到小依次为坡顶、坡中、坡下[33]，而西南紫色土区坡耕地、人工林、荒草地的土壤氮含量均表现为坡下 >坡中＞坡上[34]。本研究中，横坡垄作坡面土壤氮沿坡位的变化与垄面、垄沟等特征点位相关（图3）。可见，土壤氮的空间分布受控于土地利用方式。同时，垄面与垄沟间土壤氮含量空间分布的差异主要反映在侵蚀强度较低的坡上、坡下和泥沙易沉积的坡脚部位（图3），即垄面坡上、坡脚 TN 含量较垄沟减少 21.34% 和 17.42%～37.47%，但坡下较垄沟增加 21.56%～38.16%。这是因为坡上垄面土壤氮易随径流、泥沙发生迁移，且雨水的汇聚将部分土壤氮携带进入垄沟，从而增加了垄沟内土壤氮的含量。坡下侵蚀强度低于坡上、坡中[35-36，23]，且坡下地势平缓，使雨水极易在垄沟内发生汇聚，加之坡上部径流的汇入，从而增加了入渗深度，减少了土壤水的无效蒸发[37]，进而提高了垄沟内土壤含水量。这可加速垄沟内 NO3--N 和 NH4 +-N 的淋溶损失[38-39]、提高土壤氮素的矿化速率、促进土壤氮的侧向迁移，从而降低了垄沟内土壤氮的含量。坡脚部位易发生泥沙沉积，尤其在垄沟部位，且沉积泥沙中土壤团聚体粒径较小，易造成土壤氮的富集[4，40-41]，进而提高了垄沟内土壤氮含量。因此，可选择在侵蚀部位的垄面、泥沙沉积区垄沟内种植作物，从而提高横坡垄作措施下氮肥利用效率、增加作物产量。
-
3.2 横坡垄作坡面土壤氮的垂直分布
-
土壤氮在土壤剖面上的盈余分布，直接影响氮向地下水的迁移[42]。土壤氮输入和输出的平衡决定了其剖面分布盈余。动植物残体的分解、氮沉降、根系分泌物的分解是土壤氮的主要来源，而输出量主要涉及有机物质的分解与土壤侵蚀[43]。所以，一般认为，平坡坡面土壤氮主要集中于土壤表层，土壤氮含量随土层深度增加呈减少趋势[13-14，44]。本文研究中横坡垄作坡面垄面、垄沟土壤氮随土层深度的增加均呈波动变化，尤其在垄面 0～30 cm 和垄沟 0～20 cm（图4），这可能与花生根系分布一般不超过 30 cm 土层深度有关。说明土壤氮的剖面分布与耕作方式密切相关。进一步分析发现，垄面与垄沟间土壤氮的剖面分布存在明显差异。土壤氮含量在 10～15 cm 土层深度出现峰值，而垄沟土壤氮含量在 0～5 和 15～20 cm 出现极值。但是垄面与垄沟内土壤氮的来源相同，说明二者的差异主要与土壤氮的输出方式有关。本研究中横坡垄作措施中垄高约为 15 cm，表明垄沟内的土壤氮可随土壤入渗水侧向淋溶至垄面 10～15 cm 土层深度。同时，垄面表层土壤氮可随雨水的汇集、径流而迁移至邻近垄沟。然而，垄沟土壤氮含量随土层深度增加而减小的趋势较垄面更为明显（表2）。由此可见，垄沟内土壤氮侧向迁移于垄面的含量大于垄面土壤氮的径流流失量。由于本研究中垄间距较小（10 cm），导致土壤氮侧向淋溶范围有限。可见，增大垄间距在一定程度上可提高土壤氮的利用效率，减少氮流失，进而提高作物产量。
-
4 结论
-
基于样带法采样，利用经典统计、地统计学方法，解析了横坡垄作坡面 TN、NO3--N 和 NH4 +-N 的空间变异性及其影响因素。主要结论如下：
-
（1）北方土石山区土壤氮含量属于中低水平， TN、NO3--N 和 NH4 +-N 变异均处于中等变异水平 （17.64%～33.60%），但垄面土壤氮的变异高于垄沟。
-
（2）TN 含量表现为坡下 >坡脚 >坡中 >坡上，且在垄面其含量沿坡位的变异更明显；垄沟 NO3--N、NH4 +-N 含量沿坡位的变化较垄面明显，表现为坡脚较坡上、坡中、坡下显著增加 66.30%～83.85%。
-
（3）垄面与垄沟间土壤氮含量的差异主要反映在坡上、坡下和坡脚。其中，垄面坡下土壤氮含量较垄沟增加 21.56%～38.16%，坡上、坡脚相应的分别减少 21.34% 和 17.42%～37.47%。
-
（4）TN 含量的剖面分布受坡位和特征点位的影响，而 NO3--N、NH4 +-N 的剖面分布主要受控于特征点位，但 TN、NO3--N、NH4 +-N 的剖面分布符合指数函数关系。
-
参考文献
-
[1] Jin Y H，Zhou D W，Jiang S C．Comparison of soil water content and corn yield in furrow and conventional ridge sown systems in a semiarid region of China［J］．Agricultural Water Management，2010，97（2）：326-332．
-
[2] 郭乾坤，刘宝元，朱少波，等．中国主要水土保持耕作措施因子［J］．中国水土保持，2013，10：22-26．
-
[3] 宋轩，李立东，寇长林，等．黄水河小流域土壤养分分布及其与地形的关系［J］．应用生态学报，2011，22（12）：3163-3168．
-
[4] 魏孝荣，邵明安．黄土高原沟壑区小流域坡地土壤养分分布特征［J］．生态学报，2007，27（2）：603-612．
-
[5] Gouri S B，Pravat K S，Rabindranath C．Assessment of spatial variability of soil properties using geostatistical approach of lateritic soil（West Bengal，India）［J］．Annals of Agrarian Science，2018，16：436-443．
-
[6] 冯雪瑾，张志华，杨喜田，等．太行山低山丘陵区人工林表层土壤有机碳和全氮分布特征［J］．应用生态学报，2019，30（2）：511-517．
-
[7] 陈伟，王红阳，王志坚，等．黄土丘陵区香梨园土壤水分、养分分布特征及其与产量的关系［J］．应用生态学报，2021，32（9）：3159-3166．
-
[8] 郭琦，卜崇峰，李宜坪，等．区域尺度生物结皮下伏土壤养分的空间分布特征-以毛乌素沙地为例［J］．土壤学报，2022，59（3）：699-708．
-
[9] 杜可，王乐，张淑香，等．黑土区县域土壤养分空间分布特征及其影响因子［J］．植物营养与肥料学报，2018，24（6）：1465-1474．
-
[10] 梁月明，苏以荣，何寻阳，等．喀斯特地区不同坡位条件下优势灌木根际与非根际土壤养分与pH的分布特征［J］．中国岩溶，2018，37（1）：53-58．
-
[11] 张明礼，杨浩，徐从安，等．沂蒙山区土壤侵蚀的 137Cs 示踪法初步研究［J］．土壤学报，2010，47（3）：408-414．
-
[12] 李龙，姚云峰，秦富仓，等．黄花甸子流域土壤全氮含量空间分布及其影响因素［J］．应用生态学报，2015，26（5）：1306-1312．
-
[13] 温延臣，李燕青，袁亮，等．长期不同施肥制度土壤肥力特征综合评价方法［J］．农业工程学报，2015，31（7）：91-99．
-
[14] 张燕江，王俊鹏，王瑜，等．农牧交错带典型区土壤氮磷空间分布特征及其影响因素［J］．环境科学，2021，42（6）：3010-3017．
-
[15] 张伟，陈洪松，王克林，等．典型喀斯特峰丛洼地坡面土壤养分空间变异性研究［J］．农业工程学报，2008，24（1）：68-73．
-
[16] 刘旭阳，陈晓旋，陈优阳，等．福州不同农田土地利用类型土壤碳氮磷生态化学计量学特征［J］．水土保持学报，2019，33（6）：348-355．
-
[17] Wang J，Fu B，Qiu Y，et al．Soil nutrients in relation to land use and landscape position in the semi-arid small catchment on the loess plateau in China［J］．Journal of Arid Environments，2001，48：537-550．
-
[18] Jia H Y，Lei A L，Lei J S，et al．Effects of hydrological processes on nitrogen loss in purple soil［J］．Agricultural Water Management，2007，89：89-97．
-
[19] 于兴修，马骞，刘前进，等．横坡与顺坡垄作径流氮磷输出及其富营养化风险对比研究［J］．环境科学，2011，32（2）：428-436．
-
[20] 闫建梅，何丙辉，田太强．不同施肥与耕作对紫色土坡耕地土壤侵蚀及氮素流失的影响［J］． 中国农业科学，2014，20：4027-4035．
-
[21] Cheng Y T，Li P，Xu G E，et al．The effect of soil water content and erodibility on losses of available nitrogen and phosphorus in simulated freeze-thaw conditions［J］．Catena，2018，166：21-33．
-
[22] 和继军，蔡强国，王学强．北方土石山区坡耕地水土保持措施的空间有效配置［J］．地理研究，2010，29（6）：1017-1026．
-
[23] 李秀彬，马志尊，姚孝友，等．北方土石山区水土流失现状与综合治理对策［J］．中国水土保持科学，2008（1）：9-15．
-
[24] 刘前进，于兴修．北方土石山区土壤侵蚀强度垂直景观格局—以沂蒙山区为例［J］．地理研究，2010，29（8）：1471-1483．
-
[25] 潘成忠，上官周平．黄土半干旱区坡地土壤水分、养分及生产力空间变异［J］．应用生态学报，2004，15（11）：2061-2066．
-
[26] 王珂，沈掌泉，John S B，等．精确农业田间土壤空间变异与采样方式研究［J］．农业工程学报，2001（2）：33-36．
-
[27] 王婕，牛文全，张文倩，等．农田表层土壤养分空间变异特性研究［J］．农业工程学报，2020，36（15）：37-46．
-
[28] Liu Q J，An J，Zhang G H，et al．The effect of row grade and length on soil erosion from concentrated flow in furrows of contouring ridge systems［J］．Soil &Tillage Research，2016，160：92-100．
-
[29] 宁静，杨子，姜涛，等．东北黑土区不同垄向耕地沟蚀与地形耦合规律［J］．水土保持研究，2016，23（3）：29-36．
-
[30] 李建华，于兴修，刘前进，等．沂蒙山区坡耕花生地垄间植草对土壤理化性质的影响［J］．水土保持学报，2012，26（5）：108-112，117．
-
[31] 赵玉明，姜洪涛，王世界．东北黑土区不同地形垄向坡度的相对水土保持效益研究［J］．水土保持研究，2014，21（5）：28-34．
-
[32] 刘鑫，满秀玲，陈立明，等．坡位对小叶杨人工林生长及土壤养分空间差异的影响［J］．水土保持学报，2007（5）：76-81．
-
[33] 高雪松，邓良基，张世熔．不同利用方式与坡位土壤物理性质及养分特征分析［J］．水土保持学报，2005（2）：53-56，60-79．
-
[34] 张合兵，聂小军，程静霞．137Cs 示踪采煤沉陷坡土壤侵蚀及其对土壤养分的影响［J］．农业工程学报，2015，31（4）：137-143．
-
[35] An J，Zheng F L，Wang B．Using 137Cs technique to investigate the spatial distribution of erosion and deposition regimes for a small catchment in the black soil region，Northeast China［J］． Catena，2014，123：243-251．
-
[36] Liu T N，Chen J Z，Wang Z Y，et al．Ridge and furrow planting pattern optimizes canopy structure of summer maize and obtains higher grain yield［J］．Field Crops Research，2018，219：242-249．
-
[37] 姜翠玲，夏自强，崔广柏．土壤含水量与氮化合物迁移转化的相关性分析［J］．河海大学学报（自然科学版），2003，31（3）：241-245．
-
[38] 王晓维，徐健程，龙昌智，等．施氮量和土壤含水量对黑麦草还田红壤氮素矿化的影响［J］．植物营养与肥料学报，2018，24（2）：365-374．
-
[39] 魏晗梅，郑粉莉，冯志珍，等．薄层黑土区流域尺度土壤养分对侵蚀-沉积的响应［J］．水土保持学报，2021，35（4）：49-54．
-
[40] Lu J，Zheng F L，Li G F，et al．The effects of raindrop impact and runoff detachment on hillslope soil erosion and soil aggregate loss in the Mollisol region of Northeast China［J］．Soil and Tillage Research，2016，161：79-85．
-
[41] 邹刚华，赵凤亮，单颖．典型红壤区不同土地利用方式下氮素垂直分布及其影响因素［J］．生态与农村环境学报，2019，35（5）：644-650．
-
[42] 林英华，汪来发，田晓堃，等．三峡库区杉木马尾松混交林土壤 C、N 空间特征［J］．生态学报，2011，31（23）：326-334．
-
[43] 李亚翠，季宏兵，张涛，等．石林地区典型土壤剖面有机碳、氮含量分布特征研究［J］．地球与环境，2018，46（5）：437-443．
-
摘要
横坡垄作是坡耕地广泛采用的水土保持措施，垄沟内雨水的汇聚和泥沙沉积，使得土壤氮的空间分布特征明显有别于传统耕作坡面。以北方土石山区横坡垄作坡面为研究对象，样带法采集土壤样品，基于经典统计、地统计学、Kriging 插值等方法，分析垄面、垄沟等特征点位全氮（TN）、硝态氮（NO3 - -N）、铵态氮（NH4 + -N） 含量沿坡位变化及其剖面分布，明晰横坡垄作坡面土壤氮的空间分布特征。结果表明：北方土石山区土壤氮含量属于中低水平，变异系数介于 17.64% ～ 33.60%，但垄面土壤氮的变异水平高于垄沟。TN 含量沿坡位的变异最为显著，表现为坡下 > 坡脚 > 坡中 > 坡上；仅垄沟 NO3 - -N、NH4 + -N 含量沿坡位存在显著变化，坡脚 NO3 - -N、 NH4 + -N 含量较其他坡位显著增加 66.30% ～ 83.85%。垄面与垄沟间氮含量的差异在坡上、坡下和坡脚体现的较为明显；与垄沟相比，垄面坡上 TN 含量减少 21.34%，坡脚 NO3 - -N 和 NH4 + -N 含量减少 17.42% ～ 37.47%，坡下土壤氮含量增加 21.56% ～ 38.16%。TN 剖面分布与坡位和特征点位有关，NO3 - -N、NH4 + -N 的剖面分布主要受控于特征点位，但氮含量剖面分布符合指数函数关系。综上，横坡垄作坡面土壤氮空间异质性明显，且主要受控于特征点位和坡位。研究结果可为提高横坡垄作措施下氮肥利用效率、精准农业的顺利推广提供技术支撑。
Abstract
Contour ridging is a widely used soil conservation practice for slope farmland. Rainwater from side slopes ponding and deposition sediment in the furrow areas could induce the spatial distribution of soil nitrogen on slope within contour ridge system differ from that on traditional slope. Based on line sampling method，this study analyzed the contents of total nitrogen（TN），NO3 - -N and NH4 + -N and their profiles distribution in ridges and furrows along with slope positions in order to determine the spatial distribution of soil nitrogen within contour ridge system in rocky mountainous areas of Northern China by the classical statistical and geostatistical methods. Results showed that soil nitrogen content in topsoil of rocky mountainous areas of Northern China was at the middle level. The variation coefficients of soil nitrogen ranged from 17.64% to 33.60%， and the variation in ridges was greater than that in furrows. The variation of TN content was the most significant with the change of slope positons，and was in order of downslope>toe slope>middle slope>upslope. However，the content of NO3 - -N and NH4 + -N in furrows at toe slope statistically increased by 66.30%～83.85% relative to other three slope positions，and no significant difference in ridges was observed among slope positions. In addition，the difference in soil nitrogen contents between ridges and furrows mainly reflected at upslope，down slope，and toe slope. Compared to furrows，TN content in ridges at upslope decreased by 21.34%，and NO3 - -N and NH4 + -N in ridges at toe slope decreased by 17.42%～37.47%， but soil nitrogen content in ridges at down slope increased by 21.56%～38.16%. Furthermore，the profile distribution of TN was related to slope position and feature point，but NO3 - -N and NH4 + -N profiles distribution were mainly controlled by feature points. Meanwhile，exponent function fitted the relationship between soil nitrogen content and soil depth. Thus， the pronounced heterogeneity existed for the spatial distribution of soil nitrogen within contour ridge system，which was determined by feature points and slope positions. These findings supply technological guidance for the improvement of nitrogen use efficiency within contour ridge system and the smooth promotion of precision agriculture.
Keywords
soil nitrogen ； ridge and furrow ； slope position ； spatial distribution ； contour ridge system