紫色土中氮肥、磷肥、钾肥配施对白芷产量及有效成分的影响

2023-06-08刘仁浪蒋翼杰江美彦陈靳松肖婕妤郑全林吴卫

江苏农业学报 2023年1期

刘仁浪　蒋翼杰　江美彦　陈靳松　肖婕妤　郑全林　吴卫

摘要：在國家防止耕地“非粮化”政策下，为合理利用四川地区紫色土（坡地、荒地等）进行白芷种植，本研究采用3因素5水平二次正交回归旋转组合设计，探究紫色土中氮肥、磷肥、钾肥配施对川芷2号白芷品种产量和有效成分含量的影响。结果表明，紫色土中氮肥、磷肥、钾肥配施对白芷产量的影响大小为钾肥＞磷肥＞氮肥，钾肥和磷肥对产量的影响都达到显著或极显著水平，氮肥对白芷产量没有显著影响；氮肥、磷肥、钾肥配施对欧前胡素含量的影响为钾肥>氮肥>磷肥，对异欧前胡素含量的影响为磷肥>氮肥>钾肥，可见钾肥显著影响欧前胡素含量，磷肥显著影响异欧前胡素含量。通过统计频数法优化出施肥配方，在氮肥施用量为233.73～283.36 kg/hm²、磷肥施用量为1 380.66～1 522.04 kg/hm²、钾肥施用量为215.14～240.03 kg/hm²时川芷2号能实现高产，在氮肥施用量为240.39～276.58 kg/hm²、磷肥施用量为1 097.69～1 197.45 kg/hm²、钾肥施用量为163.12～185.91 kg/hm²时可以均衡地提高川芷2号根中欧前胡素和异欧前胡素含量。最终利用灰色关联度和逼近理想解排序法（DTOPSIS法）综合获得在紫色土中川芷2号的最优施肥方案：氮肥施用量为193.40kg/hm²、磷肥施用量为1 449.26 kg/hm²、钾肥施用量为239.88 kg/hm²。

关键词：紫色土；白芷；配施；产量；欧前胡素；异欧前胡素

中图分类号： S567.9;S143 文献标识码： A 文章编号： 1000-4440（2023）01-0208-10

Effects of combined application of nitrogen fertilizer， phosphorus fertilizer and potassium fertilizer on yield and active components of Angelica dahurica in purple soil

LIU Ren-lang¹， JIANG Yi-jie¹， JIANG Mei-yan¹， CHEN Jin-song¹， XIAO Jie-yu¹， ZHENG Quan-lin²，WU Wei¹

（1.College of Agronomy， Sichuan Agricultural University， Wenjiang 611130， China；2.Suining Tiandi Network Chuanbaizhi Industrial Co.， Ltd.， Suining 629000， China）

Abstract： Under the Chinese policy background of preventing non-grain conversion of cultivated land， a quadratic orthogonal regression rotation combination design using three factors and five levels was adopted in this study to make rational use of purple soil （sloping field， wasteland， etc.） in Sichuan area for Angelica dahurica var. formosana planting and to explore the effects of combined application of nitrogen fertilizer， phosphate fertilizer and potassium fertilizer in purple soil on the yield and effective components content of A. dahurica var. formosana variety of Chuanzhi 2. The results showed that， the effect of combined application of nitrogen fertilizer， phosphate fertilizer and potassium fertilizer on A. dahurica var. formosana yield in purple soil was ranked as follows： potassium fertilizer > phosphate fertilizer > nitrogen fertilizer. Both phosphate fertilizer and potassium fertilizer had significant or extremely significant effect on the yield of A. dahurica var. formosana， while nitrogen fertilizer had no significant effect on the yield of A. dahurica var. formosana. The effect of combined application of nitrogen fertilizer， phosphate fertilizer and potassium fertilizer on the imperatorin content was ranked as potassium fertilizer > nitrogen fertilizer > phosphate fertilizer， while the effect on isoimperatorin content was ranked as phosphate fertilizer > nitrogen fertilizer > potassium fertilizer. It can be seen that potassium fertilizer significantly affected imperatorin content， and phosphate fertilizer significantly affected isoimperatorin content. The fertilization formula was optimized by statistical frequency method， and it was found that Chuanzhi 2 could achieve high yield when the application levels of nitrogen fertilizer， phosphate fertilizer and potassium fertilizer were 233.73-283.36 kg/hm²， 1 380.66-1 522.04 kg/hm²and 215.14-240.03 kg/hm²， respectively. The contents of imperatorin and isoimperatorin in roots of Chuanzhi 2 could equably increase when the application levels of nitrogen fertilizer， phosphate fertilizer and potassium fertilizer were 240.39-276.58 kg/hm²， 1 097.69-1 197.45 kg/hm²and 163.12-185.91 kg/hm²， respectively. Finally， the optimal fertilization scheme for Chuanzhi 2 in purple soil was obtained by comprehensively using the grey relational degree and technique for order preference by similarity to ideal solution （DTOPSIS） method， and the application levels of nitrogenous fertilizer， phosphate fertilizer and potassium fertilizer were 193.40 kg/hm²， 1 449.26 kg/hm²and 239.88 kg/hm²， respectively.

