研究区位于四川省泸州市古蔺县国有林场（E105°34′—105°52′，N27°56′—28°11′）。该区域地处大娄山褶皱带西段北侧，属四川盆地南缘山地，为亚热带湿润季风气候，≥ 0 ℃年积温5200 ℃，多年均温14.5 ℃，1月平均气温5.3 ℃，7月平均气温26.6 ℃，无霜期约200 d，年平均降雨量1100 mm。该区域温凉高湿，秋冬多雾，空气湿度大，9月—11月（出笋期）日均空气湿度 ≥ 95%。研究区内竹类资源以刺竹子和毛竹为主，灌草多为铃木冬青（Ilex suzukii）、黄毛楤木（Aralia chinensis）、马桑（Coriaria nepalensis）、悬钩子属（Rubus）、蕨类（Pteridophyta）和蒿类（Artemisia）等。土壤类型为黄棕壤，土层厚度大于40 cm，呈酸性。试验林分为刺竹子纯林，郁闭度在0.9以上，林下无灌草。平均立竹密度为10万株·hm⁻²。试验林分基本情况如表1。

2017年6月在研究区海拔1780 ± 50 m处选择具有代表性的刺竹子纯林为对象，在林内建立15个5 m×5 m的观测样方并编号，样方间设置大于3 m的缓冲带。试验设置了5种密度水平：即10 万株·hm⁻² （CK）、4万株·hm⁻²（D1）、5万株·hm⁻²（D2）、6万株·hm⁻²（D3）和7万株·hm⁻²（D4），每个处理3个重复，随机分布于15个样方内。

本研究中密度调控采用间伐手段进行，间伐时间为2017年6月、2018年6月和2019年6月，每次间伐时砍伐密集生长、竞争激烈区域的成竹，保留长势良好的成竹，使保留下的刺竹子数量达到设计密度并均匀分布在样方内，将间伐剩余物清除至样地外。

在2019年9月15日—10月30日期间（根据前期调查，刺竹出笋期在9月中旬至10月底），每间隔2 d进入试验样地采集每个样方中高于土面30 cm~60 cm的竹笋（结合当地采笋习惯，林地中 > 60 cm的竹笋为留竹），并使用电子计价秤（型号：HY-799 ACS，精度：0.01 kg）就地称取带皮鲜笋质量，采集笋产量数据。

2019年6月下旬进行叶片样品的采集。在各样方中随机选择3株2~3年生，无病虫害的刺竹子为采样株，采集植株上、中、下部位，东、南、西、北4个方向的叶片各50 g左右，混合均匀后置于冰盒中尽快带回实验室。将带回的叶片样品经去离子水冲洗后置于105 ℃的烘箱中杀青0.5 h，将杀青后的叶片在80 ℃下烘干至恒重，然后粉碎，研磨过0.15 mm筛用于氮、磷含量的测定。

叶片氮含量采用硫酸-高氯酸消煮法消煮后在全自动间断化学分析仪（法国 Alliance Smart Chem 200）上分析测定；叶片磷含量采用硫酸-高氯酸-钼锑抗比色法测定。所有样品均进行水分系数的转换，测定结果以单位质量的元素含量（g·kg⁻¹）表示^[18]。

运用Excel 2019进行数据整理和绘图，应用IBM SPSS 25进行数据统计分析。采用单因素方差分析法（one-way ANOVA）分析密度调控对刺竹子叶片养分特征和笋产量进行方差分析，采用Pearson（双侧）相关性分析对叶片养分特征与笋产量进行相关性分析。

由图1可知，各密度处理下刺竹子叶片氮、磷含量和氮磷比分别介于21.49~26.85 g·kg⁻¹、0.81~1.23 g·kg⁻¹和21.42~28.90之间。与CK相比，D1和D2处理显著降低了叶片氮含量（P < 0.05），对叶片磷含量和氮磷比无显著影响（P > 0.05）；D3和D4处理显著增加了叶片磷含量（P < 0.05），显著降低了叶片氮磷比（P < 0.05），对叶片氮含量无显著影响（P > 0.05）。

图  1  密度调控对刺竹子叶片氮、磷养分特征的影响

Figure 1.  Effects of density regulation on leaf nitrogen and phosphorus nutrient characteristics of Chimonobambusa pachystachys

由图2可知，各处理下刺竹子笋产量介于2761.42~3690.05 kg·hm⁻²·a⁻¹之间。与CK相比，D1和D2处理的笋产量分别降低了273.15·kg·hm⁻²·a⁻¹和81.40·kg·hm⁻²·a⁻¹，减幅分别为9.00%和2.68%；D3和D4处理的笋产量分别提高了655.48·kg·hm⁻²·a⁻¹和225.54·kg·hm⁻²·a⁻¹，增幅分别为21.60%和7.43%。方差分析表明，D3处理显著提高了刺竹子笋产量（P < 0.05），其余处理对刺竹子笋产量无显著影响（P > 0.05）。

