土壤有机碳(SOC)和全氮(STN)是土壤肥力和土壤质量的重要指标，是维持农业土壤生产力，保持土壤可持续利用的关键因子^[1]。我国是世界苹果第一生产大国，栽培面积和产量分别占世界的42%和54%^[2]。我国苹果园分布跨越不同气候带，在内陆和滨海形成了各具特色的主产区。然而，我国苹果园总体上管理粗放，片面追求产量，不够重视均衡施肥^[3]，导致了土壤有机质含量低、氮磷含量过高等问题^[4-6]。因此，探究苹果园土壤碳氮的时空分布与变化特征对指导果园土壤管理，促进果业可持续发展具有重要意义。

前人对我国不同区域苹果园碳氮水平及分布做了大量研究工作。Ge等^[7]研究指出，我国苹果园SOC和STN含量普遍偏低，不同区域差异较大。黄土高原苹果园SOC和STN具有中等空间相关性^[8]，并随着种植年限呈先增后降的变化趋势^[9-10]。然而，我国苹果园的气候条件、土壤类型和管理方式差异很大，果园土壤碳氮的时空分布特征可能具有区域性，滨海苹果园土壤碳氮空间分布及动态变化规律尚未系统研究，其与内陆果园的异同尚缺乏了解。为此，本研究以地处黄渤海之滨的胶东半岛苹果园为研究区，基于地统计学方法分析SOC、STN含量和密度(SOCD、STND)及碳氮比(C/N)的空间异质性，揭示SOC、STN含量和密度及C/N比值随种植年限的变化特征，探讨影响胶东苹果园SOC、STN含量和密度及C/N分布的影响因素，进而通过与内陆地区苹果园的比较，进一步阐明滨海果园碳氮时空分布与变化特征，为胶东滨海苹果园乃至全国果园土壤可持续管理提供基础数据和科学依据。

胶东半岛地处东经119°45’-122°75’，北纬35°50’-38°85’之间(图 1)，包括青岛、烟台和威海三市，总面积30 825 km²，人口约1900万。胶东半岛属暖温带半湿润季风气候，四季分明，光照充足，年降雨量约为650~850 mm，年均温约12℃，平均无霜期约210 d。胶东半岛属丘陵地貌，海拔介于0~1 132 m之间，地带性土壤类型为棕壤。胶东苹果种植面积约为22.85万km²，主要分布在平原和相对缓平的丘陵之上。

图 1 研究区位置与土壤采样点分布 Fig. 1 Map of the study area and distribution of soil sampling sites

样品采集于2014年4月至6月，共采集苹果园土壤488件(图 1)。果园种植年限信息通过咨询果农获取。采样深度为0~20 cm，在每个样点以10 m为半径，采集5个子样，充分混合后，用四分法取出1 kg的样品装于聚乙烯塑料袋，注明编号，带回实验室，自然风干。采样中心位置采用Garmin GPS 60CSx记录。

样品风干后，剔除动植物残体、石块及可见侵入体等，用陶瓷研钵研磨至分别完全通过2 mm(10目)，和0.149 mm(100目)的尼龙筛。SOC和STN含量分析采用100目土样，碳氮分析仪(Vario MAX CN Macro Elemental Analyzer)测定。对石灰性土壤样品，分析前先采用0.5 molL^-1盐酸去除碳酸盐，并用超纯水洗涤，冷冻干燥保存。土壤pH分析采用10目土样，土水比为1︰2.5，用pH计(Mettler Toledo Five Easy Plus FE20)测定。具体分析方法参考文献[11]。分析时，设置10%样品重复，以检验分析质量。

土壤有机碳密度(SOCD，kg m^-2)、全氮密度(STND，kg m^-2)按照如下公式^[12]进行计算：

SOCD=SOC×BD×H×(1-δ)/100

STND=STN×BD×H×(1-δ)/100

式中，BD为土壤容重(g cm^-3)，采用Song等^[13]建立的回归模型进行估算：BD=1.377×e^{(-0.004 8×SOC)}。H为土层厚度(20cm)，δ为大于2mm砾石的体积比(%)，依据《烟台市土壤》^[14]进行估算。