Key words： purple soil;Angelica dahurica;combined application;yield;imperatorin;isoimperatorin

中药白芷[Angelica dahurica （Fisch. ex Hoffm） Benth. et Hook. f.]药食两用，其原植物为伞形科当归属植物白芷（A. dahurica）和杭白芷（A. dahurica var. formosana），具有驱风散寒、清热排脓、消肿止痛等功效^[1^-2]。川白芷（四川）、祁白芷（河北）、禹白芷（河南）、杭白芷（浙江）是目前市场上主要流通的白芷品种^[3]。据《遂宁县志》和《遂宁白芷志》记载，川白芷栽培历史至今已有400～600年，是著名的川产道地药材^[4^-5]，一般种植于涪江两岸的灰潮土上，但目前随着城镇化进程的加快，白芷种植区域从涪江两岸扩张于四川内地的紫色土区域^[6]。《全国土地利用总体规划纲要（2006-2020年）》提出，要确保18亿亩耕地红线，而四川紫色土耕地集中分布于四川盆地丘陵区和海拔800 m以下的低山区，占全省耕地面积的68%^[7^-8]。本课题组调研发现，在国家防止耕地“非粮化”政策下，川白芷种植于紫色土区域（坡地、荒地等）的面积越来越大，且市场认可度和售价都超过用灰潮土种植的白芷。

施肥、产地土壤条件、生长环境、播种时期等因素会影响白芷的生长发育及次生代谢^[9^-14]。前人对白芷施肥的相关研究中大多探究施肥措施对白芷抽薹苔率和产量的影响以及不同肥料种类对川白芷产量的影响，并未见关于紫色土中白芷配方施肥的研究报道^[9^-10，15]。贾全全等^[16]在江西省红壤丘陵区引种4大白芷品种（川白芷、杭白芷、禹白芷、祁白芷），发现杭白芷最适宜栽培于江西省红壤丘陵区。因此应依据白芷产地土壤和品种的特性制定科学合理的施肥措施，本研究采用3因素5水平二次正交回歸旋转组合设计方法，探讨紫色土中氮肥、磷肥、钾肥配施对白芷品种川芷2号产量及品质的影响，以期为川芷2号在紫色土上种植时合理施肥提供科学依据。

1 材料与方法

1.1 试验材料及试验点概况

本试验于2019-2020年在川白芷道地产区四川省遂宁市船山区余健村（30°50′N，105°57′E）进行。试验地为平地，前茬为白芷，土质为紫色土，土壤pH为8，有机质含量13.41 g/kg，全氮含量0.96 g/kg，全磷含量0.73 g/kg，全钾含量17.79 g/kg，碱解氮含量76.03 mg/kg，速效磷含量5.83 mg/kg，速效钾含量144.10 mg/kg。

供试材料为四川农业大学吴卫教授选育的川芷2号白芷品种，原植物为杭白芷，其丰产性、稳定性和适应性良好^[17]。

供试肥料：氮肥为尿素（氮含量≥46.40%），磷肥为过磷酸钙（含P₂O₅≥15%），钾肥为硫酸钾（含K₂O≥50%）。

1.2 试验设计

采用3因素5水平二次正交回归旋转组合设计进行田间试验，设置氮肥（x₁）、磷肥（x₂）、钾肥（x₃）3个因素，各5个水平，共23个小区随机排列，每个小区面积为30 m²（长×宽=5 m×6 m）。二次正交回归旋转组合设计试验因素水平编码见表1，设计与试验方案见表2。

1.3 田间管理

