李海泉，邱开阳*，谢应忠，何毅，苗虎，刘王锁，2，黄业芸
（1.宁夏大学农学院，宁夏 银川 750021；
2.宁夏葡萄酒与防沙治沙职业技术学院，宁夏 银川 750199）
我国是世界上土地沙漠化较为严重的国家，沙漠化已成为制约我国经济社会发展的重要因素，我国在沙漠化治理方面形成了诸多技术［1］。草方格协同种植不同固沙灌木作为沙漠化防治过程中最为常见的工程措施，被称为草方格造林［2］。在实践中采用的固沙灌木，有时为一种，有时是两种灌木进行混交，具有一定的随意性。有研究表明，不同灌木对降水的再分配［3］对地上-地下生物量分配［4］会产生不同的影响。同时，不同固沙灌木的凋落物含量及分解速率也不相同［5］，进而，会对土壤养分的积累产生不同的影响。罗雅曦等［6］在宁夏风沙区不同人工固沙灌丛林的研究发现，相比其他灌木类型，油蒿林地对土壤养分的提高最为显著；
而张立欣等［7］对库布齐沙漠不同人工固沙灌木林的研究发现，几种灌木类型下油蒿+杨柴混交林对土壤养分的提高最为显著；
董雪等［8］对科尔沁沙地4种典型灌木灌丛下土壤碳、氮、磷化学计量特征研究表明，东北木蓼和小叶锦鸡儿更有利于土壤养分的积累。这些研究表明，不同固沙灌木对土壤养分的提高存在一定的差异。
沙地松软，沙粒之间排列疏松，孔隙较大，因此沙地是个良好的碳库。沙地经常掩埋各种枯落物而增加其碳储量。同时，沙地土壤易渗漏，可溶性的氮磷养分也经常由于渗漏或风蚀而流失。因此，沙地土壤的碳、氮、磷含量易于变化，而这方面的研究较少，尤其是沙地不同固沙灌木配置模式对土壤碳、氮、磷含量及其化学计量比的影响尚不明确。
生态化学计量学是通过结合多学科的基本原理，形成对生态系统C、N、P等化学元素平衡以及C、N、P等化学元素对生态系统交互作用影响研究的一种理论方法［9］，C、N、P元素之间互相作用是土壤系统中的重要过程［10］，同时土壤C、N、P 作为生态化学计量学领域研究的主要对象，土壤C、N、P生态化学计量比是评价土壤质量的重要指标之一［11］。近年来生态化学计量学方面的研究主要集中在农田［12］、森林［13］、草原［14］等陆地生态系统，相比而言沙漠地区的相关研究较少。许多研究表明土壤C∶N∶P的生态化学计量比受温度［15］和降水［16］等气候因素、海拔［17］和纬度［18］等地形因子的影响。但在沙地生态恢复过程中不同固沙灌木配置模式的影响尚不清楚。
宁夏是我国沙漠化最为严重的地区之一，草方格造林作为该地区主要的固沙措施，已有多年的应用。宁夏也是全国最早开始应用草方格耦合灌木进行固沙的地区，利用固沙灌木进行沙地恢复在该区域的应用也非常广泛。白芨滩国家级自然保护区是中国草方格固沙的典范。但是，草方格灌木固沙过程中种植不同灌木对土壤化学计量特征，尤其是C∶N∶P比值的影响尚不明确。因此，本研究以宁夏白芨滩国家级自然保护区使用最广泛的5种固沙灌木配置模式：柠条纯林、杨柴纯林、柠条×花棒混交、柠条×沙拐枣混交、花棒×沙拐枣混交为研究对象，研究不同固沙灌木配置模式对土壤C、N、P含量及其生态化学计量比的影响。5种固沙灌木中柠条作为深根性树种，耐寒、耐旱、耐瘠薄、抗高温，根系发达、具根瘤、保水性能好，这些特性使得柠条可以高度适应干旱的沙地，并且具有较强的再生能力，被沙埋之后，依然可以产生不定根，生长出新的枝条，有很好的固沙作用。花棒为沙生、耐旱、喜光树种，耐干旱，抗风蚀，喜适度沙压，生长迅速，根生长快，侧根发达，根幅大，易繁殖，种子萌发率高、枝条及根易萌生新芽，形成无性繁殖体，可平茬复壮，是防风固沙的优良树种。花棒全株几乎无小叶或小叶极少，完全以绿色的叶轴代替小叶进行光合作用。杨柴是多年生落叶半灌木，具有耐寒、耐旱、耐贫瘠、抗风沙的特点，适应性强。杨柴是多年生落叶半灌木，具有耐寒、耐旱、耐贫瘠、抗风沙的特点。柠条、花棒和杨柴都属于豆科固沙灌木，根系发达，有根瘤菌，能固定空气中的游离氮，增加土壤氮含量，从而使土壤肥力得到提高和改善。沙拐枣是蓼科沙拐枣属灌木，具有抗风蚀、耐沙埋、抗干旱、耐瘠薄等特点，能适应条件极端严酷的干旱荒漠区环境，是荒漠区典型的沙生植物。不同固沙灌木的生活习性不同，凋落物的数量和质量也存在不同的差异，柠条比较高大，枝叶茂盛，凋落物比较多，而花棒、杨柴和沙拐枣相对矮小一些，枝叶也较为稀疏，凋落物相对较少。
探究宁夏东部风沙区不同固沙灌木配置模式下土壤生态化学计量变化规律，寻求最佳的固沙灌木种植模式，为宁夏沙地恢复提供理论依据。
1.1 研究区概况
研究区位于宁夏灵武白芨滩国家级自然保护区，该地区位于宁夏东部风沙区，属于宁夏毛乌素沙地南缘，地 理 坐 标E 106°21′33″～106°37′00″，N 37°48′ 28″～38°20′12″，研究区气候类型属于典型的大陆性季风气候，降水的季节变化和年际变化均较大，多年平均降水量230.49 mm，降水量季节分布不均，以7-9月最多，占全年降水量的60%以上。该地区多年平均气温为10.24 ℃，≥10 ℃的时间段分布在4-10月，年平均无霜期为154 d。研究区位于我国内陆高原地区，属于多风地区，全年大风（17 m/s以上）日数为63 d；
年沙暴日数为35 d。土壤类型主要为灰钙土和风沙土。研究区主要固沙灌木为柠条 （Caragana korshinskii）、花棒（Corethrodendronscoparium）、沙拐枣（Calligonum mongolicum）、杨柴（C.fruticosumvar.mongolicum等，植被多数是沙生植被，如黑沙蒿（Artemisia ordosica）、沙木蓼（Atraphaxi sbracteata）、沙蓬（Agriophyllum squarrosum）、烛台虫实（Corispermum candelabrumIljin）等。研究区5种典型固沙灌林基本信息见表1。
表1 研究区5种典型固沙灌木样地基本信息Table 1 Basic information of five typical sand-fixing shrub plots in the study area
1.2 试验设计
于2019年7月在白芨滩国家级自然保护区内以2002年建立的固沙恢复区为研究对象，分别选取5种不同固沙灌木配置模式（表2）样地各3个，每个样地面积100 m×100 m，各样地之间具有相同的地形条件。在每个样地内设置1条50 m长的样线，在每条样线上，每10 m设置1个10 m×10 m的灌木样方，5种固沙灌木配置模式共设置75个样方。在每个样方内采用5点法分别采集0～10 cm、10～20 cm和20～30 cm的土样，把同一个样方内5个点采集的同一土层的样品均匀混合形成一个样品，装入自封袋并带回实验室。将土壤样品放置在阴凉干燥通风且无特殊气体和尘土污染的室内，经过自然风干后用镊子拣出枯枝落叶和植物根等杂质，然后采用四分法对土样进行过筛保存，过筛具体步骤：将风干土壤过2 mm筛后，取出一半继续过1 mm筛，之后再取出一半继续过0.15 mm筛，无法通过的土壤进一步研磨，直至全部土样通过0.15 mm筛。将通过0.15 mm筛的土壤样品充分混匀装入自封袋保存，本研究中土壤有机碳、全氮和全磷均采用过0.15 mm筛后的土样进行测定。每一种固沙灌木配置模式的土壤有机碳、全氮和全磷含量均通过测定15 个重复求其平均值。
表2 5种固沙灌木配置模式概况Table 2 Overview of 5 different sand-fixing shrub configuration patterns
1.3 土壤样品分析方法
土壤有机碳含量采用重铬酸钾容量法测定，全氮含量采用凯氏定氮法测定，全磷含量采用氢氧化钠熔融-钼锑抗比色法测定［19］。
1.4 数据处理与分析
分别使用土壤有机碳含量、全氮含量和全磷含量代表土壤的C、N、P含量。采用Excel 2016对数据进行录入和初步整理，剔除异常数据。利用SPSS 22.0进行单因素方差分析（One-Way ANOVA）和Duncan 多重检验，对各个土层下5种固沙灌木配置模式的土壤样品C、N、P 含量及C∶N、C∶P和N∶P生态化学计量比分别进行差异性检验，显著性水平为α=0.05；
运用Pearson Correlation分析土壤C、N、P含量与其化学计量比的相关性；
根据相关分析结果，进一步采用一元线性回归模型分析化学计量比与土壤C、N、P含量之间的关系；
利用Matlab软件进行作图，图中数据均为平均值±标准误。
2.1 不同固沙灌木配置模式对土壤C、N、P含量的影响
在0～10 cm土层，柠条×花棒混交与柠条×沙拐枣混交固沙灌木配置模式下土壤全氮含量显著高于其他3种固沙灌木配置模式（P＜0.05）（图2），土壤有机碳含量和土壤全磷含量在不同固沙灌木配置模式间差异均不显著（P＞0.05）（图1、3）。在10～20 cm土层，柠条和杨柴固沙灌木配置模式有机碳含量显著高于其他3种固沙灌木配置模式（P＜0.05）（图1），柠条×沙拐枣混交固沙灌木配置模式全氮含量显著低于其他4种固沙灌木配置模式（P＜0.05），而杨柴固沙灌木配置模式全氮含量显著高于其他4种固沙灌木配置模式（P＜0.05）（图2），5种固沙灌木配置模式全磷含量差异不显著（P＞0.05）（图3）。在20～30 cm土层，柠条×沙拐枣混交固沙灌木配置模式有机碳含量显著低于其他4种固沙灌木配置模式（P＜0.05）（图1），柠条固沙灌木配置模式全氮含量显著高于其他4种固沙灌木配置模式（P＜0.05），同时柠条×花棒混交固沙灌木配置模式全氮含量显著高于另外3种固沙灌木配置模式（P＜0.05）（图2），杨柴固沙灌木配置模式全磷含量显著低于其他4种固沙灌木配置模式（P＜0.05）（图3）。
图1 不同固沙灌木配置模式土壤有机碳含量Fig.1 Soil organic carbon content in surface and near surface layers of different sand-fixing shrub configurations
图2 不同固沙灌木配置模式土壤全氮含量Fig.2 Soil total nitrogen content of different sand-fixing shrub configurations
图3 不同固沙灌木配置模式土壤全磷含量Fig.3 Soil total phosphorus content of different sand-fixing shrub configurations
2.2 不同固沙灌木配置模式对土壤C、N、P生态化学计量比的影响
C∶N在0～10 cm土层，柠条×花棒混交和柠条×沙拐枣混交固沙灌木配置模式显著低于其他3种固沙灌木配置模式；
在10～20 cm土层，柠条和柠条×沙拐枣混交固沙灌木配置模式显著高于其他3种固沙灌木配置模式；
在20～30 cm土层中，花棒×沙拐枣混交固沙配置模式最大，显著高于其他4种固沙灌木配置模式，而柠条和柠条×花棒混交固沙灌木配置模式最小，显著低于其他3种固沙灌木配置模式（图4）。
图4 不同固沙灌木配置模式土壤碳氮比Fig.4 Soil C∶N ratios of different sand-fixing shrub configurations
C∶P在0～10 cm土层，柠条和花棒×沙拐枣混交固沙灌木配置模式显著低于其他3种固沙灌木配置模式（图5）；
在10～20 cm土层，柠条固沙灌木配置模式显著高于其他4种固沙灌木配置模式，柠条×沙拐枣混交固沙灌木配置模式显著低于其他4种草方格造林树种配置模式（图5）；
在20～30 cm土层中，柠条×沙拐枣混交固沙灌木配置模式显著低于其他4种固沙灌木配置模式（图5）。
图5 不同固沙灌木配置模式土壤碳磷比Fig.5 Soil C∶P ratios of different sand-fixing shrub configurations
N∶P在0～10 cm土层，柠条×花棒混交和柠条×沙拐枣混交固沙灌木配置模式显著高于其他3种固沙灌木配置模式（图6）；
在10～20 cm土层中，杨柴固沙灌木配置模式最高，显著高于其他4种固沙灌木配置模式，而柠条×沙拐枣混交固沙灌木配置模式最小，显著低于其他4种固沙灌木配置模式（图6）；
在20～30 cm土层，柠条固沙灌木配置模式显著高于其他4种固沙灌木配置模式，柠条×沙拐枣混交和花棒×沙拐枣混交固沙灌木配置模式显著低于其他3种固沙灌木配置模式（图6）。
图6 不同固沙灌木配置模式土壤氮磷比Fig.6 Soil N∶P ratios of different sand-fixing shrub configurations
2.3 不同固沙灌木配置模式下土壤C、N、P含量与其生态化学计量比之间的关系
柠条固沙灌木样地C∶P与C含量呈极显著正相关（P＜0.001），与P呈极显著负相关（P＜0.001），C∶N与N呈极显著负相关（P＜0.001），C∶P与N呈极显著正相关（P＜0.001）（表2）；
杨柴固沙灌木样地C∶N、C∶P均与C含量呈极显著正相关（P＜0.001），N∶P与N呈极显著正相关（P＜0.001），C∶N分别与N含量、P含量呈极显著正相关（P＜0.001）（表3）；
柠条×花棒混交固沙灌木样地C∶N与N含量呈显著负相关（P＜0.01），C∶P与C含量呈极显著正相关（P＜0.001），N∶P与N含量呈极显著正相关（P＜0.001）（表3）；
柠条×沙拐枣混交固沙灌木样地C∶N与N含量呈极显著负相关（P＜0.001），C∶P与C含量呈极显著正相关（P＜0.001），C∶P与P含量呈极显著负相关（P＜0.05），N∶P与N含量呈极显著正相关（P＜0.001）（表3）；
花棒×沙拐枣混交固沙灌木样地C∶N、C∶P均与C含量呈极显著正相关（P＜0.001），N∶P与C含量呈显著负相关（P＜0.01），C∶N与N含量呈显著负相关（P＜0.01），C∶P分别与N、P含量呈显著负相关（P＜0.01），N∶P与N含 量 呈 极 显 著 正 相 关（P＜0.001）（表3）。
表3 不同固沙灌木配置模式下土壤C、N、P含量与其化学计量比的相关性Table 3 Correlation between soil C，N，P contents and their stoichiometric ratios under different sand-fixing shrub configuration patterns
土壤C、N、P化学计量比与其含量之间存在线性关系。C∶N和C含量的拟合方程是y=16.363x+12.048（R2=0.139，P＜0.001），和N含量的拟合方程是y=-584.809x+69.558（R2=0.504，P＜0.001），拟合强度为N＞C（图7）。C∶P与C含量的拟合方程为y=5.477x+0.123（R2=757，P＜0.001），和P的拟合 方 程 是y=-40.295x+16.857（R2=0.168，P＜0.001），拟合强度为C＞P（图7）。N∶P和N含量的拟合方程是y=4.794x+0.041（R2=880，P＜0.001），和P含量的拟合方程为y=-0.224x+0.333（R2=0003，P=0.456），拟合强度为N＞P（图7）。
图7 宁夏东部风沙区不同固沙灌木配置模式下0-30 cm土壤C、N、P化学计量比与其含量的关系Fig.7 Relationship between stoichiometric ratios of soil C，N，P and their contents in 0-30 cm soil under different sand-fixing shrub configurations in the wind-swept area of eastern Ningxia
3.1 不同固沙灌木配置模式下土壤C、N、P含量特征
草方格耦合固沙灌木有利于土壤养分的积累。一方面，草方格扎设增加了地表粗糙度，使沙面得以固定，从而逐渐形成土壤结皮［20］。并且草方格本身材料多以小麦秸秆［21］、稻草等为主，长期分解后养分归还大地，会提高土壤肥力。另一方面，在草方格中种植柠条、杨柴、花棒、沙拐枣等固沙灌木，不同灌木树种的生活习性及各类配置模式下草本层地上生物量等都对土壤C、N、P含量产生影响，灌木生长过程产生的枝叶枯落物，也会增加土壤养分［22-23］。本研究对5种固沙灌木配置模式下土壤C、N、P 含量进行比较发现，土壤C、N、P 含量在不同固沙灌木配置模式之间存在一定差异。在0～10 cm土层，柠条×花棒混交与柠条×沙拐枣混交固沙灌木配置模式全氮含量最大，显著高于其他3种固沙灌木配置模式，这可能是由于混交林有利于改善土壤物理性质［24］，从而加快枯落物的分解，有利于土壤养分积累，花棒×沙拐枣混交林全氮含量低于含柠条的混交林固沙灌木配置模式。可能是由于花棒和沙拐枣接近地面部分分枝少，通风大，不利用枯落物的积累，同时在研究区内的花棒和沙拐枣死亡率高，导致其全氮含量低于含柠条的混交林，在这3种固沙灌木配置模式中，沙拐枣×花棒混交固沙灌木配置模式的全氮含量高于另外2种纯林，虽然差异不显著。这些结果进一步表明，表层土壤混交林相比纯林更有利于土壤全氮含量的提高，从而提高表层土壤全氮的含量。在20～30 cm土层，柠条固沙灌木配置模式全氮含量最高，显著高于其他4种固沙灌木配置模式；
柠条×花棒混交固沙灌木配置模式全氮含量显著高于另外3种固沙灌木配置模式。这可能是由于研究区内柠条、花棒等固沙灌木细根主要分布在20～30 cm土层区域［25］，豆科灌木根系具有固氮作用，且柠条对土壤全氮的积累高于花棒和杨柴等其他豆科固沙灌木［6］，这些结果表明柠条纯林或柠条与其他豆科灌木混交更有利于较深土壤全氮含量的积累。
土壤有机碳在碳循环中起着至关重要的作用［26］。在0～10 cm土层，柠条×花棒混交、柠条×沙拐枣混交固沙灌木配置模式的有机碳含量高于柠条纯林，这是由于混交林更有利于土壤有机质的积累，从而使表层土壤的有机碳含量较高，这与马云波等［24］的研究结果一致。在20～30 cm土层，柠条×沙拐枣混交固沙灌木配置模式有机碳含量显著低于其他4种固沙灌木配置模式，这可能是由于不同固沙灌木的凋落物和根系分泌物的质量和数量各不相同［27］，土壤有机碳含量受枯落物影响较大，柠条×沙拐枣混交固沙灌木配置模式的地表枯落物含量虽多但分解速率缓慢。
不同固沙灌木配置模式的全磷含量在0～20 cm土层中均不显著，只有在20～30 cm土层中杨柴固沙灌木配置模式下全磷含量显著低于其他4种固沙灌木配置模式，这是由于研究区土壤主要为风沙土，风沙土的发育程度低，对磷的吸附能力也低，所以表层和近表层土壤中全磷含量不受影响，全磷含量出于一种动态平衡状态，固沙灌木对土壤全磷的提高最差，这与之前的研究结果一致［28］。
3.2 不同固沙灌木配置模式下土壤C、N、P生态化学计量比特征
从化学计量学角度阐明沙地恢复过程中，不同固沙灌木混交对土壤养分限制因子供给及养分转化特征的影响。5种固沙灌木样地位于同一个气候区，样地之间的气候差异很小，气候的影响可以忽略不计，土壤化学计量比主要受不同固沙灌木配置模式的影响。土壤碳、氮、磷化学计量特征是确定土壤碳、氮、磷平衡的一个重要指标［9］，对生态学有着重要的指示作用［29］。C∶N通常作为土壤矿化能力指标［30］，反映土壤系统中碳和氮之间的平衡关系，其值越小，表明土壤矿化能力越强，越不利于土壤有机质的积累［31］。本研究中C∶N在0～10 cm土层，柠条×花棒混交和柠条×沙拐枣混交2种固沙灌木配置模式均显著低于其他3种固沙灌木配置模式，主要原因是这2种固沙灌木配置模式下全氮含量显著大于其他3种固沙灌木配置模式，而有机碳含量与其他3种固沙灌木配置模式无显著性差异；
这说明含有柠条的混交林对表层土壤 N含量及其矿化速率影响较大，有利于土壤N释放，促进土壤N转化，不利于表层土壤 C储存［32］。而在20～30 cm土层，柠条×花棒混交、柠条纯林、柠条×沙拐枣混交3种固沙灌木配置模式的C∶N也显著低于其他2种固沙灌木配置模式，这也说明了含有柠条的混交林更有利于土壤的矿化过程，促进土壤养分循环，从而有利于沙地土壤肥力的提高，促进植被恢复。
C∶P作为指示土壤磷素矿化程度的指标，同时也是判断磷的有效性的重要指标，其比值越低表明土壤矿化能力越强［33］。C∶P在10～20 cm土层，柠条固沙灌木配置模式显著高于其他4种固沙灌木配置模式，这主要是由于柠条纯林土壤有机碳含量最大且全磷含量在5种固沙灌木配置模式下无显著差异，这个结果表明柠条纯林固沙灌木配置模式土壤P的有效性较低，这可能是由柠条固沙灌木配置模式土壤有机质的分解速率较低引起的［9］。在20～30 cm土层，柠条×沙拐枣混交固沙灌木配置模式的C∶P显著低于其他4种固沙灌木配置模式，这主要是由于其土壤有机碳含量显著低于其他4种固沙灌木配置模式，而全磷含量显著大于杨柴纯林固沙灌木配置模式且与其他3种固沙灌木配置模式差异不显著，这表明相比其他4种固沙灌木配置模式，柠条×沙拐枣混交固沙灌木配置模式的土壤磷矿化能力更强。
N∶P可以确定土壤养分的限制性阈值［34］，在0～10 cm土层，柠条×花棒混交和柠条×沙拐枣混交2种固沙灌木配置模式N∶P显著高于其他3种固沙灌木配置模式，这主要是由于这2种固沙灌木配置模式全氮含量显著大于其他3种固沙灌木配置模式且全磷含量在5种固沙灌木配置模式之间无显著差异。此结果表明，含有柠条的混交林更有利于土壤氮的矿化过程，对土壤磷含量的下降速率有一定的限制。
3.3 不同固沙灌木配置模式土壤C∶N∶P化学计量比与C、N、P含量的关系
土壤作为生态系统的重要组成部分，其C∶N∶P化学计量比与C、N、P含量之间紧密相关，是衡量土壤元素平衡的重要指标［35］。本研究不同固沙灌木配置模式之间土壤生态化学计量比与C、N、P含量的关系存在差异。C∶N与C、N含量的拟合强度为N大于C，而且，C∶N在5种固沙灌木配置模式下与全氮含量均呈极显著负相关（P＜0.001），而杨柴和花棒×沙拐枣混交2种固沙灌木配置模式的C∶N与有机碳含量呈极显著正相关（P＜0.001），这些结果表明这两种固沙灌木配置模式的土壤C∶N由土壤碳和氮的限制同时决定，而另外3种固沙灌木配置模式的土壤C∶N主要由氮的限制引起。土壤C∶P在5种固沙灌木配置模式下与有机碳间均呈极显著正相关（P＜0.001），而只有柠条纯林、柠条×沙拐枣混交、花棒×沙拐枣混交3种固沙灌木配置模式的C∶P与全磷含量呈显著负相关（P＜0.05），而且土壤C∶P与C、P含量的拟合强度为C大于P。这些结果表明，在不同固沙灌木配置模式下，对土壤C∶P起主导作用的是土壤有机碳含量，而柠条纯林、柠条×沙拐枣混交、花棒×沙拐枣混交3种固沙灌木配置模式的土壤全磷含量也在一定程度上影响土壤C∶P比值。土壤N∶P在5种固沙灌木配置模式下与全氮间均呈极显著正相关（P＜0.001），而与全磷均差异不显著（P＞0.05），而且土壤N∶P与N、P的拟合强度为N大于P。这些结果表明，在土壤元素平衡过程中土壤N作为主要的限制因素，这与陶冶等［36］在古尔班通古特沙漠的研究结果一致。在前人的研究中，C∶N∶P化学计量比与C、N、P含量之间均有一定的线性或非线性关系［30，36-38］。由于研究区域不同，其关系也有不同的具体表现。在本研究中 C∶N与C含量呈线性正相关关系，与N含量呈线性负相关关系，这与张书齐［39］、董雪［8］、邱岭军［40］等研究结果相同；
C∶P与C含量呈线性正相关关系，与P含量不存在线性相关关系。综上所述，本研究C∶N和N∶P主要限制因素均为全氮含量，C∶P主要限制因素为有机碳含量，这一结果与赵力兴等［41］在科尔沁沙地对二龄紫花苜蓿草地研究结果一致。利用土壤化学计量关系来判断土壤养分的限制还需要进一步的野外施肥实验加以印证［42-43］。
宁夏灵武白芨滩国家级自然保护区典型的5种固沙灌木配置模式对土壤碳、氮、磷及其生态化学计量比的影响不同，在表层土壤（0～10 cm），柠条×花棒混交与柠条×沙拐枣混交显著提高土壤全氮含量，2种模式土壤C∶N显著低于其他3种固沙灌木配置模式，土壤N∶P显著高于其他3种固沙灌木配置模式。柠条×花棒混交与柠条×沙拐枣混交固沙灌木配置模式有利于表层土壤的矿化过程，促进土壤N释放；
10～20 cm土层，柠条固沙灌木配置模式C∶P显著高于其他4种，土壤P的有效性较低；
20～30 cm土层，柠条固沙灌木配置模式全氮含量显著高于其他4种，柠条×沙拐枣混交固沙灌木配置模式有机碳含量和C∶P均显著低于其他4种。不同固沙灌木模式下碳、氮、磷及其生态化学计量比具有一定的线性关系，土壤有机碳和全氮是影响化学计量比的主要限制因子，不同固沙灌木土壤元素平衡主要受到土壤氮限制。因此，未来在本区域进行固沙恢复实践中应考虑种植对氮素积累有利的柠条混交林。
