茶(Camellia sinensis(L.)O. Kuntze)是世界上最普遍的非酒精饮料作物，属于山茶科山茶属，其中，中国种(Camellia sinensis var. sinensis)包括云南亚种和武夷亚种^[1-3]。岩茶是武夷山最具代表性的茶类，属于半发酵型的闽北乌龙茶^[4]。武夷岩茶的品质具有明显的地域特性(俗称：山场)，民间根据茶叶品质的高低，大致将武夷山分为三个茶区，包括正岩区、半岩区和洲茶区^[5]。正岩区主要是指武夷山景区内“三坑两涧”的茶园，品质最高；半岩区为九曲溪内、“三坑两涧”外的茶园，品质较好；洲茶区则为景区外的茶园，品质一般^[5]。有研究利用超高效液相色谱-四级杆-飞行时间质谱(UPLC-Q-TOF-MS)分析方法，揭示了武夷岩茶的品质是一系列代谢物的综合效果，并且与产地的密切关系客观存在^[6]。

不同茶区的土壤成土母质差别很大，其矿质养分含量也不同^[7]，唐代陆羽在《茶经》中提到，不同类型的土壤上种植出来的茶叶等级不同，表明土壤养分可能影响茶叶的品质。通过对福建省茶叶主产区武夷山和安溪茶园进行调研分析发现，茶园的土壤肥力显著影响茶叶的品质，肥力不足或过量均会导致茶叶的品质下降^[8-9]。土壤中氮过高会显著降低茶园土壤的酶活性，不利于茶树养分的吸收转化，从而导致茶叶代谢物的累积减少^[10]。此外，土壤中的微量元素也会影响茶叶的生化成分^[11]，通过外施微量元素可调控茶叶的品质，铜和锌肥均可提高茶叶中茶多酚、咖啡碱和游离氨基酸等主要代谢物的含量^[12]，施加适量锰能够降低茶叶中茶多酚的含量^[13]。以上研究表明，茶园土壤的养分状况显著影响茶叶的品质。

同时，茶青中的矿质养分影响代谢物的合成和转化，进而决定茶叶的品质。通过大量分析安溪与武夷山茶青中代谢物与养分的关系，发现茶青中的氮磷钾存在一定的比例关系，养分过高或过低均不利于代谢物的积累，其中对茶氨酸浓度的影响最大^[14]。茶叶中的氮直接参与茶氨酸和儿茶素的代谢^[15-16]，钾影响着氨基酸与咖啡因的代谢^[17-18]，锰、钙、镁则与总糖存在一定关系^[19]。因此，茶青中的矿质元素对茶叶的品质形成具有重要影响。本研究根据武夷山的不同茶叶等级划分的传统区域(正岩区、半岩区、洲茶区)，选取代表性茶园，分析不同产地茶青中矿质元素的差异，寻找土壤肥力与茶青养分的关系，旨在为茶园土壤养分管理提供科学依据。

采样概况：本研究以武夷山主栽茶树品种肉桂为主要研究对象，采样时间为2018年春季，主要采集茶青样品和茶园对应的土壤样品。取样茶园共计112个，涉及武夷山正岩区、半岩区和洲茶区三类茶区，茶园数分别为41、21和50个。具体采样点包括正岩区：牛栏坑、慧苑坑、倒水坑、流香涧、悟源涧等；半岩区：岩茶村、乞丐岩、白云岩、一线天、燕子窠等；洲茶区：兴田、黄柏、仙店、历屯、太庙等(图 1)。

图 1 武夷山取样分布图 Fig. 1 Distribution of the sampling spots in the Wuyi Mountain

茶叶样品的处理：按照当地的采茶标准采摘三叶一新(新叶为全开面)，样品采集后24 h内采用热固样法固定，105 ℃烘箱加热30 min，75 ℃烘干至恒重，研磨成粉末后称样消解，用于测定12种主要矿质元素。每个茶园取10个茶青样品，共计1 120个茶青样品。

茶园土壤的处理：取土前，先清除土壤表层的杂草及落叶；取样时，在茶树滴水缘垂直向下、靠近茶树根系的部位，使用钻土器采集0~20 cm土层的土壤；取土后，土壤摊放于阴凉处风干、研磨过筛后，用于检测土壤理化性质及养分指标。每个茶园取3个土壤样品，共计336个土壤样品。

茶青养分检测^[20]：浓硫酸-双氧水消解法消解茶青获得待测溶液，使用流动分析仪(Skalar，荷兰)测定大量元素氮(N)和磷(P)浓度，使用火焰光度计(410型，Sherwood，英国)测定钾(K)浓度；硝酸消解法获得茶叶样品待测溶液，使用电感耦合等离子体质谱仪ICP-MS7900(Agilent，美国)检测钙(Ca)、镁(Mg)、铁(Fe)、锰(Mn)、硼(B)、锌(Zn)、铜(Cu)、钼(Mo)和镍(Ni)等9种中微量元素。