经典统计分析、相关分析和单因素方差分析利用IBM SPSS 20软件完成。在SOC、STN和C/N随年限的动态变化分析中，对各年限数据取平均值后再进行制图和拟合。运用地统计学手段进行SOC、SOCD、STN、STND和C/N空间异质性分析，半方差函数的拟合在GS+ 9.0软件上进行。将GS+中获得的最优拟合参数导入ArcGIS 10.2软件的地统计模块(Geostatistical Analyst)中，运用普通克里格插值法(Ordinary Kriging)获得SOC、STN和C/N的空间分布图。海拔、坡度、坡向等均基于30m×30m分辨率的DEM数据，采用ArcGIS 10.2软件计算得到。

由表 1可见，胶东苹果园SOC、STN的平均含量为10.78~1.42 g kg^-1，平均密度为2.81、0.37k gm^-2，平均C/N为7.70。与我国其他典型苹果园相比，胶东苹果园SOC、STN的含量和密度处于相对较高水平，但C/N偏低，其含量和密度均高于辽西和渭北果园，但低于京郊和伊犁果园。胶东苹果园处于温带湿润丘陵地区，地势起伏幅度较小，土壤质地较疏松；受渤黄海海洋性气候影响，水热相对充足，但蒸发量约是降雨量的2倍，水土流失较弱，果树生长较快，总体上有利于土壤碳氮的积累。目前，SOC含量处于较高水平；就氮而言，由于常年氮肥施用量大，而有机肥使用少，STN的积累速率快于SOC。葛顺峰等^[4]认为土壤C/N比处于21~23有利于提高氮素的利用率，并促进苹果植株的生长。胶东果园C/N远低于此范围，说明该区土壤碳氮比例失调，可供微生物利用的碳源较少，微生物活性可能较低，可能会影响氮素的转化。

表 1 Table 1

表 1 胶东苹果园SOC、STN含量和密度及C/N的统计特征及其与我国内陆主要果园的比较 Table 1 Statistics of SOC, SOCD, STN, STND and C/N in the orchardsof Jiaodong Peninsula and their comparison with those in the typical inland orchards of China 样地 Sample site SOC（g kg-1） SOCD（kg m-2） STN（g kg-1） STND（kg m-2） C/N 胶东苹果园① 10.78±3.38 2.81±0.82 1.42±0.43 0.37±0.11 7.7±1.23 变异系数CV 0.31 0.29 0.30 0.29 0.16 辽西果园②[15] 7.77 1.97 1.45 0.37 5.36 京郊果园③[16] 11.84 3.46 1.53 0.45 7.74 渭北果园④[17] 6.14 1.54 0.72 0.18 8.53 伊犁果园⑤[18] 11.83 2.55 — — — ①Orchards in Jiaodong Peninsula; ②Orchards in West Liaoning; ③Orchards in suburbs of Beijing; ④Orchards in Weibei, Shaanxi; ⑤Orchard in Ili, Xinjiang

表 1 胶东苹果园SOC、STN含量和密度及C/N的统计特征及其与我国内陆主要果园的比较 Table 1 Statistics of SOC, SOCD, STN, STND and C/N in the orchardsof Jiaodong Peninsula and their comparison with those in the typical inland orchards of China

空间分析结果(表 2)显示，胶东果园SOC、SOCD、STN、STND和C/N均符合指数模型，块基比分别为0.432、0.340、0420、0.387和0.391，均介于0.25~0.75之间，属中等强度的空间异质性^[19]。这表明五者均受到结构性因素(如土壤类型、成土母质、地形条件等)和随机性因素(如肥料施用、翻耕灌溉、种植年限等)的共同影响^[19]。这与郭宏^[8]和张彬^[20]等在相同尺度下对典型黄土高原果园的分析结果一致，说明滨海果园与内陆果园SOC和STN具有相似的空间异质性，且控制空间异质性的因素也是一致的。胶东果园SOC、SOCD、STN、STND和C/N的空间分布如图 2所示，SOC和STN分布趋势一致，反映在两者的显著相关性上(r=0.800，p < 0.01)，同时暗示着C/N具有相对的稳定性^[21]。胶东果园SOC分布以10.00~12.00 g kg^-1区间为主，SOCD以2.40~3.00 k gm^-2为主，STN以1.20~1.60 g kg^-1为主，STND以0.30~0.42 k gm^-2为主，C/N以7.0~8.0区间为主。

表 2 Table 2