2019年10月采用条播方式将川芷2号播种于遂宁市船山区余健村，当苗高4～7 cm时，进行1～2次间苗，并按10 cm×40 cm株行距定苗。氮肥（尿素）作为苗肥（12月）、冬肥（次年2月）、春肥（次年3月）施用，施用比例为2∶5∶3；磷肥（过磷酸钙）与钾肥（硫酸钾）作为底肥（9月）、春肥（次年3月）施用，施用比例为5∶5。其余田间栽培措施，均与当地农户保持一致。

1.4 测定指标与测定方法

1.4.1 农艺性状及产量指标的测定 2020年8月10日收获时采用5点取样法，每个小区取10株样品，将白芷植株从茎基部分离为地上部分和地下部分。用清水洗干净地下部分（根）后称鲜质量，用直尺测量根长，用游标卡尺测量根直径。将白芷根放置于105 ℃烘箱内杀青15～30 min后，45 ℃烘干至恒质量后称其干质量。各小区挖取3垄白芷称量总的根鲜质量，再折算出产量。

1.4.2 有效成分含量的测定参照《中华人民共和国药典》^[1]（简称《中国药典》）中的方法测定白芷根部欧前胡素及异欧前胡素的含量。

1.5 数据分析与处理

采用Excel 2017、 DPS 7.05对数据进行初步整理与分析，并建立回归模型，使用Excel 2017和GraphPad Prism 8制图。

根折干率=根干质量（g）/根鲜质量（g）×100%

产量（kg/hm²）=挖取白芷总鲜质量（kg）/挖取面积（m²）×10 000

2 结果与分析

2.1 氮肥、磷肥、钾肥配施对川芷2号产量的影响

2.1.1 模型的建立与检验使用DPS7.05统计软件对产量（Y₁）、根鲜质量（Y₂）、根干质量（Y₃）、折干率（Y₄）、根长（Y₅）、根直径（Y₆）、欧前胡素含量（Y₇）、异欧前胡素含量（Y₈）的数据进行回归分析，得到氮肥（x₁）、磷肥（x₂）、钾肥（x₃）施肥量的回归模型。通过F检验分别验证方程各项回归系数、总回归系数和失拟度，得到方差分析结果（表3）。

由表3可知，产量、折干率的方程总回归系数FR（F1R=6.370，F4R=6.244）均达到极显著水平，失拟度F1Lf未达到显著水平（F1Lf=2.117），失拟度F4Lf达到显著水平（F4Lf=3.272），说明产量的二次回归方程模型拟合得较好，其可靠程度较高。根鲜质量、根干质量、根长、根直径、欧前胡素含量、异欧前胡素含量方程的总回归系数FR和失拟度FLf均不显著，说明其二次回归方程拟合效果不好，不可以进行模型决策。在产量的二次回归方程中，一次项（x2、x3）、二次项（x22、x32）、交互项（x2x3）均达到显著或极显著水平。因此，剔除在α=0.05显著水平下不显著的项后，可得到简化后的产量（Y1′）二次回归方程：Y1′=1 390.01+200.08x2+148.56x3-137.77x22-164.42x32+174.27x2x3。

2.1.2 因子主效应分析根据表3中各项回归系数的F值，使用以下贡献率计算公式得出各因素对因变量的贡献率：

在紫色土中氮肥、磷肥、钾肥对川芷2号产量的贡献率分别为：△氮肥=0.346，△磷肥=1.812，△钾肥=1.842，其贡献率大小排序为钾肥＞磷肥＞氮肥，其中钾肥和磷肥对产量的影响最大，氮肥对产量的影响最小。

2.1.3 单因素效应分析由表3可知，单因素项中，氮肥对川芷2号产量没有显著影响，而磷肥和钾肥对其产量的影响达到显著或极显著水平，其中磷肥的影响达到极显著水平。通过降维法，得出磷肥、钾肥对产量的单因素效应方程，由方程作单因素效应图（图1）。

由图1可知，随着磷肥和钾肥施用水平的升高，川芷2号产量先增加后降低，符合肥料报酬定律，说明磷肥和钾肥的施用范围设计合理。磷肥和钾肥在0～1水平时川芷2号的产量达到最高值，说明在紫色土中，当磷肥施用量为1 099.44～1 449.26 kg/hm²、钾肥施用量为179.91～239.88 kg/hm²时有利于川芷2号获得高产。