图  2  密度调控对笋产量的影响

Figure 2.  Effect of density regulation on Chimonobambusa pachystachys bamboo shoot yields

由表2可知，刺竹子笋产量与叶片氮、磷含量呈显著正相关（P < 0.01），与叶片氮磷比呈显著负相关（P < 0.01）。

林分密度调节着植被生境的光照、水分和土壤养分等环境因子，是林木叶片养分和地上生产力的重要影响因素^{[17, 19]}。本研究中，刺竹子叶片氮、磷含量随立竹密度的增加整体表现出先上升后下降的趋势（见图1A和图1B）。其原因是密度调控直接改变了林地内的立竹数量和生物量分配格局^[6]，改变了林地的凋落物分解量和养分归还过程，进而对叶片氮、磷含量产生了影响。有研究表明，森林植物吸收的养分中约90%来自凋落物分解后归还给土壤的养分^[19]，低密度下竹林凋落物较少，养分归还量低，植物可利用的养分来源减少，因而植株叶片氮、磷养分含量降低；而高密度下立竹数量增加会加剧养分的竞争和消耗，竹林总体上干物质积累较多，在一定程度上引起养分稀释^[20]，并且立竹密度增大也会引起土壤速效养分需求量的增加，导致土壤养分的可获得性和生物有效性降低，从而引起植株养分吸收限制和获取量的减少^[21]，因而叶片氮、磷含量呈下降趋势。值得注意的是，本研究中D3和D4处理下叶片氮含量与CK差异不显著（见图1A），其原因是相较于低密度处理（D1和D2），高密度（D3、D4和CK）处理的竹林中植株可能同化更多的氮素来增强其羧化能力以竞争更多光照和生存空间^[22]，光合羧化能力的增强可能是导致高密度下竹株叶片氮含量未呈现显著差异的原因。

植物叶片的氮磷比值不仅反映着植物群落的养分利用情况，也能对林地生产力的限制性元素提供指示^[9]。本研究中各处理下刺竹子叶片氮磷比均大于16（见图1C），表明研究区刺竹子生长易受到磷素的限制^[9]，这与本课题组在同区域的前期研究^[23]和大多数亚热带区域植物养分限制的研究结果^[24]一致。因此，可对研究区刺竹子林分适当施用磷肥。本研究还发现，密度调控显著改变了刺竹子叶片的氮磷比。郭子武等对不同密度下四季竹养分重吸收特征的研究结果也表明，为适应不同立竹密度下养分和空间资源的差异，植物会调节其营养元素的重吸收量^[5,21]，进而改变植株叶片的氮磷比。生长速率假说指出，生物体较低的氮磷比有利于核糖体快速合成蛋白质^[25]，促进植被生长。本研究中，D3和D4立竹密度下刺竹子叶片氮磷比显著低于其他处理，表明将立竹密度控制在D3和D4水平时氮、磷养分较为均衡，有利于研究区刺竹子的生长。

林分结构和植株养分含量是刺竹子笋产量的关键影响因子^[2]。本研究中，刺竹子的竹笋产量随林分密度的增大呈先增加后减少的规律，当立竹密度为D3时竹笋产量达到最大值（见图2）。一方面是密度调控直接改变了立竹的空间结构和母竹数量，进而影响了刺竹子林分的笋产量。高密度下虽然母竹数量较多，但竹株对林地空间与养分的竞争加剧^[4]会导致竹笋减产；而在低密度下虽然林地空间充足，立竹竞争小，但母竹数量减少也会导致出笋数量的下降，引起笋产量降低^[2]。另一方面，密度调控通过调节竹株的氮、磷养分状况改变了竹笋产量（见图1）。一般而言，植物体内养分含量越高，养分越均衡，其地上生产力也就越高^[26]。本研究中，相关性分析结果也表明，刺竹子笋产量与叶片氮、磷含量呈正相关，与叶片氮磷比呈负相关关系（见表2），D3处理下刺竹子叶片氮、磷含量高，氮磷比低，因此立竹密度为D3时笋产量较高。因此，研究区刺竹子立竹密度保留在6万株 hm⁻²时有利于笋产量的提高。

研究区刺竹子生长易受到磷素限制，可对竹林适当施磷肥。密度调控显著改变了刺竹子叶片氮、磷养分特征，当立竹密度为6~7万株·hm⁻²时养分限制情况得到缓解。刺竹子叶片氮、磷养分含量与笋产量呈正相关关系，氮磷比与笋产量呈负相关关系。当立竹密度调控为6万株·hm⁻²时，刺竹子的笋产量得到了显著提高。竹笋产量除了受到立竹密度、植株养分的影响外，也与气候、土壤肥力等环境因子密切联系，在今后的研究中可进一步探讨这些因子与笋产量的关系，为刺竹子笋用林的可持续经营提供依据。