表 2 胶东苹果园SOC、STN含量和密度及C/N的最优半方差函数拟合参数 Table 2 Fitting parameters of the optimal semivariograms of SOC, SOCD, STN, STND and C/N in the orchards of Jiaodong Peninsula 项目 Item 块金值 Nugget 基台值 Sill 块基比 Nugget/Sill 变程 Range (m) R2 残差 RSS LnSOC 0.042 0.096 0.432 4200 0.618 1.22E-04 LnSOCD 0.030 0.087 0.340 4 200 0.499 1.05E-04 LnSTN 0.040 0.096 0.420 5400 0.705 7.85E-05 LnSTND 0.034 0.089 0.387 5 400 0.706 6.78E-05 C/N 0.593 1.517 0.391 6300 0.601 3.63E-02   注：胶东果园SOC、SOCD、STN和STND不满足正态分布，故采用对数变换后的数据进行半方差函数拟合Note: SOC, SOCD, STN and STND in the orchard soils of Jiaodong Peninsula were not in normal distribution, thus logarithmictransformed data were utilized in fitting of semivariogram

表 2 胶东苹果园SOC、STN含量和密度及C/N的最优半方差函数拟合参数 Table 2 Fitting parameters of the optimal semivariograms of SOC, SOCD, STN, STND and C/N in the orchards of Jiaodong Peninsula

图 2 胶东苹果园SOC、STN含量和密度及C/N的空间分布 Fig. 2 Spatial distribution of SOC, SOCD, STN, STND and C/N in the orchards of Jiaodong Peninsula

因胶东苹果园SOC、STN和SOCD、STND的动态变化规律基本一致，故本文仅以SOC和STN指示土壤碳氮的动态变化。胶东苹果园SOC、STN和C/N随种植年限的变化趋势如图 3所示。SOC含量随种植年限的增加呈现上升趋势，但上升的过程是有起伏的，可能与土壤类型及果园管理措施的差异有关。在种植约30年时SOC呈现下降趋势，此时间节点与苹果树的“枯果期”基本一致。我国果园多处在瘠薄土壤上，立地条件较差^[22]。且果园改建对土壤扰动强烈，使得有机质暴露，有利于微生物对有机质的矿化分解^[23]，因此果园利用初期SOC含量较低。随着经营年限的延长，以凋落物和有机肥形式输入的碳素，使微生物呼吸作用增强，在释放新增碳的同时，也释放出已贮存的有机质^[24]，因而SOC含量不断上升。SOC含量的升高促进了团聚体的形成，进而也增强了有机质的稳定性^[25]。对于树龄超过30年的果园土壤，处于退化阶段，果树细根的生物量减小^[26]，根系分泌物减少，果树的光合生理功能也变弱^[9]。此外，老果园常因产量低而疏于施肥^[27]，这些因素均可导致老龄果园SOC含量的下降。

图 3 胶东苹果园SOC、STN和C/N比值随种植年限的变化 Fig. 3 Variation of SOC, STN and C/N in the orchards of Jiaodong Peninsula with age

胶东果园STN含量随着年限的延长而呈上升趋势，但增长率不断减小。氮肥、有机肥的施用及凋落物输入是果园STN的主要来源，其含量是矿化与积累的平衡结果。果树生长需要吸收大量氮素，尤其在盛果期，但因施肥量大，使得STN含量并没有降低。在我国苹果生产中，氮肥利用率仅为30%~50%^[28]。盈余的氮素一部分被淋洗或挥发外，其余在土壤中渐渐积累，使得STN含量不断增加。随着年限延长，长期单一、过量施用氮肥导致土壤蓄积养分能力下降^[29]，因此STN增长率逐渐下降。此外，老果园施肥不受重视，果树归还量降低，导致STN含量不再显著提高。

SOC和STN积累并非同步，因此C/N变化规律与SOC和STN不同。随着年限的延长，土壤C/N呈先下降、后上升、再下降的趋势。幼龄果园因生产需要而氮肥投入大，氮素积累，因而C/N降低。随着年限进一步延长，SOC不断积累，土壤微生物数量会增加，氮素矿化率提高，因而C/N增加^[30]。而在枯果期，果园SOC含量开始下降，土壤微生物数量相应降低，氮素利用率下降，故而C/N再次下降^[30]。