2.1.4 互作效应分析由表3的回归系数显著性分析可知，磷肥-鉀肥互作对紫色土中川芷2号产量的影响达到了显著水平，而氮肥-磷肥、氮肥-钾肥互作对其产量影响不显著。使用降维法固定氮肥因子使其为0，可得磷肥-钾肥互作与产量的方程，根据方程可作出互作效应图（图2）。

由图2可知，不同磷肥、钾肥配比对紫色土中种植的川芷2号产量影响差异较大，当磷肥、钾肥水平处于低磷高钾和低钾高磷时都对白芷产量有明显的抑制作用。随着磷肥和钾肥施用水平的增加，产量也随着增加。磷肥水平为1.68，钾肥水平为1.34（即磷肥施用量为1 679.14 kg/hm²，钾肥施用量为260.33 kg/hm²）时川芷2号产量最高。说明磷肥-钾肥互作有利于川芷2号产量的增加。

2.1.5 川芷2号施肥方案的优选使用统计频数法进行分析可以得到一个合适的高产范围。由当地土壤、气候条件可知小区产量高于 27 347.68 kg/hm²即为高产。因此，以产量方程作为相关函数目标，采用统计频数法优选出川芷2号产量高于 27 347.68 kg/hm²的氮肥、磷肥、钾肥配施方案。由表4可知，在紫色土中种植的川芷2号产量高于27 347.68 kg/hm²的方案有40个，通过计算95%的置信区间，可优化出稳定的施肥范围，即氮肥施用量为233.73～283.36 kg/hm²、磷肥施用量为1 380.66～1 522.04 kg/hm²、钾肥施用量为215.14～240.03 kg/hm²时，在紫色土中种植的川芷2号产量高于27 347.68 kg/hm²。

2.2 氮肥、磷肥、钾肥配施对川芷2号有效成分的影响

2.2.1 因子主效应分析白芷有效药用成分主要为欧前胡素、异欧前胡素等香豆素类化合物^[18]。利用2.1.2中贡献率公式计算各因素对川芷2号有效成分指标的贡献率。由表5可知，各因素对川芷2号根中的有效成分的影响不同，对欧前胡素含量的影响为钾肥>氮肥>磷肥，对异欧前胡素含量的影响为磷肥>氮肥>钾肥。钾肥对根中欧前胡素含量影响较大，而对根中异欧前胡素含量的影响最小，氮肥处于中间调控地位，磷肥对根中欧前胡素含量影响最小，而对根中异欧前胡素含量的影响较大。

2.2.2 单因素效应分析由表3可知，单因素项中，钾肥显著影响川芷2号根中欧前胡素含量，磷肥显著影响川芷2号根中异欧前胡素含量，通过降维法得到钾肥对欧前胡素含量及磷肥对异欧前胡素含量影响的单因素方程，由方程得到单因素效应图（图3）。

由图3可知，钾肥对欧前胡素含量及磷肥对异欧前胡素含量的影响各不相同。由图3A可知，钾肥对欧前胡素含量的影响整体呈下降的趋势，即随着钾肥施用水平加大，欧前胡素含量呈下降趋势，说明施用过量的钾肥不利于欧前胡素的积累，并对欧前胡素含量的积累有负面影响，钾肥施用量与欧前胡素含量呈负相关。由图3B可知，磷肥对异欧前胡素含量的影响呈先下降后上升趋势，随着磷肥施用水平的加大，异欧前胡素含量先降低再升高。磷肥施用量在-1.68～1.68水平时异欧前胡素含量先降低后再升高，1.68水平时异欧前胡素含量远远大于-1.68水平时异欧前胡素的含量，说明高磷水平比低磷水平更有利于异欧前胡素的积累，磷肥施用量与异欧前胡素含量呈正相关。

2.2.3 提高川芷2号有效成分含量的施肥方案的优选《中国药典》规定，白芷根中欧前胡素含量不得少于0.08%^[1]。根据《中国药典》的规定和当地生产实际，将紫色土上种植的川芷2号根中的欧前胡素含量≥0.28%，异欧前胡素含量≥0.15%视为优质白芷标准。通过统计频数法得到欧前胡素含量大于0.28%的方案有55个（表6）；异欧前胡素含量大于0.15%的方案有75个（表7）。