土壤肥力检测^[20]：土壤酸碱度(pH)测定采用电位法；有机质(OM)测定采用高温外热重铬酸钾氧化-容量法；碱解氮(AN)测定采用碱解扩散吸收法；有效磷(AP)测定采用钼锑抗比色法；速效钾(AK)测定采用乙酸铵浸提-火焰光度法。

采用SPSS 19.0对数据进行单因素方差分析(One-way ANOVA)并选用邓肯(Duncan)法进行多重比较，使用R语言中的Factoextra和FactoMineR程序包进行主成分分析，RandomForest程序包进行随机森林分析，运用ggplot2程序包绘制小提琴图，corrplot和RColorBrewer程序包进行皮尔森相关性分析及可视化。

对武夷山112个茶园采集的1 120个茶青样品进行12种主要矿质元素的测定分析，发现茶青中各矿质元素的含量差异很大(表 1)。大量元素中，氮的含量最高，为23.48 g·kg^-1；钾和磷次之，分别为14.41和3.13 g·kg^-1，氮：磷：钾的比例约为(7~8)：1：(4~5)，表明大量元素中磷的含量相对较低(表 1)。中量元素中，钙和镁的含量分别约为磷的82.0%和57.8%，表明茶树对钙和镁的需求量较大(表 1)。微量元素中，锰的含量最高，达到了551.77 mg·kg^-1，约为铁的4.3倍，但其变异系数也是所分析元素中最大的，为51.11%，表明茶树需要大量的锰，且其阈值范围较广(表 1)。

表 1 Table 1

表 1 茶青中矿质元素含量 Table 1 Contents of mineral elements in tea leaves 矿质元素 Mineral elements 平均值1) Mean 标准误 SE 变异系数 CV/% 氮N 23.48 0.22 9.96 磷P 3.13 0.05 16.58 钾K 14.41 0.11 7.63 钙Ca 2.57 0.07 28.70 镁Mg 1.81 0.02 11.47 铁Fe 127.48 2.82 22.81 锰Mn 551.77 27.26 51.11 硼B 17.34 0.59 35.13 锌Zn 70.01 1.13 16.69 铜Cu 10.90 0.32 30.66 钼Mo 0.87 0.07 84.70 镍Ni 7.70 0.15 20.28 注：1)茶青中氮、磷、钾、钙和镁含量的单位为g·kg-1，其余元素含量单位为mg·kg-1。Note：The unit for content of N，P，K，Ca and Mg is g·kg-1，while for contents of the other elements are mg·kg-1.

表 1 茶青中矿质元素含量 Table 1 Contents of mineral elements in tea leaves

针对12种矿质元素含量进行主成分分析，仅前两个主成分就可解释茶青养分整体差异的52%(第一主成分和第二主成分可分别解释差异的38.2%和13.8%)(图 2)。基于第一主成分和第二主成分，不同产地之间的茶青样品差异较大，尤其是正岩区与洲茶区之间，仅存在少许叠加；而半岩区的茶青样品之间差异较大，明显分隔成两个部分，其中一部分完全与正岩重叠，另一部分则趋近于洲茶区(图 2)。上述结果可以说明正岩区与洲茶区的茶青养分含量差异明显，而半岩茶则介于正岩茶与洲茶之间。

图 2 茶青中矿质元素的主成分分析 Fig. 2 Principal component analysis(PCA)of the mineral elements in tea leaves

为了筛选形成产地差异的主要贡献元素，以产地作为分类依据进行随机森林分析，得到每个元素的平均精度下降值(mean decrease accuracy)和平均基尼指数降低值(mean decrease Gini)。结果显示磷的平均精度下降值和平均基尼指数降低值均为最高，铜和氮次之，表明对产地差异贡献最大的三种矿质元素依次为磷、铜和氮(图 3)。因此，在分析的12种矿质元素中，这三种元素可能对不同产地茶叶的品质差异起到主要贡献。

图 3 不同产地茶青中矿质元素的随机森林分析 Fig. 3 Random forest analysis of the mineral elements in tea leaves relative to area

同时，本研究发现磷、铜、氮含量在正岩区、半岩区以及洲茶区之间确实存在着显著差异(图 4)。正岩区茶青中这三种元素的含量总体低于洲茶，表明高品质的正岩茶中含有较低含量的磷、铜、氮；而半岩茶居于正岩茶与洲茶之间，半岩茶中的磷和氮含量与正岩茶无显著差异，铜则与洲茶无差异(图 4)。以上结果表明，茶青中过高的磷、铜、氮含量可能不利于高品质茶叶的形成。

注：不同字母表示不同产地间差异显著(P < 0.05)。  Note: Different letters indicate significant differences between areas at 0.05 level. 图 4 不同产地茶青中磷、铜、氮的含量 Fig. 4 Contents of P, Cu and N in the tea leaves relative to area