甘卓亭等^[9]研究报道，内陆苹果园碳氮表现出随种植年限动态变化。然而，不同于胶东苹果园，可能由于受水分条件限制，内陆苹果园的生产周期一般较短，如渭北塬区20年限苹果园基本处于退化阶段，其碳氮含量也随之减少^[9-10]。黄土高原和北京昌平、通州等地区果园SOC也随种植年限变化^[31]，但不如胶东地区果园显著，这可能与不同区域果园管理措施、区域气候条件和土壤特征有关。胶东果园水热相对充足，植被生产力高，土壤微生物活动更旺盛，因此土壤的碳氮动态变化会强于内陆果园，其碳氮周转的速率会更快。

在现代果园管理措施中，常采用“清园”管理，导致地表凋落物残体很少。因而，施肥被认为是影响土壤碳氮含量的主要因素^[32]。不同地区苹果园施肥量及肥料种类差异较大^[33-34]，增大了土壤碳氮的空间变异性，使得果园土壤碳氮的空间分布和动态变化变得更加复杂。据报道，胶东果园施用化学氮肥水平在612.6 kghm^-2，局部地区高达1 000 kg hm^-2以上，而有机肥投入却明显不足^[33]。渭北果园氮肥用量在672 kg hm^-2，北京昌平果园氮肥用量约为1 000 kg hm^-2，同样有机肥投入比例较低^{[27, 34]}。这意味着不同地区果园的集约化程度不同，共同的是大量使用化肥，有机肥投入量过少^{[27, 33-34]}。这影响着不同区域果园碳氮的均衡性和时空特征。此外，灌溉、地膜覆盖、生草和翻耕等措施也可能在不同程度上影响果园土壤的碳氮含量水平及时空变化^[32]。

种植年限是衡量果园管理强度的重要指标，代表了果园的开垦历史、管理强度，因而会对SOC、STN含量和密度及C/N有显著影响。同时，不同果园因种植年限不同而使其SOC、STN含量和C/N不同，进而在空间上呈现出异质性。由图 4可见，2014年胶东果园典型土壤(以土属为单元)SOC和STN含量较1980年增长约1倍，说明长期施肥等管理措施提高了果园的碳氮水平。在不同种植年限之间，胶东苹果园碳氮水平也表现出显著的差异(图 2)，同样也印证了施肥等田间管理对碳氮含量及密度的影响。

图 4 胶东苹果园典型土壤类型SOC和STN含量在2014年与1980年之间的对比 Fig. 4 Comparison between the years of 1980 and 2014 in SOC andSTN content in the orchards of Jiaodong Peninsula 注：1麻砂质棕壤，2麻砂质棕壤性土，3中性粗骨土，4麻砂质中性粗骨土，5暗泥质中性粗骨土Note：1Hapli-Udic Argosols, 2Hapli-Udic Cambosols, 3~5Lithic Usti-Orthic Primosols

胶东果园SOC和STN含量在不同土属之间表现出显著性差异(p < 0.05)(图 5)，F值分别为2.06和1.85，说明土壤类型影响碳氮的时空分布。在土壤分类体系中，土属主要根据成土母质及区域水文控制的盐分类型等地方性因素划分，主要反映出母质和地形的影响。胶东果园土壤的成土母质主要包括残积坡积物、酸性岩坡积洪积物和洪积冲积物等类型，果园SOC和STN含量既继承于母质，又受母质控制的土壤性质(如质地)的影响，从而表现出与成土母质的关联性。

图 5 基于土属统计的胶东苹果园SOC、STN含量和C/N比值的分布 Fig. 5 Distribution of SOC, STN content and C/N in the orchards of Jiaodong Peninsula on a soil family basis 注：1麻砂质棕壤，2麻砂质棕壤性土，3灰泥质棕壤性土，4洪冲积壤质潮棕壤，5褐土，6淋溶褐土，7中性粗骨土，8麻砂质中性粗骨土，9暗泥质中性粗骨土，10钙质粗骨土，11砂质潮土，12壤质潮土，13砂壤质潮土，14黑姜土 Note:1Hapli-Udic Argosols, 2~3Hapli-Udic Cambosols, 4Mottlic Hapli-Udic Argosols, 5~6Hapli-Ustic Argosols, 7~9 Lithic Usti-Orthic Primosols, 10 Carbonatic Udi-Orthic Primosols, 11~13 Ochri-Aquic Cambosols, 14 Shajiang Calci-Aquic Vertosols