由表6、表7可知，通過计算95%的置信区间，得出当氮肥施用量为233.73～277.06 kg/hm²，磷肥施用量为1 001.35～1 197.45 kg/hm²，钾肥施用量为144.83～185.91 kg/hm²时川芷2号根中的欧前胡素含量可达到0.28%及以上；当氮肥施用量为240.39～276.58 kg/hm²，磷肥施用量为1 097.69～1 317.05 kg/hm²，钾肥施用量为163.12～196.70 kg/hm²时川芷2号根中的异欧前胡素含量可达到0.15%及以上。再进行交集优化后可得出当氮肥施用量为240.39～276.58 kg/hm²，磷肥施用量为1 097.69～1 197.45 kg/hm²，钾肥施用量为163.12～185.91 kg/hm²时可以均衡地提高紫色土中种植的川芷2号根中欧前胡素和异欧前胡素含量。

2.3 基于灰色关联度和DTOPSIS法综合评价川芷2号的产量及有效成分含量

为了精准探究在紫色土中种植川芷2号的施肥配方，本研究综合了产量、欧前胡素含量、异欧前胡素含量，通过灰色关联度法得到其关联度和权重（表8）。通过关联度和权重进行DTOPSIS法处理，根据DTOPSIS法中的C_i（相对接近度）进行排序，并计算出待评材料的γ_i（关联度）（表9）。

由表9可知，灰色关联度法中最大关联度差异为45.15%，而DTOPSIS法中最大C_i差异为71.03%，说明DTOPSIS法相比灰色关联度法放大了处理间的差异性，更有利于区分待评材料间的优劣。处理5、处理12和处理21等处理采用灰色关联度和DTOPSIS法得到的综合评价排序结果一致，其中综合表现最优的是处理5，即当氮肥施用量为193.40 kg/hm²，磷肥施用量为1 449.26 kg/hm²，钾肥施用量为239.88 kg/hm²（氮肥∶磷肥∶钾肥=0.75∶1.81∶1.00）时在紫色土中川芷2号的产量和有效成分含量综合最优。

3 讨论

遂宁市位于四川盆地中部丘陵低山地区，其土壤以紫红色砂土和泥岩为主^[19]。由于川白芷种植区域由涪江冲积坝上（灰潮土）逐渐向丘陵低山地区扩张，导致白芷种植在紫色土上的面积越来越大。本课题组在川芷2号的多点试验中发现紫色土中种植川芷2号的产量与灰潮土基本一致，但其欧前胡素和异欧前胡素含量略高于灰潮土，可能是由于紫色土母岩疏松，易崩解，所以矿质养分含量高，肥力好，有机质含量较低，磷、钾等含量丰富^[20]。本研究中在紫色土中种植的川芷2号氮肥、磷肥、钾肥的配比为0.75∶1.81∶1.00，其产量和欧前胡素含量最高。陈郡雯等^[21]研究发现灰潮土中白芷氮肥、磷肥、钾肥配比为1.3∶1.9∶1.0时白芷产量和欧前胡素含量最高。这可能是由于紫色土与灰潮土的土壤母质、理化性质、养分状况及分布形态存在差异，从而对白芷产量和品质产生不同影响^[22^-24]；也可能是不同白芷品种需肥量有所差异。紫色土与灰潮土对川芷2号产量和品质的影响因素还需进一步探究。

土壤耕作层速效氮含量低，速效磷和速效钾含量较为丰富是白芷高产的基础^[11]。从本研究单因素分析中可知，氮肥对紫色土中种植的川芷2号产量影响不显著，磷肥对产量有极显著影响，钾肥对产量有显著影响，同时磷肥-钾肥对紫色土中种植的川芷2号产量还有交互作用。翟娟园等^[13]研究发现，土壤有效磷及速效钾的含量与白芷根质量呈现正相关，这与本研究的结果一致。

白芷根中总香豆素主要为欧前胡素和异欧前胡素，欧前胡素和异欧前胡素含量能反映川白芷的品质^[25^-26]。贾蕾等^[27]研究发现，磷肥和钾肥有利于提高白芷根中香豆素的含量。陈郡雯等^[28]研究发现，川白芷在生长过程中异欧前胡素含量与根部氮、磷、钾的积累量呈显著正相关。在本研究中川芷2号根中的欧前胡素和异欧前胡素含量均大于0.1%，符合《中国药典》的规定。本研究中氮肥、磷肥、钾肥配施对白芷欧前胡素含量的影响为钾肥>氮肥>磷肥，对异欧前胡素含量的影响为磷肥>氮肥>钾肥，磷肥显著影响异欧前胡素含量，钾肥显著影响欧前胡素含量，磷肥施用量与异欧前胡素含量呈正相关，而钾肥施用量与欧前胡素含量呈负相关。陈郡雯等^[21]在氮肥、磷肥、钾肥配施对川白芷欧前胡素和异欧前胡素含量影响的研究中发现，钾肥施用量与欧前胡素含量呈负相关，磷肥施用量与异欧前胡素含量呈正相关，本研究结果与其一致。因此在四川白芷主产区紫色土上种植白芷，只要选择合理的氮肥、磷肥、钾肥施用量及配比，就可以保障其产量和品质，白芷产量最高可达27 347.68 kg/hm²，欧前胡素和异欧前胡素含量分别可达0.28%、0.15%。这为新产区白芷生产提供了科学施肥指导，为合理利用坡地、荒地等开展白芷生产奠定了基础。

