Soil Nutrient Status and Comprehensive Evaluation of Soil Fertility Quality of Tea Garden in Anxi
-
摘要:目的 茶叶产量和品质与茶园土壤养分状况及肥力质量密切相关。然而，目前针对茶园肥力状况的调查评价依旧不足。方法 本研究对福建省安溪县13个主要产茶乡镇共38个铁观音生产茶园进行土壤样品采集与理化指标检测，并通过对土壤pH值以及关键营养元素分析，建立土壤肥力综合指数（IFI）。结果 结果表明，茶园土壤pH值平均值为4.0，有近55%样点茶园土壤酸化严重；茶园土壤的有机质、全氮、全磷、全钾、碱解氮、有效磷、速效钾等指标相对优质标准较高；综合各评价指标，得到IFI值范围在0.52 ~ 1.00，平均值为0.84，有84%样点的茶园处于 Ⅰ 级优质标准。结论 安溪县铁观音茶园土壤大半已严重酸化，但安溪茶园土壤整体养分状况优良，茶园土壤肥力综合指数IFI较高，表明安溪铁观音茶园土壤总体较健康。Abstract:Objective The nutrient status and fertility quality of tea garden soil directly affect the quality of tea. However, the investigation and assessment of soil health in tea garden is still insufficient.Method The physical and chemical indices of soil samples from 38 Tieguanyin tea gardens in 13 main tea-producing towns in Anxi were detected, and a comprehensive evaluation with the integrated fertility index (IFI) was constructed by measuring soil pH value and key nutrient elements.Result The results showed that the average pH value of tea garden soil was 4.0, and nearly 55% of tea garden soil was strongly acidic. The indies of soil organic matter, total nitrogen, total phosphorus, total potassium, alkali hydrolyzable nitrogen, available phosphorus and available potassium in tea garden soil were relatively high. According to the comprehensive evaluation indices, the IFI value ranged from 0.52 to 1.00, with an average value of 0.84. The percentages of 84 in the tea gardens were in class I high-quality standard.Conclusion Most of the Tieguanyin tea garden soil in Anxi has been seriously acidified, but the overall nutrient status of Anxi tea garden soil is still in good condition, and the comprehensive index of soil fertility (IFI) of the tea garden is high, indicating that the soil is generally healthy in Anxi Tieguanyin tea garden.
-
【研究意义】中国是茶叶的发源地，茶叶中因富含对人体健康有益的物质而备受欢迎[1]。铁观音是中国十大名茶之一，安溪县作为铁观音茶叶主要生产区[2]，安溪铁观音中国特色农产品优势区也被认定为第三批中国特色农产品优势区。【前人研究进展】随着生活水平日益提高，人们对茶叶品质要求也越来越高。茶园土壤肥力优良与否是影响茶叶品质好坏的关键因素[3]，施肥方式影响茶园土壤地力。近几十年来，安溪铁观音茶园经历了长期单一种植和不合理施肥，土壤酸化、重金属与潜在风险元素污染、养分不平衡等问题接踵而至，这不仅影响了茶园土壤健康，更影响了茶叶品质，导致经济效益下滑，严重可能危害人体健康[4] 。然而，目前对安溪县茶园土壤调查局限在某一养分是否丰缺，对综合肥力因子的茶园土壤综合评价较少[5]，导致对安溪县茶园土壤认识不足。
【本研究切入点】健康土壤生态产品助力乡村振兴，目前土壤健康面临污染加剧的风险，化肥、农药过度施用造成土壤污染、酸化、退化严重[6]，低耕地质量与全球气候变化导致土壤变质、有机质含量下降[6-9]。有机肥因为能改善土壤的板结性状，提高土壤有机质含量[5]，长期施用对土壤和茶叶质量都有好处而受到大力推广[7-9]。安溪县域内各茶园积极响应国家政策，施用有机肥替代无机化肥，但茶园施用有机肥料后其土壤养分状况是否改善尚未明确，土壤肥力质量也无法据此得知。【拟解决的问题】因此，我们对安溪县铁观音主要产区的茶园土壤进行了采集与分析，以期全面了解当前茶园土壤养分状况及肥力水平，旨在为铁观音茶园土壤养分管理和肥力提升提供更加科学的指导，保障茶叶的营养健康，助力乡村振兴，促进共同富裕。
1. 材料与方法
1.1 土壤样品采集
安溪县（24°50′ ~ 25°26′ N，117°36′ ~ 118°17′ E）是福建省泉州市辖县，共有15个镇、9个乡，由于地形地貌的差异，东部外和西部内安溪有着明显不同的气候特点，东部外安溪属南亚热带，农作物一年三熟；西部内安溪四季较分明，农作物一般一年两熟。
土壤样品采集自安溪县铁观音茶园（如图1），分布在感德镇、虎邱镇、龙涓乡等共13个乡镇，共38个茶园土壤采样点。所有样点均采用棋盘式的采样方法，采集耕层（0 ~ 20 cm）的土壤样品，每个采样点采集16份土样，用四分法最终混合成4份平行样。土壤样品采集时间为2020年9月下旬至10月初。采集好的土壤经自然风干、研磨后，分别过10目、100目、200目筛，用于测定不同理化指标。
1.2 土壤肥力相关理化指标的测定
本实验测定了茶园土壤的酸碱度（pH）、有机质（OM）、全氮（TN）、全磷（TP）、全钾（TK）、碱解氮（AN）、有效磷（AP）和速效钾（AK）。
土壤pH值测定采用电位测定法，称取4 g过筛后（10目）的土壤与10 mL无菌水混合成土壤悬液，振荡2 h后静置，使用pH计（IS126C Icon，上海仪迈仪器科技有限公司）读取pH值。
土壤全磷、全钾含量测定使用X射线荧光光谱仪（Axios-MAX，荷兰帕纳科公司），取4 g过200目筛孔的土壤制成压片，进行元素的分析；土壤全氮含量测定使用CNS元素分析仪（Vario MAX，Elementar Analysensysteme GmbH），取200 mg过200目筛孔的土壤在机器上测定。
土壤有机质含量测定参照农业部发布的《中华人民共和国农业行业标准》土壤检测第6部分（NY/T 1121.6—2006），采用重铬酸钾容量法进行测定；土壤有效磷含量测定参照土壤检测第7部分（NY/T 1121.7—2014），采用紫外/可见分光光度计法进行测定。土壤速效钾含量测定参照《中华人民共和国国家标准》森林土壤速效钾的测定（GB 7856—87），采用1 mol L–1乙酸铵浸提－火焰光度法进行测定；土壤碱解氮含量测定参照《中华人民共和国林业行业标准》森林土壤水解性氮的测定（LY/T 1229—1999），采用碱解-扩散法进行测定。以上土壤样品均过100目筛孔。
1.3 茶园土壤养分评价标准
土壤养分状况评价标准参照《中华人民共和国农业行业标准》茶叶产地环境技术条件（Y/T 853—2004）中茶园土壤肥力分级标准和国家绿色食品产地环境质量标准（NY/T 391—2013），将茶园土壤养分分为 Ⅰ、Ⅱ、Ⅲ 3个等级，茶园土壤养分含量分级标准如表1：
表 1 茶园土壤养分分级标准Table 1. Classification standard of soil nutrients in tea garden采样点
Sampling pointⅠ Ⅱ Ⅲ pH 4.5 ~ 5.5 4.0 ~ 4.5或5.5 ~ 6.5 < 4.0或 > 6.5 有机质 （g kg–1） > 15 10 ~ 15 < 10 全氮（g kg–1） > 1.0 0.8 ~ 1.0 < 0.8 全磷（g kg–1） > 0.6 0.4 ~ 0.6 < 0.4 全钾（g kg–1） > 10 5 ~ 10 < 5 碱解氮（mg kg–1） > 100 50 ~ 100 < 50 有效磷（mg kg–1） > 10 5 ~ 10 < 5 速效钾（mg kg–1） > 120 80 ~ 120 < 80 1.4 土壤肥力综合指数计算
1.4.1 土壤单项肥力指标隶属度的确定
本研究选取pH、有机质、全氮、全磷、全钾、碱解氮、有效磷、速效钾共8个指标确定安溪县土壤肥力综合指标值。常用的隶属度函数有S型和抛物线型两种[10-13]，其中，pH属于抛物型隶属函数，其余都属于S型隶属函数。函数表达式如下：
抛物线型隶属函数：
$$ {A}_{i}=f\left(x\right)=\left\{\begin{array}{*{20}{c}} 1.0-\dfrac{0.9\left(x-{x}_{3}\right)}{{x}_{4}-{x}_{3}}&{x}_{3} < x < {x}_{4}\\ 1.0&{x}_{2}\le x\le {x}_{3}\\ 0.1 + \dfrac{0.9\left(x-{x}_{1}\right)}{{x}_{2}-{x}_{1}}&{x}_{1} < x < {x}_{2}\\ 0.1 &x\ge {x}_{4}或x\le {x}_{1} \end{array}\right. $$ (1) S型隶属函数：
$$ {A}_{i}=f\left(x\right)=\left\{\begin{array}{*{20}{c}} 1.0 &x\ge {x}_{2}\\ \dfrac{0.9\left(x-{x}_{1}\right)}{{x}_{2}-{x}_{1}} + 0.1 & {x}_{1} < x < {x}_{2}\\ 0.1&x\le {x}_{1} \end{array}\right. $$ (2) 式中：Ai 为评价因子i的隶属值，x为实测值，隶属值大小为0.1 ~ 1.0。各肥力指标的隶属度函数曲线转折点如表2：
表 2 各肥力指标的隶属度函数曲线转折点的取值Table 2. Value of turning point of membership function curve of each fertility index转折点取值
Value of turning pointpH 有机质
OM全氮
TN全磷
TP全钾
TK碱解氮
AN有效磷
AP速效钾
AKx1 4.0 10 0.8 0.4 5 50 5 80 x2 4.5 15 1.0 0.6 10 100 10 120 x3 5.5 − − − − − − − x4 6.5 − − − − − − − 1.4.2 土壤单项肥力指标权重系数的确定
采用主成分分析法，选取累计百分率大于80% 的主成分来确定土壤肥力因子的权重值[10-13]，计算公式如下：
$$ {B}_{i}=\frac{{r}_{i\mathrm{均}}}{\displaystyle \sum {r}_{\mathrm{均}}} $$ (3) 式中：ri均 为评价因子i与其他各因子间相关系数的平均值，
${\displaystyle \sum {r}_{\mathrm{均}}} $ 表示所有因子相关系数平均值和。1.4.3 土壤肥力综合指数IFI的计算
土壤肥力综合指数IFI的计算公式如下：
$$ \mathrm{I}\mathrm{F}\mathrm{I}=\sum _{i=1}^{n}{B}_{i}\times {A}_{i} $$ (4) 式中：n 为评价指标数量，Bi 和Ai 分别为第i个评价指标的权重系数和隶属值。将土壤肥力综合指数IFI进行等级划分[10-13]，如表3：
表 3 土壤肥力综合指数划分表Table 3. Division of comprehensive index of soil fertility等级
LevelⅠ Ⅱ Ⅲ IV 土壤肥力综合指数IFI IFI ≥ 0.75 0.50 ≤ IFI ≤ 0.75 0.25≤IFI < 0.50 IFI < 0.25 1.5 数据分析
采用Excel 2019进行土壤养分含量数据处理并求算平均数、变异系数及肥力综合指数，使用MATLAB 2014计算土壤单项肥力指标隶属度，使用IBM SPSS 22.0进行主成分分析确定权重系数，使用Arcmap10.7绘制养分含量点状图与柱状图等。
2. 结果与分析
2.1 茶园土壤养分状况
2.1.1 茶园土壤pH值
安溪各乡镇铁观音茶园土壤pH值在3.3 ~ 5.3之间，平均值为4.0，变异系数为12.09 %，属于低变异，空间差异性较小，说明茶园土壤普遍酸化。如图2-1所示，共7个样点的茶园土壤酸碱度达到 Ⅰ 级标准，占总茶园样点数的19%；10个样点的茶园土壤酸碱度达到 Ⅱ 级标准，占总茶园样点数的26%；21个样点的茶园土壤酸碱度处于 Ⅲ 级标准，占总茶园样点数的55%，茶园土壤酸化较为严重。
![]()
图 2 茶园土壤养分状况含量图茶园土壤养分状况含量图分为点状图与柱状图。点状图用大、中、小三个点表示，分别代表 Ⅰ 级、Ⅱ 级、Ⅲ 级标准，柱状图用三条柱体表示，其中中间柱体代表测量平均值，用Av表示，两边柱体为固定值，代表了分级标准中的两个临界值。如果中间柱体高于右边柱体，代表该处土壤处于 Ⅰ 级标准，低于左侧则代表该处土壤处于 Ⅲ 级标准，处于两侧之间则代表该处土壤处于 Ⅱ 级标准。Figure 2. Content map of soil nutrients in tea garden安溪各乡镇铁观音茶园土壤全钾含量在7 ~ 37 g kg–1之间，平均值为15 g·kg–1 ，变异系数为45.52%，属于中等变异。如图2-7，共29个样点的茶园土壤达到 Ⅰ 级标准，占总茶园样点数的76%；9个样点的茶园土壤达到 Ⅱ 级标准，占总茶园样点数的24%。
安溪各乡镇铁观音茶园土壤有效钾含量在64 ~ 336 mg kg–1之间，平均值为125 mg kg–1，变异系数为45.71%，属于中等变异。如图2-8，有16个样点的茶园土壤达到 Ⅰ 级标准，占总茶园样点数的42%；有14个样点的茶园土壤达到 Ⅱ 级标准，占总茶园样点数的37%；8个样点的茶园土壤达到 Ⅲ 级标准，占总茶园样点数的21%。
2.1.2 茶园土壤有机质含量
安溪各乡镇铁观音茶园土壤有机质含量在20 ~ 59 g kg–1之间，平均值为35 g kg–1，变异系数为26.36%，属于中等变异。如图2-2所示，所有茶园有机质含量都达到了 Ⅰ 级优质茶园标准，茶园土壤的有机质含量丰富。
2.1.3 茶园土壤氮素含量
安溪各乡镇铁观音茶园土壤全氮含量在0.6 ~ 2.1 g kg–1之间，平均值为1.5 g·kg–1，变异系数为26.07%，属于中等变异。如图2-3所示，共31个样点的茶园土壤达到 Ⅰ 级标准，占总茶园样点数的81%；6个样点的茶园土壤达到 Ⅱ 级标准，占总茶园样点数的16%；1个样点的茶园土壤达到 Ⅲ 级标准，占总茶园样点数的3%。
安溪各乡镇铁观音茶园土壤碱解氮含量在44 ~ 288 mg·kg–1之间，平均值为162 mg·kg–1，变异系数为38.49%，属于中等变异。如图2-4所示，共有30个样点的茶园土壤达到 Ⅰ 级标准，占总茶园样点数的79%；7个样点的茶园土壤达到 Ⅱ 级标准，占总茶园样点数的18%；1个样点的茶园土壤达到 Ⅲ 级标准，占总茶园数的3%。
2.1.4 茶园土壤磷素含量
安溪各乡镇铁观音茶园土壤全磷含量在0.3 ~ 2.2 g kg–1之间，平均值为1.0 g·kg–1，变异系数为43.71%，属于中等变异。如图2-5所示，共31个样点的茶园土壤达到 Ⅰ 级标准，占总茶园样点数的81%；6个样点的茶园土壤达到 Ⅱ 级标准，占总茶园样点数的16%；1个样点的茶园土壤为 Ⅲ 级标准，占总茶园样点数的3%。
安溪各乡镇铁观音茶园土壤有效磷含量在7 ~ 254 mg kg–1之间，平均值为84 mg kg–1，远高于10 mg kg–1，变异系数为78.84%，属于强变异，说明有效磷含量分布空间差异大。如图2-6所示，共34个样点的茶园土壤达到 Ⅰ 级标准，占总茶园数的89%；4个样点的茶园土壤达到 Ⅱ 级标准，占总茶园数的11%。
2.1.5 茶园土壤钾素含量
2.2 茶园土壤肥力综合评价
2.2.1 单项土壤肥力指标的隶属度与权重系数
根据1.4.1公式（1）（2），计算各指标的隶属度，如表4：
表 4 安溪茶园土壤肥力指标的隶属度Table 4. Membership degree of soil fertility index in Anxi tea garden乡镇
Township村落/地点
PointpH 有机质
OM全氮
TN全磷
TP全钾
TK碱解氮
AN速效磷
AP速效钾
AK桃舟乡 吾培村 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.775 产业园 0.1 1.0 1.0 1.0 0.46 1.0 1.0 1.0 莲山村 0.1 1.0 1.0 1.0 0.46 1.0 1.0 0.91 福田乡 双按村 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.1 尾洋村 0.1 1.0 1.0 1.0 1.0 1.0 1.0 1.0 丰都村 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.73 祥华乡 新寨村1 0.1 1.0 1.0 1.0 1.0 0.208 1.0 1.0 新寨村2 1.0 1.0 1.0 1.0 1.0 0.982 1.0 1.0 新寨村3 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.865 龙涓乡 内灶大桥头 1.0 1.0 1.0 1.0 1.0 0.532 1.0 1.0 楼梯垵 0.1 1.0 1.0 1.0 1.0 1.0 1.0 1.0 水云波 0.28 1.0 0.55 0.1 1.0 1.0 1.0 0.37 感德镇 槐植村 0.1 1.0 1.0 1.0 0.46 1.0 1.0 0.3475 福德村 1.0 1.0 1.0 1.0 0.82 1.0 1.0 0.1 霞云村 0.1 1.0 0.55 0.1 1.0 1.0 1.0 1.0 长坑乡 南斗村1 0.1 1.0 1.0 1.0 0.64 1.0 1.0 0.1 南斗村2 0.1 1.0 1.0 1.0 1.0 0.604 0.82 0.1 蓝田乡 后清村1 0.1 1.0 1.0 1.0 0.46 0.496 1.0 0.4375 后清村2 0.28 1.0 1.0 1.0 1.0 1.0 1.0 1.0 芦田镇 芦田村1 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 芦田村2 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.2575 剑斗镇 御屏村 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.8425 古石庵 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.1 柿树仔 0.64 1.0 1.0 1.0 0.82 1.0 1.0 0.3925 湖上乡 湖上乡1 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.4825 湖上乡2 0.82 1.0 1.0 1.0 1.0 1.0 0.46 1.0 湖上乡3 0.28 1.0 0.1 1.0 1.0 1.0 1.0 0.19 西坪镇 盖竹村 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 朱杨村 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.4375 宝山村 0.1 1.0 1.0 0.1 1.0 0.604 0.82 1.0 大坪乡 萍州村 0.46 1.0 1.0 1.0 1.0 1.0 1.0 1.0 福美村 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.1 双美村 0.1 1.0 0.1 1.0 1.0 0.1 1.0 0.1 虎邱镇 双格村 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.1 林东村 0.1 1.0 1.0 1.0 0.46 1.0 1.0 1.0 高村村 0.1 1.0 1.0 0.55 1.0 1.0 1.0 0.64 不同山1 1.0 1.0 0.1 0.55 1.0 1.0 1.0 1.0 不同山2 1.0 1.0 1.0 1.0 1.0 0.496 1.0 1.0 采用主成分分析计算相关矩阵，如表5，pH与AK呈显著正相关，证明pH提高可以使有效钾含量显著升高，即需要适宜pH来提高钾素利用率，TK与OM呈显著负相关，缺钾是限制有机质提升的重要原因。TP与AP间显著正相关，说明铁观音茶树对磷肥的吸收较好，OM与TN、TP呈显著正相关且TN与TP显著正相关，表明土壤有机质含量提升有助于提高土壤养分库容，土壤各养分间关系密切。
表 5 安溪茶园土壤肥力指标间相关矩阵Table 5. Correlation matrix between soil fertility indices of Anxi tea garden指标
IndexpH 有机质
OM全氮
TN全磷
TP全钾
TK碱解氮
AN速效磷
AP速效钾
AKpH 1.000 −0.088 −0.162 0.069 0.011 0.089 0.065 0.373* 有机质 −0.088 1.000 0.779** 0.322* −0.477** 0.019 −0.050 0.049 全氮 −0.162 0.779** 1.000 0.499** −0.309 0.272 0.213 0.063 全磷 0.069 0.322* 0.499** 1.000 −0.241 0.224 0.689** 0.150 全钾 0.011 −0.477** −0.309 −0.241 1.000 0.045 0.030 0.046 碱解氮 −0.089 0.019 0.272 0.224 −0.045 1.000 0.145 0.204 速效磷 0.065 −0.050 0.213 0.689** 0.030 0.145 1.000 0.119 速效钾 0.373* 0.049 0.063 0.150 0.046 0.204 0.119 1.000 注：“**”表示极显著相关（P < 0.01），“*”表示显著相关（P < 0.05）。 采用主成分分析法选取累计百分率大于80%的主成分（表6）来确定肥力因子的权重系数（表7）。
表 6 各肥力指标主成分总方差解释Table 6. Explanation of total variance of principal components of each fertility index成分
Component初始特征值
Initial eigenvalue提取载荷平方和
Extract sum of squares of load总计Total 方差百分比
Variance percentage累积( %)
Cumulative percentage总计
Total方差百分比
Variance percentage累积( %)
Cumulative percentage
1 2.564 32.055 32.055 2.564 32.055 32.055 2 1.628 20.353 52.408 1.628 20.353 52.408 3 1.219 15.240 67.648 1.219 15.240 67.648 4 1.014 12.678 80.326 1.014 12.678 80.326 5 0.715 8.943 89.269 6 0.503 6.290 95.559 7 0.222 2.776 98.335 8 0.133 1.665 100.000 表 7 土壤各养分指标权重系数Table 7. Weight coefficients of soil nutrient indexes指标
IndexpH 有机质
OM全氮
TN全磷
TP全钾
TK碱解氮
AN速效磷
AP速效钾
AK权重系数 0.085 0.113 0.171 0.197 0.001 0.126 0.161 0.146 如表6，选取累计百分率大于80%的前4个主成分，如表7，各养分权重系数中，全氮，全磷、有效磷占权重相对较高，说明磷素与氮素含量及利用率均较高，而pH与全钾权重较低，即对最终肥力贡献较低，说明酸化与钾素流失较为严重。权重系数结果表明，影响土壤最终肥力的因素主要为氮肥与磷肥的施用。
2.2.2 茶园土壤肥力综合评价
根据1.4计算安溪各乡镇铁观音茶园土壤肥力综合指数IFI（表8），IFI值在0.52 ~ 1.00之间，平均值为0.84，符合 Ⅰ 级优质茶园标准，变异系数为13.33%，属于弱变异，说明各茶园土壤肥力普遍较优良。其中有32个茶园土壤达到 Ⅰ 级标准，6个样点土壤达到 Ⅱ 级标准，优良率达100%，安溪铁观音茶园土壤适宜种植铁观音茶树。
表 8 铁观音茶园土壤肥力综合指数IFITable 8. Comprehensive index of soil fertility IFI of Tieguanyin tea garden乡镇
Township村落/地点
PointIFI 乡镇
Township村落/地点
PointIFI 桃舟乡 吾培村 0.89 芦田镇 芦田村1 1.00 产业园 0.92 芦田村2 0.81 莲山村 0.91 剑斗镇 御屏村 0.90 福田乡 双按村 0.79 古石庵 0.79 尾洋村 0.92 柿树仔 0.88 丰都村 0.88 湖上乡 湖上乡1 0.85 祥华乡 新寨村1 0.82 湖上乡2 0.90 新寨村2 1.00 湖上乡3 0.67 新寨村3 0.90 西坪镇 盖竹村 1.00 龙涓乡 内灶大桥头 0.94 朱杨村 0.84 楼梯垵 0.92 宝山村 0.67 水云波 0.59 大坪乡 萍州村 0.95 感德镇 槐植村 0.83 福美村 0.79 福德村 0.87 双美村 0.52 霞云村 0.67 虎邱镇 双格村 0.79 长坑乡 南斗村1 0.79 林东村 0.92 南斗村2 0.71 高村村 0.78 蓝田乡 后清村1 0.78 不同山1 0.76 后清村2 0.94 不同山2 0.94 使用克里金插值法，预测安溪县域内铁观音茶园土壤肥力综合指数IFI值空间分布。如图3，采样区域乡镇茶园土壤IFI值较高，土壤肥力较高；西部内安溪区域土壤肥力相对东部外安溪区域土壤肥力更高，且安溪西部边界处土壤更肥沃。内安溪是铁观音茶叶的主产区，因其独特的气候特点，铁观音茶叶品质也相对更好，土壤肥力综合指数IFI值空间分布图也证实了这点。
![]()
图 3 茶园土壤肥力综合指数IFI值空间分布图Figure 3. Spatial distribution map of IFI value of comprehensive index of soil fertility in the tea gardens3. 讨论
3.1 茶园土壤养分状况与肥料选择
茶树适宜生长于在酸性土壤环境中，但长期茶树种植会使茶园土壤酸化严重[14]。一般来说，适宜茶树生长pH值为4.5 ~ 5.5，在过去几十年间，茶园土壤酸化程度远大于种植水果、蔬菜和谷物土壤[14]。从近些年研究结果上看[3,15-18]，茶园土壤pH一直显强酸性（pH < 4.5），明显低于周边土壤，并且呈现出随时间递减的趋势。强酸性土壤占比（尤其pH < 4.0的茶园土壤）越来越高，茶园土壤过酸已经是一个普遍性问题。
造成这种现象的可能原因之一是土壤的自然酸化，即内源酸化因子，铁观音茶树自身物质循环与根系分泌物在土壤中的运动导致了茶园土壤的自然酸化[17]。另一个重要原因，即外源酸化因子，是无机氮肥的过度使用[19]。茶树是喜铵植物，茶树生长幼芽需要更多的氮，茶园中经常施用氮肥来提高铁观音茶叶的产量与品质[20]，而铵态氮肥增加了氮的淋失率，导致了土壤酸化。因此，长期的单一茶树种植和不合理施肥可能是造成茶园土壤酸化严重的重要原因。
从养分含量结果中可以看到，所有铁观音茶园土壤有机质含量均处于最适宜水平，与前人测量结果相比较[3-6,15-17, 21]，养分状况也有了较明显提升，但也不能忽略pH过低的问题，土壤长期过酸对铁观音茶园土壤、茶树状况及茶叶品质都会产生反作用。施肥是茶叶采前处理的重要手段，但是单一的施肥不一定可以按比例提高茶叶的产量与质量，反而可能加剧如土壤酸化、养分淋失等问题，肥料的选择至关重要。
近年来茶树种植者们响应国家号召，使用有机肥替代无机肥，取得了较好的成效，土壤中有机质含量大幅提升[7-9]，安溪铁观音茶园土壤养分状况得到了一定的改善，事实证明有机肥适合长期施用。研究发现，生态友好型肥料或土壤改良剂对茶园土壤改善具有较大潜力。例如，施用腐殖基质可以提高养分的有效性[22]，施用饼肥是既可以保持土壤肥力，又可以减少农业上面源污染[19]。复合肥与尿素按比例调和施用不仅能提高茶叶浸出物与茶汤浓度，也能改善土壤环境[23]。生物炭作为一种土壤改良剂，可以提高有效氮磷钾含量与pH值，提高茶叶叶面积和光合作用效率，促进植物的生长[23-25]。生物炭与肥料相结合，即生物炭基肥料，可以改善低养分酸性土壤[26]。土壤改良剂与菌剂联合使用比单独使用更能抑制土壤酸化[27]，将蜡样芽胞杆菌这类根际促生菌接种到制备的生物炭中，可以得到高效的生物肥料[28]。
目前，还需长期的田间试验、国家的补贴政策等，研究出更高效有益的肥料来改善土壤品质，实现茶叶品质与经济共增，助力乡村振兴，实现“三农梦”。
3.2 茶园土壤肥力状况与茶园管理
土壤肥力综合评价可以去除主观因素影响，体现土壤肥力综合水平。安溪县茶园土壤肥力综合评价结果表明，铁观音茶园土壤总体处于优良状态。整体来说，安溪茶园土壤肥力有了较高提升[3,6,16]。
监测茶园土壤肥力状况并提出适宜改良方案可以促进绿色茶园的发展。因此，我们要注重养分管理[29]，改善施肥方式，改善土壤养分循环，保持土壤的结构与肥力[30]，综合养分管理是一种强有力的策略[31]。通过建立养分专家系统（Nutrient Expert system，NE），既能减少化肥使用量，提高其利用率，又能降低温室气体排放，降低环境风险。对于已经不适宜种植的地区，实施“退茶还林”，转向原生态。密集耕作会导致养分大量流失，疏植会减少这种流失[32]。温室条件下修剪枯枝落叶回放到茶园土壤上，能提高生物量，提高茶树的株高[22]。在茶树间间作或套种大豆、花生等绿肥有利于提高pH、有机质、有效碳氮磷钾、土壤酶活性与资源利用率[33]。通过疏植留高、茶草共生、茶林混合等方式进行合理有效的生态茶园管理[34]，可以促进铁观音茶园生态系统可持续发展，管理后的茶园也会具有一定的观赏价值，未来或可实现茶园“公园化”。
4. 结论
安溪县13个乡镇38个铁观音茶园土壤pH为3.3 ~ 5.5，超过半数铁观音茶园的土壤已酸化严重；茶园土壤有机质、全氮、全磷、全钾、碱解氮、有效磷、速效钾整体含量丰富，平均值均达到 Ⅰ 级优质标准，各个养分指标含量存在着一定的空间差异；绝大多数茶园土壤肥力综合指数IFI值达到了 Ⅰ 级标准，安溪茶园土壤总体适宜种植铁观音茶树。
-
![]()
图 2 茶园土壤养分状况含量图
茶园土壤养分状况含量图分为点状图与柱状图。点状图用大、中、小三个点表示，分别代表 Ⅰ 级、Ⅱ 级、Ⅲ 级标准，柱状图用三条柱体表示，其中中间柱体代表测量平均值，用Av表示，两边柱体为固定值，代表了分级标准中的两个临界值。如果中间柱体高于右边柱体，代表该处土壤处于 Ⅰ 级标准，低于左侧则代表该处土壤处于 Ⅲ 级标准，处于两侧之间则代表该处土壤处于 Ⅱ 级标准。
Figure 2. Content map of soil nutrients in tea garden
![]()
图 3 茶园土壤肥力综合指数IFI值空间分布图
Figure 3. Spatial distribution map of IFI value of comprehensive index of soil fertility in the tea gardens
表 1 茶园土壤养分分级标准
Table 1 Classification standard of soil nutrients in tea garden
采样点
Sampling pointⅠ Ⅱ Ⅲ pH 4.5 ~ 5.5 4.0 ~ 4.5或5.5 ~ 6.5 < 4.0或 > 6.5 有机质 （g kg–1） > 15 10 ~ 15 < 10 全氮（g kg–1） > 1.0 0.8 ~ 1.0 < 0.8 全磷（g kg–1） > 0.6 0.4 ~ 0.6 < 0.4 全钾（g kg–1） > 10 5 ~ 10 < 5 碱解氮（mg kg–1） > 100 50 ~ 100 < 50 有效磷（mg kg–1） > 10 5 ~ 10 < 5 速效钾（mg kg–1） > 120 80 ~ 120 < 80 
下载: 导出CSV
表 2 各肥力指标的隶属度函数曲线转折点的取值
Table 2 Value of turning point of membership function curve of each fertility index
转折点取值
Value of turning pointpH 有机质
OM全氮
TN全磷
TP全钾
TK碱解氮
AN有效磷
AP速效钾
AKx1 4.0 10 0.8 0.4 5 50 5 80 x2 4.5 15 1.0 0.6 10 100 10 120 x3 5.5 − − − − − − − x4 6.5 − − − − − − −
下载: 导出CSV
表 3 土壤肥力综合指数划分表
Table 3 Division of comprehensive index of soil fertility
等级
LevelⅠ Ⅱ Ⅲ IV 土壤肥力综合指数IFI IFI ≥ 0.75 0.50 ≤ IFI ≤ 0.75 0.25≤IFI < 0.50 IFI < 0.25
下载: 导出CSV
表 4 安溪茶园土壤肥力指标的隶属度
Table 4 Membership degree of soil fertility index in Anxi tea garden
乡镇
Township村落/地点
PointpH 有机质
OM全氮
TN全磷
TP全钾
TK碱解氮
AN速效磷
AP速效钾
AK桃舟乡 吾培村 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.775 产业园 0.1 1.0 1.0 1.0 0.46 1.0 1.0 1.0 莲山村 0.1 1.0 1.0 1.0 0.46 1.0 1.0 0.91 福田乡 双按村 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.1 尾洋村 0.1 1.0 1.0 1.0 1.0 1.0 1.0 1.0 丰都村 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.73 祥华乡 新寨村1 0.1 1.0 1.0 1.0 1.0 0.208 1.0 1.0 新寨村2 1.0 1.0 1.0 1.0 1.0 0.982 1.0 1.0 新寨村3 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.865 龙涓乡 内灶大桥头 1.0 1.0 1.0 1.0 1.0 0.532 1.0 1.0 楼梯垵 0.1 1.0 1.0 1.0 1.0 1.0 1.0 1.0 水云波 0.28 1.0 0.55 0.1 1.0 1.0 1.0 0.37 感德镇 槐植村 0.1 1.0 1.0 1.0 0.46 1.0 1.0 0.3475 福德村 1.0 1.0 1.0 1.0 0.82 1.0 1.0 0.1 霞云村 0.1 1.0 0.55 0.1 1.0 1.0 1.0 1.0 长坑乡 南斗村1 0.1 1.0 1.0 1.0 0.64 1.0 1.0 0.1 南斗村2 0.1 1.0 1.0 1.0 1.0 0.604 0.82 0.1 蓝田乡 后清村1 0.1 1.0 1.0 1.0 0.46 0.496 1.0 0.4375 后清村2 0.28 1.0 1.0 1.0 1.0 1.0 1.0 1.0 芦田镇 芦田村1 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 芦田村2 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.2575 剑斗镇 御屏村 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.8425 古石庵 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.1 柿树仔 0.64 1.0 1.0 1.0 0.82 1.0 1.0 0.3925 湖上乡 湖上乡1 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.4825 湖上乡2 0.82 1.0 1.0 1.0 1.0 1.0 0.46 1.0 湖上乡3 0.28 1.0 0.1 1.0 1.0 1.0 1.0 0.19 西坪镇 盖竹村 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 朱杨村 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.4375 宝山村 0.1 1.0 1.0 0.1 1.0 0.604 0.82 1.0 大坪乡 萍州村 0.46 1.0 1.0 1.0 1.0 1.0 1.0 1.0 福美村 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.1 双美村 0.1 1.0 0.1 1.0 1.0 0.1 1.0 0.1 虎邱镇 双格村 0.1 1.0 1.0 1.0 1.0 1.0 1.0 0.1 林东村 0.1 1.0 1.0 1.0 0.46 1.0 1.0 1.0 高村村 0.1 1.0 1.0 0.55 1.0 1.0 1.0 0.64 不同山1 1.0 1.0 0.1 0.55 1.0 1.0 1.0 1.0 不同山2 1.0 1.0 1.0 1.0 1.0 0.496 1.0 1.0
下载: 导出CSV
表 5 安溪茶园土壤肥力指标间相关矩阵
Table 5 Correlation matrix between soil fertility indices of Anxi tea garden
指标
IndexpH 有机质
OM全氮
TN全磷
TP全钾
TK碱解氮
AN速效磷
AP速效钾
AKpH 1.000 −0.088 −0.162 0.069 0.011 0.089 0.065 0.373* 有机质 −0.088 1.000 0.779** 0.322* −0.477** 0.019 −0.050 0.049 全氮 −0.162 0.779** 1.000 0.499** −0.309 0.272 0.213 0.063 全磷 0.069 0.322* 0.499** 1.000 −0.241 0.224 0.689** 0.150 全钾 0.011 −0.477** −0.309 −0.241 1.000 0.045 0.030 0.046 碱解氮 −0.089 0.019 0.272 0.224 −0.045 1.000 0.145 0.204 速效磷 0.065 −0.050 0.213 0.689** 0.030 0.145 1.000 0.119 速效钾 0.373* 0.049 0.063 0.150 0.046 0.204 0.119 1.000 注：“**”表示极显著相关（P < 0.01），“*”表示显著相关（P < 0.05）。
下载: 导出CSV
表 6 各肥力指标主成分总方差解释
Table 6 Explanation of total variance of principal components of each fertility index
成分
Component初始特征值
Initial eigenvalue提取载荷平方和
Extract sum of squares of load总计Total 方差百分比
Variance percentage累积( %)
Cumulative percentage总计
Total方差百分比
Variance percentage累积( %)
Cumulative percentage
1 2.564 32.055 32.055 2.564 32.055 32.055 2 1.628 20.353 52.408 1.628 20.353 52.408 3 1.219 15.240 67.648 1.219 15.240 67.648 4 1.014 12.678 80.326 1.014 12.678 80.326 5 0.715 8.943 89.269 6 0.503 6.290 95.559 7 0.222 2.776 98.335 8 0.133 1.665 100.000
下载: 导出CSV
表 7 土壤各养分指标权重系数
Table 7 Weight coefficients of soil nutrient indexes
指标
IndexpH 有机质
OM全氮
TN全磷
TP全钾
TK碱解氮
AN速效磷
AP速效钾
AK权重系数 0.085 0.113 0.171 0.197 0.001 0.126 0.161 0.146
下载: 导出CSV
表 8 铁观音茶园土壤肥力综合指数IFI
Table 8 Comprehensive index of soil fertility IFI of Tieguanyin tea garden
乡镇
Township村落/地点
PointIFI 乡镇
Township村落/地点
PointIFI 桃舟乡 吾培村 0.89 芦田镇 芦田村1 1.00 产业园 0.92 芦田村2 0.81 莲山村 0.91 剑斗镇 御屏村 0.90 福田乡 双按村 0.79 古石庵 0.79 尾洋村 0.92 柿树仔 0.88 丰都村 0.88 湖上乡 湖上乡1 0.85 祥华乡 新寨村1 0.82 湖上乡2 0.90 新寨村2 1.00 湖上乡3 0.67 新寨村3 0.90 西坪镇 盖竹村 1.00 龙涓乡 内灶大桥头 0.94 朱杨村 0.84 楼梯垵 0.92 宝山村 0.67 水云波 0.59 大坪乡 萍州村 0.95 感德镇 槐植村 0.83 福美村 0.79 福德村 0.87 双美村 0.52 霞云村 0.67 虎邱镇 双格村 0.79 长坑乡 南斗村1 0.79 林东村 0.92 南斗村2 0.71 高村村 0.78 蓝田乡 后清村1 0.78 不同山1 0.76 后清村2 0.94 不同山2 0.94
下载: 导出CSV
-
[1] Karak T, Bhanat R M. Trace elements in tea leaves, made tea and tea infusion: A review[J], Food Research International, 2010, 43(9): 2234 ~ 2252.
[2] 黄华斌, 林承奇, 于瑞莲, 等. 安溪铁观音茶园土壤重金属分布及污染评价[J]. 环境化学, 2018, 37(5): 994 − 1001. [3] 杨如兴, 张 磊, 王文建, 等. 安溪铁观音茶园土壤肥力分析[J]. 中国农学通报, 2010, 26(21): 160 − 166. [4] 张阿惠. 茶园土壤肥力综合评价及其绿色可持续发展思考[J]. 低碳世界, 2021, 11(1): 213 − 215. doi: 10.3969/j.issn.2095-2066.2021.01.106 [5] 叶伟娇. 安溪县茶园质量等级评价指标体系构建研究[D]. 福建农林大学, 2020. [6] 李 娟, 章明清, 尤志明, 章赞德, 姚建族, 苏火贵, 刘辽源. 福建乌龙茶园土壤速效氮磷钾丰缺指标与肥力状况评价[J]. 茶叶学报, 2018, 59(01): 19 − 25. doi: 10.3969/j.issn.1007-4872.2018.01.005 [7] 耿赛攀, 马立锋, 杨向德, 等. 有机肥替代化肥对茶园土壤质量及新梢养分吸收的影响[J]. 中国茶叶, 2021, 43(5): 52 − 57. doi: 10.3969/j.issn.1000-3150.2021.05.008 [8] 郭 龙, 李 陈, 刘佩诗, 等. 牛粪有机肥替代化肥对茶叶产量、品质及茶园土壤肥力的影响[J]. 水土保持学报, 2021, 35(6): 264 − 269. [9] 龚雪蛟, 秦 琳, 刘 飞, 等. 有机类肥料对土壤养分含量的影响[J]. 应用生态学报, 2020, 31(4): 1403 − 1416. doi: 10.13287/j.1001-9332.202004.025 [10] 任艳芳, 何俊瑜, 张艳超, 等. 贵州省开阳茶园土壤养分状况与肥力质量评价[J]. 土壤, 2016, 48(4): 668 − 674. [11] 柳书俊, 姚新转, 赵德刚, 等. 湄潭茶园土壤养分特征及肥力质量评价[J]. 草业学报, 2020, 29(11): 33 − 45. doi: 10.11686/cyxb2020061 [12] 刘 娟, 张乃明, 邓 洪. 勐海县茶园土壤养分状况及肥力质量评价[J]. 农业资源与环境学报, 2021, 38(1): 79 − 86. [13] 王育军, 周冀衡, 孙书斌, 等. 云南省罗平县烟区土壤肥力适宜性评价及养分时空变异特征[J]. 土壤, 2015, 47(3): 515 − 523. [14] Yan P, Wu L Q, Wang D H, et al. Soil acidification in Chinese tea plantations[J]. Science of The Total Environment, 2020, 715: e136963. doi: 10.1016/j.scitotenv.2020.136963
[15] 张炳铃. 安溪县茶园土壤养分特性评价[J]. 福建农业科技, 2012, (7): 40 − 42. doi: 10.3969/j.issn.0253-2301.2012.07.023 [16] 苏火贵. 安溪县茶园肥力状况与合理利用研究[J]. 中国农学通报, 2015, 31(28): 132 − 135. doi: 10.11924/j.issn.1000-6850.casb15040049 [17] 穆 聪, 黄梓璨, 苏 达, 等. 铁观音茶园土壤酸化与交换性镁的现状[J]. 中国土壤与肥料, 2021, (2): 39 − 43. doi: 10.11838/sfsc.1673-6257.20061 [18] 林 诚, 陈子聪, 吴一群, 等. 林地转变为茶园的土壤pH及养分变化特征[J]. 茶叶科学, 2020, 40(2): 186 − 193. doi: 10.3969/j.issn.1000-369X.2020.02.005 [19] 王海斌, 陈晓婷, 丁 力, 等. 福建省安溪县茶园土壤酸化对茶树产量及品质的影响[J]. 应用与环境生物学报, 2018, 24(6): 1398 − 1403. [20] Xie S W, Yang F, Feng H X, et al. Organic fertilizer reduced carbon and nitrogen in runoff and buffered soil acidification in tea plantations: Evidence in nutrient contents and isotope fractionations[J]. Science of The Total Environment, 2021, 762: e143059. doi: 10.1016/j.scitotenv.2020.143059
[21] 吴志丹, 江福英, 张 磊, 等. 福建省安溪县铁观音茶园土壤氮素状况[J]. 土壤, 2020, 52(1): 16 − 24. [22] Pramanik P, Safique S, Jahan A. Humic substrates facilitated nitrogen uptake in tea bushes by influencing biochemical and microbiological properties in soil of Northeast India[J]. Journal of Plant Nutrition, 2018, 41(18): 2376 − 2385. doi: 10.1080/01904167.2018.1510510
[23] Gao S L, Hu S S, He P, et al. Effects of reducing chemical fertilizer on the quality components of Tieguanyin tea leaves[J]. IOP Conference Series Earth and Environmental Science, 2020, 559: e012020. doi: 10.1088/1755~1315/559/1/012020
[24] Wang Z T, Geng Y B, Liang T. Optimization of reduced chemical fertilizer use in tea gardens based on the assessment of related environmental and economic benefits[J]. Science of The Total Environment, 713: e136439. doi: 10.1016/j.scitotenv.2019.136439.
[25] Yan P, Shen C, Zou Z H, et al. Biochar stimulates tea growth by improving nutrients in acidic soil[J]. Scientia Horticulturae, 2021, 283: e110078. doi: 10.1016/j.scienta.2021.110078
[26] 金桂梅, 李昱航, 郑向群, 等. 不同土壤管理与施肥模式对茶园土壤环境及茶叶产量的影响[J]. 土壤通报, 2020, 51(1): 152 − 158. doi: 10.19336/j.cnki.trtb.2020.01.20 [27] Wang L, Butterly C R, Chen Q H, et al. Surface Amendments Can Ameliorate Subsoil Acidity in Tea Garden Soils of High ~ Rainfall Environments[J]. Pedosphere, 2016, 26(2): 180 − 191. doi: 10.1016/S1002-0160(15)60031-5
[28] Azeem M, Hassan T U, Tahir M I, et al. Tea leaves biochar as a carrier of Bacillus cereus improves the soil function and crop productivity[J]. Applied Soil Ecology, 2021, 157: e103732. doi: 10.1016/j.apsoil.2020.103732
[29] 刘 扬, 孙丽莉, 廖 红. 养分管理对安溪茶园土壤肥力及茶叶品质的影响[J]. 土壤学报, 2020, 57(4): 917 − 927. [30] He S Q, Zheng Z C, Zhu R H. Long ~ term tea plantation effects on composition and stabilization of soil organic matter in Southwest China[J]. Catena, 2021: 199. doi: 10.1016/j.catena.2020.105132
[31] Tang S, Zheng N, Ma Q X, et al. Applying Nutrient Expert system for rational fertilisation to tea (Camellia sinensis) reduces environmental risks and increases economic benefits[J]. Journal of Cleaner Production, 2021, 305: e127197. doi: 10.1016/j.jclepro.2021.127197
[32] Chu L, Hennayake H M K D, Sun H. Biochar Effectively Reduces Ammonia Volatilization from Nitrogen ~ Applied Soils in Tea and Bamboo Plantations[J]. Phyton, 2019, 88(3): 261 − 267. doi: 10.32604/phyton.2019.07791
[33] Ma Y H, Fu S L, Zhang X P, et al. Intercropping improves soil nutrient availability, soil enzyme activity and tea quantity and quality[J]. Applied Soil Ecology, 2017, 119: 171 − 178. doi: 10.1016/j.apsoil.2017.06.028
[34] 刘金龙, 高水练. 生态茶园管理模式探索与实践−以安溪县举源茶叶专业合作社为例[J]. 中国茶叶加工, 2018, (2): 28 − 30. doi: 10.15905/j.cnki.33-1157/ts.2018.02.010 -
期刊类型引用(9)
1. 龙亚芹，罗梓文，陈洪云，陈龙，曲浩，龙丽雪，陈强，杨向德，陈林波. 福贡县茶园土壤养分状况分析与评价. 中国茶叶. 2025(02): 54-60 . 
百度学术
2. 杨兰珍，陈辉，唐越岳，赵福林，何腾兵. 基于GIS和模糊数学法的凤冈县田坝村茶园土壤肥力质量评价. 扬州大学学报(农业与生命科学版). 2025(01): 100-109 . 百度学术
3. 吕思航. 茶产业可持续发展助力乡村振兴对策探究——以福建省安溪县为例. 广东蚕业. 2024(02): 115-117 . 百度学术
4. 庞英华，顾万帆，金鑫，张舟娜，金雅慧，郑洪福. 杭州余杭径山茶园土壤养分状况及肥力质量评价. 中国农学通报. 2024(09): 97-105 . 百度学术
5. 郭灿，段学艺，高秀兵，曹雨，匡模，王校常. 湄潭县茶园土壤主要养分变化特征浅析. 茶叶. 2024(04): 221-226 . 百度学术
6. 王育平. 安溪县北部茶园土壤钾素含量分布状况研究. 茶叶学报. 2024(06): 48-54 . 百度学术
7. 谢辉煌. 植物工厂生产方式下铁观音茶树栽培技术. 中国茶叶加工. 2023(03): 69-73 . 百度学术
8. 高丽云. 安溪铁观音主产区茶园土壤环境条件与茶叶质量提升建议. 福建茶叶. 2023(11): 30-33 . 百度学术
9. 沈玥，李倩怡，赵文慧，王桂良，钱晓晴. 东海县农田土壤养分状况与肥力评价. 河南农业科学. 2023(12): 77-87 . 百度学术
其他类型引用(5)
计量
- 文章访问数: 187
- HTML全文浏览量: 37
- PDF下载量: 61
- 被引次数: 14