土壤分布具有地带性特征，而果园种植也受到地带性土壤类型的影响。胶东果园和辽西果园等滨海果园以棕壤为主，渭北果园以黄绵土为主，伊犁果园多以潮土为主。这几类土壤性质差别很大，即使在同等的管理强度下，肥力水平也各不相同。土壤类型也影响土壤有机质的稳定性，进而可影响碳氮矿化的速率^[35]。

气候在土壤碳氮的积累过程中扮演重要角色。一方面，气候影响果园植被的生产力和生物量，控制输入土壤的碳氮量；另一方面，气候通过土壤水分、温度等来影响微生物对碳氮的分解和转化，因此温度和降雨的综合作用决定了果园土壤碳氮分布的地带性特征^[32]。胶东苹果园与我国其他内陆果园存在着显著的气候条件差别。胶东苹果园深受渤黄海海洋性气候的影响，水热条件相对充足，碳氮的周转速度较快。胶东半岛自身由于海陆、陆地下垫面分布不均匀，造成了冷热水平和垂直分布的差异，使得滨临渤海果园比滨临黄海果园每年有100 mm降水量的差异，这可能会影响碳氮的时空分布特征，需要进一步研究。

胶东果园SOC、STN和C/N与地形因子的相关性见表 3。总体而言，SOC、STN和C/N比值与胶东地形因子有关。SOC和STN与海拔呈极显著正相关(p < 0.01)，相关系数分别为0.161和0.172；SOC和TN与坡度呈显著负相关(p < 0.05)。与张彬等^[20]对黄土高原果园的研究结果相比，胶东果园碳氮受地形因子的影响较弱，这主要与胶东果园丘陵、低山为主，地形起伏不大以及水土流失较弱有关。

表 3 Table 3

表 3 胶东苹果园SOC、STN和C/N与地形因子、土壤pH的相关性 Table 3 Pearson correlation of SOC, STN and C/N with terrain indices and soil pH 项目 Item 海拔 Elevation 坡度 Slope 坡向 Aspect pH SOC 0.161** -0.107* -0.045 0.102* STN 0.172** -0.103* 0.042 -0.124** C/N -0.020 0.003 -0.024 0.265**   注：*在0.05水平（双侧）上显著相关，**在0.01水平（双侧）上显著相关Note: *significant correlationat the 0.05 level, **significant correlation at the 0.01 level

表 3 胶东苹果园SOC、STN和C/N与地形因子、土壤pH的相关性 Table 3 Pearson correlation of SOC, STN and C/N with terrain indices and soil pH

胶东苹果园土壤pH平均为5.88，且随种植年限快速降低；相比于第二次土壤普查时的数值，胶东果园土壤严重酸化^[7]。胶东果园土壤SOC与pH存在显著正相关(r=0.102，p < 0.05)(表 3)，但土壤STN与pH呈极显著负相关关系(r=-0.124，p < 0.01)，可能与化学氮肥在硝化过程中产生大量H⁺和NO₃^-导致土壤pH降低有关^[34]。土壤C/N比值和pH呈极显著正相关(r=0.265，p < 0.01)。土壤pH的高低会影响土壤微生物的活动，因为大多数土壤微生物都对酸敏感^[36]。土壤的细菌和放线菌适宜生活在中性至微碱性的土壤环境中，当土壤pH过低时，它们的活性会受到严重影响，使得土壤有机质矿化速率下降^[37]。低pH可限制硝化微生物的生长，抑制土壤硝化作用，进而影响氮素的周转^[38]。与胶东滨海果园不同，黄土高原等内陆果园土壤pH多为碱性，碳酸盐对酸缓冲能力强，因而在现阶段，内陆苹果园土壤酸化对土壤碳氮含量及密度分布的影响可能不显著。

胶东滨海果园SOC和STN含量和密度较高，总体上果园土壤肥力质量较高。胶东果园SOC、STN含量和密度及C/N表现出中等强度的空间分布变异性，受到果树种植和管理年限的显著影响。在种植年限30年之内，SOC和STN含量持续增加。土壤类型、酸碱性、地形、施肥等自然和人为因素共同影响着SOC和STN含量和密度及其变化。滨海与内陆苹果园的土壤碳氮空间分布和动态变化存在异同性。