4 结论

本试验通过3因素5水平二次正交回归旋转组合设计，建立了以氮肥、磷肥、钾肥施用量为自变量，在紫色土中种植的白芷产量和有效成分含量为因变量的函数模型，通过模型分析得出，磷肥和钾肥施用量对川芷2号产量有显著影响，磷肥-钾肥对产量还有交互作用。利用频数分析法优化出紫色土中种植的川芷2号高产的施肥配方：氮肥施用量为233.73～283.36 kg/hm²、磷肥施用量为1 380.66～1 522.04 kg/hm²、钾肥施用量为215.142～240.03 kg/hm²。钾肥施用量对欧前胡素含量有显著作用，磷肥施用量对异欧前胡素含量有显著作用，使用频数分析法进行交集优化后可得出，当氮肥施用量为240.39～276.58 kg/hm²，磷肥施用量为1 097.69～1 197.45 kg/hm²，钾肥施用量为163.12～185.91 kg/hm²时可以均衡地提高紫色土中种植的川芷2号根中欧前胡素和异欧前胡素含量。最后通过灰色关联度和DTOPSIS法综合评价同时兼顾紫色土中川芷2号的产量和有效成分含量的最优施肥方案：氮肥施用量为193.40kg/hm²，磷肥施用量为1 449.26 kg/hm²，钾肥施用量为239.88 kg/hm²。

参考文献：

[1] 中华人民共和国国家药典委员会. 中华人民共和国药典 [M]. 北京：中国医药科技出版社， 2020.

[2] WEI W L， ZENG R， GU C M， et al. Angelica sinensis in China-A review of botanical profile， ethnopharmacology， phytochemistry and chemical analysis[J]. Journal of Ethnopharmacology， 2016， 190： 116-141.

[3] 张庆芝，王开疆，刀莉芳. 中药白芷的品种论述 [J]. 云南中医学院学报， 2000， 23（2）： 22-24.

[4] 王梦月，贾敏如. 白芷本草考证 [J]. 中药材， 2004， 27（5）： 382-385.

[5] 遂宁市地方志编篡委员会. 遂宁县志 [M]. 成都：巴蜀书社， 1992.

[6] 郑利，邓聪，冯亮，等. 遂宁川白芷产业发展现状与分析 [J]. 农村经济与科技， 2020， 31（20）： 181-182.

[7] 施业家，吴贤静. 生态红线概念规范化探讨 [J]. 中南民族大学学报（人文社会科学版）， 2016， 36（3）： 149-153.

[8] 中国科学院成都分院土壤研究室. 中国紫色土（上篇）[M]. 北京：科学出版社， 1991.

[9] 童文，孙佩，杨晓，等. 施用腐植酸肥对白芷产量和质量的影响 [J]. 西南农业学报， 2011， 24（3）： 1236-1238.

[10]张亚琴，雷飞益，陈雨，等. 锌硼钼配施对川白芷药材农艺性状与产量的影响 [J]. 植物营养与肥料学报， 2018， 24（3）： 769-778.

[11]陈兴福，丁德蓉，刘岁荣，等. 白芷生长土壤的研究 [J]. 中国中药杂志， 1994， 19（10）： 591-593.

[12]陈兴福，丁德蓉，刘岁荣，等. 白芷生态环境和土壤理化特性的研究 [J]. 中草药， 1996， 27（8）： 489-591.

[13]翟娟园，吴卫，廖凯，等. 土壤环境对川白芷产量和品质的影响研究 [J]. 中草药， 2010， 41（6）： 984-988.

[14]陈兴福，卢进，丁德蓉. 播种期对白芷早期抽苔影响的研究 [J]. 中国中药杂志， 1999， 24（4）： 211-212.