为了探究土壤肥力与茶青养分的关系，本研究将土壤肥力(pH、有机质、碱解氮、有效磷和速效钾)与茶青矿质元素含量(磷、铜和氮)进行皮尔森相关性分析。结果发现，土壤pH与有效磷、速效钾不相关，而与有机质、碱解氮呈很弱的负相关，表明茶园土壤酸碱度与土壤养分关系不明显；土壤速效钾分别与碱解氮、有效磷的关系较为密切，这可能与茶农施用等配比氮磷钾肥的习惯有关(图 5)。茶青中，磷与氮含量呈现较强的正相关，说明茶叶中磷和氮的吸收存在协同效应(图 5)。而就土壤与茶青的关系而言，相关系数的绝对值最高为0.11，表明二者几乎不相关，即土壤肥力对茶青中主要贡献元素磷、铜、氮的水平影响不明显(图 5)。

注：pH、OM、AN、AP、AK分别表示土壤酸碱度、有机质、碱解氮、有效磷和速效钾。  Note: pH, OM, AN, AP and AK stands for soil pH, organic matter, alkaline nitrogen, available phosphorus and available potassium, respectively. 0.01 < *P < 0.05;**P < 0.01. 图 5 土壤肥力与茶青中磷、铜和氮含量的关系 Fig. 5 Relationships between soil fertility and contents of P, Cu and N in the tea leaves

武夷山的茶叶根据产地所在区域划分为正岩区(高品质)、半岩区(中等品质)和洲茶区(一般品质)^[5]。本研究在三个不同品质区域分别选取代表性茶园(图 1)，分析不同产地茶青中矿质元素的差异。通过主成分分析方法，发现正岩区与洲茶区的茶青养分含量差异明显(图 2)；此外，利用随机森林分析方法中的平均精度下降值和平均基尼指数降低值作为判断指标，来反映变量对决策树节点纯度的影响。一般而言，两个指标的数值越大，代表该变量在分类中的重要性越大。从图 3结果可发现对产地差异贡献最大的三种矿质元素，其中就包括大量元素氮和磷。氮和磷是茶树生长发育过程中需求量较大的两个矿质元素。Sun等^[14]在对福建武夷山与安溪的茶叶样品进行养分测定后，认为茶青中氮磷比例范围为6:1~13:1，磷含量显著低于氮。本研究发现，武夷山肉桂茶青样品的总体氮磷比例较为恒定，约为7.5:1(表 1)，符合前期研究比例范围。进一步对武夷山不同品质茶区洲茶区、半岩区和正岩区的茶青样品进行分析后，得出这三个区域的氮磷比例分别为6.8:1、7.7:1和8.0:1，虽较为相近，但存在一定的趋势，随着品质的提高而逐渐增大，即高品质的茶青中磷与氮的比例较低；并且，正岩茶青中的氮和磷含量均显著低于洲茶茶青，分别为洲茶的86.8%与74.2%(图 4)。综合以上研究得出，高品质的茶青中磷与氮的含量及其比例均较低。

植物的养分主要通过根系从土壤中吸收获得，通常情况下，土壤养分状况能够影响植物体内的养分，比如，高磷处理能够分别提高柱花草与大豆中全磷或可溶性磷的含量^[21-22]。因此，本研究对土壤主要肥力指标与茶青中对产地差异贡献最大的三种矿质元素进行皮尔森相关性分析，结果发现，土壤肥力对茶青中主要贡献元素磷、铜、氮的水平影响并不明显(图 5)。周志等^[6]研究发现，相比洲茶区土壤，正岩区土壤为高磷低氮；同样地，早在20年前，姚月明^[23]认为正岩茶园土壤的磷较洲茶区高，而氮要低；然而，品质更优的名岩区土壤中的铜含量较品质一般的丹岩区低^[7]。本研究中正岩区茶青中的磷、铜、氮含量均低于洲茶区(图 4)，说明土壤肥力部分影响茶青的矿质元素组成，但不是决定因素。

武夷岩茶中的矿质养分与土壤养分的关系形成可能受两个因素的影响，一是茶树自身的养分需求，即不同茶区的茶树品种对矿质养分的需求存在差异。孙威江等^[24]研究发现，不同产区的武夷岩茶鲜叶中营养元素存在显著差异，这与本身的土壤营养环境之间并不存在很好的相关性。因此，植物对养分的吸收和需求以植物自身的响应为判断标准，土壤状况仅能作为参考。二是不同岩区的地理环境差异。武夷山作为中国著名的风景旅游区，具有典型的丹霞地貌特征^[25]。与洲茶区相比，位于武夷山核心地段的正岩区，具有降水量多，湿度大，雾日长等气候特点，这些因素将综合影响茶树中离子组和代谢组的变化，并最终形成其独特的口感和品质^[5]。

不同产地的茶青中12种主要矿质元素存在一定的差异，尤其是磷、铜和氮元素，对不同岩区茶青中矿质元素的整体差异贡献最大，且均在正岩区茶青中浓度最低，洲茶中最高，而土壤肥力5项指标(pH、碱解氮、有效磷、速效钾和有机质)与茶青养分的累积相关性较弱。