Printer Friendly

Stable Isotope Ratio and Elemental Profile Combined with Support Vector Machine for Provenance Discrimination of Oolong Tea (Wuyi-Rock Tea).

1. Introduction

Oolong tea is a traditional beverage favoured by consumers all over the world for its pleasurable aroma and taste. In addition, Oolong tea is a rich source of antioxidants, such as tea polyphenol and tea polysaccharide, so it is also reported as a functional drink that combats obesity, hypoglycaemia, and oral bacterial infection [1-4].

For tea products, their aroma and savor are influenced by many aspects, such as geographical origin [5,6], tea specie, cultivation, and processing method [7,8]. Among these aspects, the geographical and natural conditions in which the tea trees grow are widely perceived to be a key factor. Therefore, in China, the majority of famous teas are named for their provenance, such as the Anxi-Tieguanyin tea, the West Lake-Longjing tea, and the Anji-White tea.

Wuyi-Rock tea is originally cultivated in a mountain in the north of Fujian Province (Wuyi Mountain). Contributed to by the unique climate and edatope of Wuyi Mountain, Wuyi-Rock tea (WRT) is recognized as one of the most prestigious Oolong teas for its special savor and long-lasting fragrance. Therefore, Wuyi-Rock tea has been awarded a protected geographical indication (PGI) and exported to more than 30 countries. However, the actual yield of WRT is limited, and it cannot satisfy the needs of consumers. In various markets, many teas labeled as Wuyi-Rock tea were actually cultivated outside the protected production area; some of them are not even cultivated in Fujian Province. Although the taste and aroma of non-Wuyi-Rock teas (NWRT) are inferior to authentic WRT, teas planted in different geographical origins still have a similar appearance, and they can hardly be distinguished just by the naked eye. In traditional sensory analysis, WRT was tasted by professional tea tasters, and then counterfeits were identified based on a series of sensory scores. However, the result of sensory analysis depends a great deal on subjective decisions by tea tasters, and a well-trained tea taster is hard to find. Therefore, an urgent demand exists for developing a more effective and stable technique to discriminate the provenances of tea products.

In recent years, analytical methods based on instrumental technology have been widely applied in food quality control [9,10]. Isotope ratio mass spectrometry (IRMS) is a technique for measuring the isotope content, which is a highly indicative parameter in provenance discrimination. For example, the stable isotope ratio of hydrogen ([[delta].sup.2]H) in plants is influenced by the latitude and altitude of the production site. The concentration of deuterium in the water decreases when clouds form above the ocean. Then, as rainwater falls and the clouds move inland and gain altitude, the content of [[delta].sup.2]H in rainwater decreases gradually [11]. As a result, an isotopic gradient exists in groundwater from coast to inland. Moreover, the variation of oxygen-18 ([[delta].sup.18]O) follows the same pattern as hydrogen in the hydrosphere [12]. The isotope ratio of carbon ([[delta].sup.13]C) is strongly environmentally dependent; plants cultivated in humid environments have a lower [[delta].sup.13]C than plants in arid environments. The isotope ratio of nitrogen ([[delta].sup.15]N) is influenced by agricultural practices; plants treated with organic fertilizer develop a higher [[delta].sup.15]N content than plants treated with chemical fertilizer. For provenance analysis, IRMS has been successfully used for the analysis of orange juice, fruits, cow milk, and wine [13-16]. In addition to isotopes, trace element profiling is also available. There are various metallic elements in agricultural food; some of them are easily affected by edaphic and environmental factors, such as fertilization, soil type, climate, and temperature. Both inductively coupled plasma mass spectrometry (ICP-MS) and atomic spectroscopy are tools for quantitative determination of trace elements. In provenance discrimination, element profiling has been used for honey, onion, black tea, and wine [17-21].

This paper aims to develop an automatic analytical method for discriminating geographical origin of WRT by stable isotope and trace element contents. To model the complicated relationship between the measured data and production site, nonlinear multivariate classification models are usually used [22-25]. Compared with other nonlinear models such as kernel partial least squares (PLS) and artificial neural networks (ANNs), support vector machine (SVM) analysis is more suitable for the target data in this experiment because of the small sample size [26-28]. Consequently, the classification model was built based on SVM. Then, each variable of isotope and element data was ranked by its contribution to the model.

2. Materials and Methods

2.1. Tea Samples. In total, 99 authentic WRT samples were collected from 11 main rich-producing areas in Wuyishan, and 33 NWRT samples were collected from 11 different production sites. All of the samples were made of spring teas picked in 2015. Before analysis, the samples were preserved in cold and dry storage with lightproof packaging. Detailed information about the above samples is displayed in Table 1 and Figure 1.

The origin of Wuyi-Rock tea is the administrative area of Wuyishan city, and Wuyishan has 11 subdistricts. The range of sample collection covered all over the whole Wuyishan city area. So the number of samples is adequate. All the collected tea samples were obtained directly from the local tea processing space with the help of official department. All the tea samples belonged to "Wuyi-Rock tea," and the non-Wuyi-Rock tea samples were purchased outside the protected production area such as Jianyang, Jianou, and Ganzhou.

2.2. Isotopic Ratio Determinations. [[delta].sup.13]C, [[delta].sup.15]N, [[delta].sup.18]O, and [[delta].sup.2]H were measured using a MAT-253 isotope ratio mass spectrometer (Thermo Fisher, USA) connected to a Flash-2000 organic elemental analyser (Thermo Fisher, USA). During carbon and nitrogen isotope analysis, the quartz tube of the reactor was packed with chromium trioxide, high purity copper, and silver cobalt oxide to completely oxidize the organic matter. The carbon and nitrogen element carried by helium gas entered the IRMS in the form of C[O.sub.2] and [N.sub.2], respectively. The standard C[O.sub.2] and [N.sub.2] gases were used as reference gas before and after the isotope test of organic matter, and the detection of the instrument state in the sample analysis was completed by the standard sample such as labeled urea, IAEA-600, IAEA-CH-3, and VPDB (Vienna Pee Dee Belemnite). Similarly, the instrument state was detected by benzoic acid, IAEA-601, IAEA-602, IAEA-CH-7, and VSMOW (Vienna Standard Mean Ocean Water) in the analysis process of hydrogen and oxygen isotope. Each sample was repeated three times.

The measured values ([sup.13]C/[sup.12]C, [sup.15]N/[sup.14]N, [sup.18]O/[sup.16]O, and [sup.2]H/[sup.1]H) are usually presented as isotopic deviations, [delta], defined as follows:

[delta][per thousand] = ([R.sub.s]/[R.sub.std] - 1) x 1000, (1)

where [R.sub.s] is the measured value of the sample and [R.sub.std] is the measured value of an international standard. [R.sub.s] and [R.sub.std] are the ratios of the heavier isotope and lighter isotope of an element.

In isotope profiling, calibration was conducted according to calibrated-urea and calibrated-benzoic acid standards as well as the IAEA-600, IAEA-601, IAEA-602, IAEA-CH-3, and IAEA-CH-3 standards of the International Atomic Energy Agency (IAEA, Vienna).

2.3. Metal Determinations. Before metallic element detection, all of the tea samples were pretreated with microwave assisted digestion. The samples were dried before digestion process (placed in the oven for 4 hours at 80[degrees] C). All of the tea samples were manufactured in May 2015 and simultaneously analysed. Water content of tea samples was less than 6% before drying. 0.3 g of each dried tea sample was placed into a digestion vessel. Then, 1 mL ultrapure water and 5 mL nitric acid were added to the vessel. Finally, the digestion vessel was heated in the microwave cavity (2450 MHz). The vibration of gas pressure in the digestion procedure was conducted as follows: (1) ramping from normal pressure to 0.5 Mpa and holding at 0.5 Mpa for 70 s, (2) ramping to 1.0 Mpa for 50 s, (3) ramping to 1.5 Mpa for 50 s, and (4) ramping to 2.0 Mpa for 300 s. After the digestion, we evaporated the excess nitric acid in 130[degrees]C. When the temperature of the digestion vessel had cooled to room temperature, the digested liquid was moved to a volumetric flask and diluted to 50 mL with ultrapure water.

Subsequently, concentrations of 14 metallic elements were measured. The concentrations of Ti, Cr, Co, Ni, Cu, Zn, Rb, Cd, Cs, Ba, and Sr were detected using an X Series-IIICP-MS (Thermo Fisher, The USA). The information about the working condition and parameters for ICP-MS is shown in Table 2. The concentrations of Ca, Mg, and Mn were analysed by HITACHI 180-50 flame atomic absorption spectroscopy (FAAS, HITACHI, Japan). The main parameters for FAAS are presented in Table 3. The results of ICP-MS and FAAS were calibrated by mixed standard solution (GSB04-1767-2004) and biological component analysis standard substance, tea (GBW10052), respectively.

2.4. Data Splicing. All the data analysis was performed using MATLAB 7.14.0.739 (Mathworks, Sherborn, MA). For data splicing, the data of IRMS can be described as an n x p matrix A with n rows and p columns. n represents the number of samples and p is the number of features in this paper (n = 132, p = 4). In the same way, an n x q matrix B (in this paper, q = 14) was obtained from metallic element detection. Then, the columns of matrix B were arranged behind the last column of matrix A. As a result, a union matrix C (with n rows and p + q columns) was formed that contains both the isotope and metallic element information.

Before data analysis, each variable in matrix C was normalized as follows:

[x.sup.*] = x - [x.sub.min]/[x.sub.max] - [x.sub.min], (2)

where x is the value of ith (i = 1: 132) row and jth (j = 1: 18) column, [x.sub.max] is the max value in the jth column, and [x.sub.min] is the min value in the jth column.

2.5. SVM Analysis. The support vector machine algorithm is a type of classification and regression model for supervised machine learning. The kernel function is the main factor in the SVM algorithm. Kernels have the advantage of operating in the input space, where the solution of the classification problem is a weighted sum of kernel functions evaluated at the support vectors. SVM is designed to find an optimal plane that all the sample units can be divided into two classes in a multidimensional space. The optimal plane is in the middle of the nearest points between two classes and makes the distance as far as possible. For N variables, the optimal hyperplane is of N-1 dimensions [29].

After the data splicing was performed, three SVM classification models were established, based on the isotope data (Matrix A), the metallic element data (Matrix B), and the union data (Matrix C). For all of the 132 samples (99 WRT samples and 33 NWRT samples), 88 of them (including WRT and NWRT samples) were selected as a training class at random and the other 44 were put into a prediction class.

To estimate the performance of the SVM model, the sensitivity and specificity were calculated as follows:

Sens. = TP/TP + FN,

Spec. TN/TN + FP, (3)

where TP and FN represent the number of true positives and false negatives, respectively, and TN and FP denote the number of true negatives and false positives, respectively.

2.6. Variable Ranking. For SVM models, it is obvious that each variable does not contribute equally to prediction accuracy. In addition, some useless information may even have negative influence on prediction, so the significance of each measured isotope and element was investigated in this paper. For each feature, the column was removed from the data matrix and a new SVM model based on the incomplete data was built. In this way, each feature of the isotope and element data was separately removed, and 18 models were built. Then, the models were compared, and if a model showed a lower accuracy, the missing feature was considered important in provenance discrimination. Using this method, every variable was ranked by its contribution.

3. Results and Discussion

[[delta].sup.13]C in plant samples is mainly affected by the metabolic pathway of plant photosynthesis, so [[delta].sup.13]C is significantly different between different plants, while [[delta].sup.15]N is mainly under the influence of such regional agricultural activities as fertilization [12]. [[delta].sup.2]H and [[delta].sup.18]O are affected by atmospheric water cycle. Obviously, they have dimensional and land effects through the meteorological cycle of evaporation, condensation, and precipitation. Decreasing temperatures causes a progressive heavy-isotope depletion of the precipitation when the water vapour from oceans in equatorial regions moves to higher latitudes and altitudes [11]. [[delta].sup.2]H and [[delta].sup.18]O in plants are affected by [[delta].sup.2]H and [[delta].sup.18]O in the surrounding environment, so they are well used to characterize the origin of agricultural products [12]. The characteristics of metallic elements in plants are not only related to the composition of mineral elements in the soil, but also affected by varieties, climate, and agricultural activities [30,31]. Alkaline metals, especially Cs and Rb being easily mobilised in the soil and transported into plants, are good indicators of geographical identity [12]. Cu, Zn, and Cd in the soil will be affected by agricultural activities (organic fertilizer), so these elements in plants will also be affected [32-34]. In conclusion, it is necessary to evaluate the influence of these trace elements in the identification of tea samples.

The results of isotope and metallic element profiling are presented in Tables 4 and 5. The tables demonstrate a considerable difference in isotope and metallic content, but distinguishing the provenance of the sample just by these values proved difficult. Chemometric models are powerful tools in such situations.

Three SVM models based on isotope data, element content data, and coupled data were established, and their prediction results are shown in Table 6. The accuracy of isotope-SVM model reached 0.9318, and only 3 samples were mispredicted. Although the performance of element-SVM models (0.7727) was not very satisfactory, when it was applied coupled with isotope data, the model can greatly improve predictions, and the accuracy of the coupled model reached 0.9773.

The rank of each feature is reported in Table 7. In the table, [[delta].sup.2]H, [[delta].sup.18]O, Cs, Cu, Ca, and Rb contents are ranked highest, so they are the most significant indication in provenance analysis of WRT. Moreover, in further analysis, each feature was accumulatively assembled by its rank order and the variation of accuracy was plotted in Figure 2 when a new variable was added. As shown in Figure 2, the prediction accuracy was reduced as the [[delta].sup.18]O feature was added. This result may be caused by the overlapping chemical information between [[delta].sup.2]H and [[delta].sup.18]O because they have analogous variation in the hydrosphere. Therefore, the relationship between [[delta].sup.2]H and [[delta].sup.18]O was examined, and the correlation coefficient between [[delta].sup.2]H and [[delta].sup.18]O reached 0.8634, which strongly supports this assumption. Afterwards, the model achieved better performance as the variables of [[delta].sup.15]N and [[delta].sup.13]C were added. In element profiling, the SVM model achieved the best prediction result when only Cs and Cu contents were applied. The prediction accuracy decreased as more element features were used, so some metal elements were not significant in the identification of Wuyi-Rock tea geographical origin.

4. Conclusion

The origin place of Wuyi-Rock tea is the typical Danxia landform constituting purple soil, red soil, and moist sandy soil and the microenvironment is unique and exclusive. According to the results of tea identification, the importance of [[delta].sup.2]H and [[delta].sup.18]8O exactly reflects the particular climatic environment in Wuyishan. [[delta].sup.13]C mainly reflects the difference between different plants, and [[delta].sup.15]N is easily influenced by the agricultural activities [12]; therefore, the importance of [[delta].sup.13]C and [[delta].sup.15]N is not significant in the identification. Cs, Rb, and Sr have a higher contribution to Wuyi-Rock tea discrimination than great majority of elements, which has illustrated that the special geology in Wuyishan area provides unique features of trace elements for Wuyi-Rock tea. The contents of Cu, Ca, and Zn in Wuyi-Rock tea were affected by many factors, such as the kind of soil, fertilization, and tea varieties, so it needed further investigation and analysis of the relationship between identification and those aspects mentioned above in Wuyishan tea field.

In this paper, isotope and metallic element analyses demonstrate the potential for geographical origin discrimination of Wuyi-Rock tea. As a nonlinear model, SVM was performed for classification, and the chemical information of isotopes and metallic elements is complementary in provenance discrimination. In addition, the ranks of isotope and element features were carried out using established methods. The result shows that [[delta].sup.2]H, [[delta].sup.18]O, Cs, Cu, Ca, and Rb contents are significant in provenance analysis, [[delta].sup.2]H and [[delta].sup.18]O are interrelated, and not every element is helpful in geographical origin discrimination.

https://doi.org/10.1155/2017/5454231

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors' Contributions

Yun-xiao Lou and Xian-shu Fu equally contributed to this study.

Acknowledgments

The authors are grateful for the financial support from the Public Innovation Platform Analysis and Test Program of Zhejiang Province (no. 2016C37031) and the financial support from the National Public Welfare Industry Projects of China (no. 2012104019).

References

[1] A. Moghbel, H. Farjzadeh, N. Aghel, H. Agheli, and N. Raisi, "The effect of green tea on prevention of mouth bacterial infection, halitosis, and plaque formation on teeth," Iranian Journal of Toxicology, vol. 5, no. 14, pp. 502-515, 2011.

[2] L.-K. Han, T. Takaku, J. Li, Y. Kimura, and H. Okuda, "Antiobesity action of oolong tea," International Journal of Obesity, vol. 23, no. 1, pp. 98-105, 1999.

[3] K. Wang, F. Liu, Z. Liu et al., "Analysis of chemical components in oolong tea in relation to perceived quality," International Journal of Food Science and Technology, vol. 45, no. 5, pp. 913-920, 2010.

[4] X. Su, J. Duan, Y. Jiang, X. Duan, and F. Chen, "Polyphenolic profile and antioxidant activities of oolong tea infusion under various steeping conditions," International Journal of Molecular Sciences, vol. 8, no. 12, pp. 1196-1205, 2007.

[5] W.-Y. Han, F.-J. Zhao, Y.-Z. Shi, L.-F. Ma, and J.-Y. Ruan, "Scale and causes of lead contamination in Chinese tea," Environmental Pollution, vol. 139, no. 1, pp. 125-132, 2006.

[6] D. Lin, L. Zhu, W. He, and Y. Tu, "Tea plant uptake and translocation of polycyclic aromatic hydrocarbons from water and around air," Journal of Agricultural and Food Chemistry, vol. 54, no. 10, pp. 3658-3662, 2006.

[7] H.-P. Lv, Y.-J. Zhang, Z. Lin, and Y.-R. Liang, "Processing and chemical constituents of Pu-erh tea: a review," Food Research International, vol. 53, no. 2, pp. 608-618, 2013.

[8] Y. Uchiyama, M. Yamashita, M. Kato, T. Suzuki, M. Omori, and R. Chen, "Evaluation of the taste of tea with different degrees of fermentation using a taste sensing system," Sensors and Materials, vol. 23, no. 8, pp. 501-506, 2011.

[9] M. Chudzinska and D. Baralkiewicz, "Application of ICP-MS method of determination of 15 elements in honey with chemometric approach for the verification of their authenticity," Food and Chemical Toxicology, vol. 49, no. 11, pp. 2741-2749, 2011.

[10] Y. Chen, M.-Y. Xie, Y. Yan et al., "Discrimination of Ganoderma lucidum according to geographical origin with near infrared diffuse reflectance spectroscopy and pattern recognition techniques," Analytica Chimica Acta, vol. 618, no. 2, pp. 121-130, 2008.

[11] W. Dansgaard, "Stable isotopes in precipitation," Tellus, vol. 16, no. 4, pp. 436-468, 1964.

[12] S. Kelly, K. Heaton, and J. Hoogewerff, "Tracing the geographical origin of food: the application of multi-element and multi-isotope analysis," Trends in Food Science and Technology, vol. 16, no. 12, pp. 555-567, 2005.

[13] R. Nakashita, Y. Suzuki, F. Akamatsu, Y. Iizumi, T. Korenaga, and Y. Chikaraishi, "Stable carbon, nitrogen, and oxygen isotope analysis as a potential tool for verifying geographical origin of beef," Analytica Chimica Acta, vol. 617, no. 1-2, pp. 148-152, 2008.

[14] R. G. Crittenden, A. S. Andrew, M. LeFournour, M. D. Young, H. Middleton, and R. Stockmann, "Determining the geographic origin of milk in Australasia using multi-element stable isotope ratio analysis," International Dairy Journal, vol. 17, no. 5, pp. 421-428, 2007.

[15] D. Luo, H. Dong, H. Luo, Y. Xian, X. Guo, and Y. Wu, "Multi-element (C, N, H, O) stable isotope ratio analysis for determining the geographical origin of pure milk from different regions," Food Analytical Methods, vol. 9, no. 2, pp. 437-442, 2016.

[16] F. Camin, N. Dordevic, R. Wehrens et al., "Climatic and geographical dependence of the H, C and O stable isotope ratios of Italian wine," Analytica Chimica Acta, vol. 853, no. 1, pp. 384-390, 2015.

[17] A. A. D'Archivio, A. Giannitto, A. Incani, and S. Nisi, "Analysis of the mineral composition of Italian saffron by ICP-MS and classification of geographical origin," Food Chemistry, vol. 157, no. 4, pp. 485-489, 2014.

[18] P. P. Coetzee, F. P. Van Jaarsveld, and F. Vanhaecke, "Intraregional classification of wine via ICP-MS elemental fingerprinting," Food Chemistry, vol. 164, no. 3, pp. 485-492, 2014.

[19] V. S. Selih, M. Sala, and V. Drgan, "Multi-element analysis of wines by ICP-MS and ICP-OES and their classification according to geographical origin in Slovenia," Food Chemistry, vol. 153, pp. 414-423, 2014.

[20] M. Madejczyk and D. Baralkiewicz, "Characterization of Polish rape and honeydew honey according to their mineral contents using ICP-MS and F-AAS/AES," Analytica Chimica Acta, vol. 617, no. 1-2, pp. 11-17, 2008.

[21] T. S. Pilgrim, R. J. Watling, and K. Grice, "Application of trace element and stable isotope signatures to determine the provenance of tea (Camelliasinensis) samples," Food Chemistry, vol. 118, no. 4, pp. 921-926, 2010.

[22] M. J. Martelo-Vidal and M. Vazquez, "Rapid authentication of white wines. Part 1: classification by designation of origin," Agro Food Industry Hi-Tech, vol. 25, no. 5, pp. 44-48, 2014.

[23] E. Furia, A. Naccarato, G. Sindona, G. Stabile, and A. Tagarelli, "Multielement fingerprinting as a tool in origin authentication of pgi food products: tropea red onion," Journal of Agricultural and Food Chemistry, vol. 59, no. 15, pp. 8450-8457, 2011.

[24] M. J. Martelo-Vidal and M. Vazquez, "Rapid authentication of white wines part 2: classification by grape variety," Agro Food Industry Hi Tech, vol. 25, no. 6, p. 6, 2014.

[25] J. D. Giraldi and I. B. Lopes, "The country of origin effect for Brazilian fresh fruits: a study using Partial Least Squares Procedures," Journal of Management Research, vol. 4, no. 4, 19 pages, 2012.

[26] R. M. Barbosa, B. L. Batista, R. M. Varrique, V. A. Coelho, A. D. Campiglia, and F. Barbosa, "The use of advanced chemometric techniques and trace element levels for controlling the authenticity of organic coffee," Food Research International, vol. 61, pp. 246-251, 2014.

[27] B. L. Batista, L. R. S. da Silva, B. A. Rocha et al., "Multi-element determination in Brazilian honey samples by inductively coupled plasma mass spectrometry and estimation of geographic origin with data mining techniques," Food Research International, vol. 49, no. 1, pp. 209-215, 2012.

[28] Q. Chen, J. Zhao, C. H. Fang, and D. Wang, "Feasibility study on identification of green, black and Oolong teas using near-infrared reflectance spectroscopy based on support vector machine (SVM)," Spectrochimica Acta, vol. 66, no. 3, pp. 568-574, 2007.

[29] O. Ivanciuc, "Applications of support vector machines in chemistry," Reviews in Computational Chemistry, vol. 23, pp. 291-400, 2007.

[30] J. S. Chen, B. S. Deng, M. Pan et al., "Geographical tendencies of trace element contents in soils derived from granite, basalt and limestone of eastern china," Pedosphere, vol. 3, no. 1, pp. 45-55, 1993.

[31] H. Zhao, B. Guo, Y. Wei et al., "Determining the geographic origin of wheat using multielement analysis and multivariate statistics," Journal of Agricultural and Food Chemistry, vol. 59, no. 9, pp. 4397-4402, 2011.

[32] P. J. Florijn, J. A. Nelemans, and M. L. van Beusichem, "The influence of the form of nitrogen nutrition on uptake and distribution of cadmium in lettuce varieties," Journal of Plant Nutrition, vol. 15, no. 11, pp. 2405-2416, 1992.

[33] P. Gramatica, F. Battaini, E. Giani, E. Papa, R. J. A. Jones, and R. M. Cenci, "Multivariate analysis of heavy metal concentrations in soils and mosses of two North-Italy regions," Fresenius Environmental Bulletin, vol. 15, no. 8 A, pp. 731-737, 2006.

[34] K. W. Perrott, B. F. L. Smith, and R. H. E. Inkson, "The reaction of fluoride with soils and soil minerals," Journal of Soil Science, vol. 27, no. 1, pp. 58-67, 1976.

Yun-xiao Lou, Xian-shu Fu, Xiao-ping Yu, Zi-hong Ye, Hai-feng Cui, and Ya-fen Zhang

Zhejiang Provincial Key Laboratory of Biometrology and Inspection & Quarantine, College of Life Sciences, China Jiliang University, Hangzhou 310018, China

Correspondence should be addressed to Zi-hong Ye; zhye@cjlu.edu.cn

Received 18 January 2017; Revised 10 March 2017; Accepted 23 March 2017; Published 3 April 2017

Academic Editor: Constantinos Georgiou

Caption: Figure 1: Tea samples collected in Wuyishan.

Caption: Figure 2: (a) Feature selection of isotope and prediction accuracy. (b) Feature selection of metallic element and prediction accuracy.
Table 1: Detailed information of the samples.

                                                    Size of
Number                Producing area               sample (a)

A01           Shangmei, Wuyi, Fujian Province          9
A02           Xingcun, Wuyi, Fujian Province           9
A03             Wufu, Wuyi, Fujian Province            9
A04            Langu, Wuyi, Fujian Province            9
A05        Chongan Street, Wuyi, Fujian Province       9
A06        Xinfeng Street, Wuyi, Fujian Province       9
A07          Yangzhuang, Wuyi, Fujian Province         9
A08           Xingtian, Wuyi, Fujian Province          9
A09            Xiamei, Wuyi, Fujian Province           9
A10            Wutun, Wuyi, Fujian Province            9
A11         Wuyi Street, Wuyi, Fujian Province         9
N01              Jianyang, Fujian Province             3
N02               Jianou, Fujian Province              3
N03            Zhangzhou A, Fujian Province            3
N04            Zhangzhou B, Fujian Province            3
N05              Quanzhou, Fujian Province             3
N06                  Guangxi Province                  3
N07                  Guizhou Province                  3
N08              Ganzhou, Jiangxi Province             3
N09              Wuyuan, Jiangxi Province              3
N10               Songxi, Fujian Province              3
N11              Zhenghe, Fujian Province              3

Number   Type (b)

A01         A
A02         A
A03         A
A04         A
A05         A
A06         A
A07         A
A08         A
A09         A
A10         A
A11         A
N01         N
N02         N
N03         N
N04         N
N05         N
N06         N
N07         N
N08         N
N09         N
N10         N
N11         N

(a) The number of the tea samples from the production site.

(b) A: Wuyi-Rock tea; N: non-Wuyi-Rock tea.

Table 2: Optimized instrumental conditions and parameters
of the ICP-MS.

Parameter                 Setting value

Forward power                1200 W
Scanning times              100 times
Scanning mode              Peak height
Dwell time                    10 ms
Acquisition time              20 s
Sample uptake rate           1mL/min
Coolant gas flow rate       14 L/min
Auxiliary gas flow rate    0.75 L/min
Nebulizer gas flow rate    0.92 L/min

Table 3: Main measurement parameters for Ca, Mg, and Mn.

                                   Ethyne
                                    gas         Air
                                    flow        flow
Element   Wavelength   Bandpass     rate        rate

Ca         422.7 nm     2.6 nm    2.6 L/min   9.4 L/min
Mg         285.2 nm     2.6 nm     2 L/min    9.4 L/min
Mn         279.5 nm     2.3 nm    2.3 L/min   9.4 L/min

Table 4: Results of provenance experiment: averages values of
4 stable isotopes in WRT and NWRT samples.

Number     [[delta].sup.13]C        [[delta].sup.2]H

A01       -26.77 [+ or -] 0.68    -80.90 [+ or -] 9.90
A02       -27.90 [+ or -] 0.76    -74.85 [+ or -] 7.54
A03       -27.14 [+ or -] 0.50    -77.05 [+ or -] 8.27
A04       -27.78 [+ or -] 0.96    -86.32 [+ or -] 0.83
A05       -26.79 [+ or -] 0.67    -78.34 [+ or -] 4.39
A06       -27.95 [+ or -] 0.78    -72.42 [+ or -] 5.61
A07       -27.56 [+ or -] 1.61    -86.27 [+ or -] 12.88
A08       -26.97 [+ or -] 0.56    -78.89 [+ or -] 4.15
A09       -27.50 [+ or -] 0.33    -86.99 [+ or -] 14.94
A10       -27.18 [+ or -] 0.12    -89.79 [+ or -] 6.63
A11       -26.90 [+ or -] 0.21    -73.76 [+ or -] 2.04
N01       -28.56 [+ or -] 0.19   -102.59 [+ or -] 0.89
N02       -26.96 [+ or -] 0.08    -89.39 [+ or -] 2.34
N03       -27.77 [+ or -] 0.07   -111.98 [+ or -] 1.55
N04       -28.34 [+ or -] 0.08   -111.98 [+ or -] 0.78
N05       -27.29 [+ or -] 0.15    -98.31 [+ or -] 1.83
N06       -27.84 [+ or -] 0.06   -107.15 [+ or -] 1.60
N07       -26.80 [+ or -] 0.09   -112.65 [+ or -] 0.64
N08       -28.04 [+ or -] 0.16   -103.89 [+ or -] 2.10
N09       -28.35 [+ or -] 0.11   -107.52 [+ or -] 1.83
N10       -27.88 [+ or -] 0.13   -107.75 [+ or -] 1.03
N11       -26.85 [+ or -] 0.14   -102.27 [+ or -] 1.97

Number    [[delta].sup.18]O     [[delta].sup.15]N

A01      24.44 [+ or -] 0.33   1.14 [+ or -] 1.17
A02      23.39 [+ or -] 2.49   2.26 [+ or -] 0.93
A03      23.01 [+ or -] 0.81   2.55 [+ or -] 1.54
A04      22.18 [+ or -] 1.11   1.71 [+ or -] 1.06
A05      23.96 [+ or -] 1.49   2.76 [+ or -] 0.88
A06      25.05 [+ or -] 0.37   1.06 [+ or -] 2.46
A07      21.65 [+ or -] 2.63   1.40 [+ or -] 1.38
A08      24.24 [+ or -] 0.72   4.74 [+ or -] 2.10
A09      23.36 [+ or -] 3.10   2.54 [+ or -] 1.49
A10      24.56 [+ or -] 1.57   3.71 [+ or -] 0.96
A11      24.41 [+ or -] 0.81   1.83 [+ or -] 0.52
N01      18.63 [+ or -] 0.17   6.12 [+ or -] 0.22
N02      20.41 [+ or -] 0.44   3.63 [+ or -] 0.03
N03      18.34 [+ or -] 0.16   7.06 [+ or -] 0.22
N04      18.24 [+ or -] 0.22   6.78 [+ or -] 0.10
N05      22.39 [+ or -] 0.39   4.19 [+ or -] 0.07
N06      18.71 [+ or -] 0.25   4.14 [+ or -] 0.02
N07      16.65 [+ or -] 0.30   6.50 [+ or -] 0.07
N08      23.24 [+ or -] 0.20   1.98 [+ or -] 0.24
N09      18.48 [+ or -] 0.02   0.79 [+ or -] 0.02
N10      19.11 [+ or -] 0.35   4.91 [+ or -] 0.24
N11      21.09 [+ or -] 0.59   1.30 [+ or -] 0.18

Table 5: Results of provenance experiment: averages values of
14 element concentrations ([micro]g/g) in WRT and NWRT samples.

Number               Ti                   Cr

A01         27.56 [+ or -] 1.25   5.31 [+ or -] 1.05
A02         27.06 [+ or -] 1.07   5.43 [+ or -] 0.26
A03         25.78 [+ or -] 1.39   5.04 [+ or -] 0.65
A04         23.72 [+ or -] 2.74   5.54 [+ or -] 0.16
A05         27.06 [+ or -] 3.82   5.29 [+ or -] 0.71
A06         28.22 [+ or -] 5.55   4.99 [+ or -] 0.18
A07         24.39 [+ or -] 2.75   5.01 [+ or -] 0.59
A08         24.67 [+ or -] 5.80   4.62 [+ or -] 0.29
A09         26.22 [+ or -] 3.91   4.07 [+ or -] 0.95
A10         26.94 [+ or -] 0.67   4.89 [+ or -] 0.27
A11         32.06 [+ or -] 2.84   4.68 [+ or -] 0.40
N01         20.33 [+ or -] 0.50   3.88 [+ or -] 0.21
N02         28.94 [+ or -] 0.25   5.17 [+ or -] 0.08
N03         24.06 [+ or -] 0.48   4.89 [+ or -] 0.09
N04         20.78 [+ or -] 0.42   4.63 [+ or -] 0.12
N05         22.72 [+ or -] 0.51   4.94 [+ or -] 0.09
N06         31.17 [+ or -] 2.90   4.50 [+ or -] 0.07
N07         39.61 [+ or -] 1.70   4.26 [+ or -] 0.02
N08         21.17 [+ or -] 0.44   4.74 [+ or -] 0.08
N09         29.72 [+ or -] 0.92   4.86 [+ or -] 0.06
N10         21.06 [+ or -] 0.86   4.62 [+ or -] 0.08
N11         29.17 [+ or -] 1.17   4.28 [+ or -] 0.09
Ref                27.50                 4.42
LOD (ng/L)         0.0019                0.025

Number               Co                     Ni

A01          0.88 [+ or -] 0.28    12.58 [+ or -] 2.61
A02          0.26 [+ or -] 0.17    10.19 [+ or -] 1.22
A03          0.14 [+ or -] 0.02     9.68 [+ or -] 1.28
A04          0.27 [+ or -] 0.08    12.61 [+ or -] 1.17
A05          0.24 [+ or -] 0.02    10.58 [+ or -] 0.61
A06          0.48 [+ or -] 0.15    10.09 [+ or -] 0.83
A07          0.33 [+ or -] 0.08    11.63 [+ or -] 2.02
A08          0.41 [+ or -] 0.36    10.68 [+ or -] 2.50
A09          0.30 [+ or -] 0.18     9.04 [+ or -] 2.24
A10          0.88 [+ or -] 0.17    11.56 [+ or -] 1.57
A11          0.28 [+ or -] 0.04    10.52 [+ or -] 0.20
N01          0.17 [+ or -] 0.003    9.16 [+ or -] 0.28
N02          0.93 [+ or -] 0.01    14.23 [+ or -] 0.34
N03          0.19 [+ or -] 0.02    10.79 [+ or -] 0.24
N04          0.22 [+ or -] 0        9.73 [+ or -] 0.32
N05          0.24 [+ or -] 0.01    10.08 [+ or -] 0.48
N06          0.64 [+ or -] 0.02    11.75 [+ or -] 0.17
N07          1.35 [+ or -] 0.03    17.44 [+ or -] 0.25
N08          0.48 [+ or -] 0.02    11.89 [+ or -] 0.79
N09          0.24 [+ or -] 0.01    16.34 [+ or -] 0.18
N10          0.28 [+ or -] 0.01     9.59 [+ or -] 0.29
N11          0.34 [+ or -] 0.02    10.31 [+ or -] 0.26
Ref                 0.42                  13.62
LOD (ng/L)         0.001                  0.008

Number               Cu                   Zn

A01         11.65 [+ or -] 0.95   30.56 [+ or -] 1.17
A02         11.40 [+ or -] 0.12   29.61 [+ or -] 3.51
A03         10.08 [+ or -] 0.72   27.67 [+ or -] 3.66
A04         15.78 [+ or -] 1.77   29.67 [+ or -] 3.18
A05         10.54 [+ or -] 1.89   33.11 [+ or -] 9.73
A06         11.20 [+ or -] 1.60   32.28 [+ or -] 2.60
A07         12.92 [+ or -] 1.22   27.83 [+ or -] 7.55
A08         12.10 [+ or -] 1.61   28.11 [+ or -] 9.48
A09         11.38 [+ or -] 0.66   30.94 [+ or -] 6.46
A10         13.68 [+ or -] 0.72   21.06 [+ or -] 2.51
A11         12.06 [+ or -] 0.69   52.17 [+ or -] 27.66
N01         10.39 [+ or -] 0.65   19.17 [+ or -] 0.29
N02         13.39 [+ or -] 0.31   40.28 [+ or -] 1.02
N03         11.67 [+ or -] 0.14   18.00 [+ or -] 0.58
N04         11.69 [+ or -] 0.07   18.33 [+ or -] 1.01
N05         11.86 [+ or -] 0.67   19.39 [+ or -] 0.98
N06         15.68 [+ or -] 0.18   30.33 [+ or -] 2.17
N07         23.67 [+ or -] 0.17   36.72 [+ or -] 2.12
N08         30.39 [+ or -] 0.54   35.33 [+ or -] 8.99
N09         27.89 [+ or -] 0.42   33.06 [+ or -] 0.35
N10         13.36 [+ or -] 0.32   18.56 [+ or -] 1.34
N11         15.85 [+ or -] 0.37   33.50 [+ or -] 0.44
Ref                34.00                 30.17
LOD (ng/L)         0.045                 0.091

Number               Rb                   Cd

A01         50.89 [+ or -] 10.16   0.07 [+ or -] 0.01
A02         40.11 [+ or -] 6.92    0.07 [+ or -] 0.01
A03         46.22 [+ or -] 15.16   0.04 [+ or -] 0.01
A04         66.22 [+ or -] 22.27   0.07 [+ or -] 0.01
A05         35.72 [+ or -] 15.54   0.08 [+ or -] 0.06
A06         59.39 [+ or -] 18.63   0.06 [+ or -] 0.02
A07         69.08 [+ or -] 17.80   0.10 [+ or -] 0.07
A08         80.78 [+ or -] 35.31   0.06 [+ or -] 0.05
A09         40.11 [+ or -] 19.16   0.09 [+ or -] 0.04
A10         54.94 [+ or -] 6.42    0.04 [+ or -] 0.02
A11         53.11 [+ or -] 17.23   0.04 [+ or -] 0.01
N01         43.17 [+ or -] 0.93    0.03 [+ or -] 0.01
N02        104.33 [+ or -] 0.67    0.09 [+ or -] 0.01
N03         73.83 [+ or -] 0.76    0.04 [+ or -] 0.01
N04         69.22 [+ or -] 0.10    0.04 [+ or -] 0.01
N05         45.11 [+ or -] 0.79    0.04 [+ or -] 0.01
N06        347.78 [+ or -] 5.85    0.06 [+ or -] 0.02
N07         59.17 [+ or -] 0.88    0.07 [+ or -] 0.01
N08        117.17 [+ or -] 1.36    0.36 [+ or -] 0.16
N09         76.56 [+ or -] 0.69    0.10 [+ or -] 0.01
N10         78.83 [+ or -] 1.32    0.04 [+ or -] 0.01
N11         92.56 [+ or -] 2.18    0.08 [+ or -] 0.02
Ref                60.50                  0.10
LOD (ng/L)         0.001                  0.002

Number               Cs                   Ba

A01          0.53 [+ or -] 0.26   14.24 [+ or -] 1.53
A02          0.48 [+ or -] 0.14   10.64 [+ or -] 2.63
A03          0.44 [+ or -] 0.18   12.09 [+ or -] 4.20
A04          0.32 [+ or -] 0.10   32.78 [+ or -] 13.76
A05          0.26 [+ or -] 0.19   15.01 [+ or -] 13.13
A06          0.34 [+ or -] 0.03    7.68 [+ or -] 0.34
A07          1.82 [+ or -] 1.22   22.17 [+ or -] 2.89
A08          1.26 [+ or -] 0.09   13.38 [+ or -] 3.94
A09          0.57 [+ or -] 0.28    8.10 [+ or -] 0.98
A10          0.86 [+ or -] 0.48   12.04 [+ or -] 2.37
A11          0.39 [+ or -] 0.21   12.87 [+ or -] 7.39
N01          1.08 [+ or -] 0.04   18.78 [+ or -] 0.79
N02         2.406 [+ or -] 0.03   12.46 [+ or -] 0.22
N03          1.58 [+ or -] 0.02   16.83 [+ or -] 0.17
N04          1.73 [+ or -] 0.01   22.11 [+ or -] 0.92
N05          1.64 [+ or -] 0.02    8.04 [+ or -] 0.12
N06          2.80 [+ or -] 0.07   12.59 [+ or -] 0.38
N07          1.20 [+ or -] 0.02   13.13 [+ or -] 0.21
N08          3.93 [+ or -] 0.03   31.83 [+ or -] 0.60
N09          0.83 [+ or -] 0.02   41.56 [+ or -] 0.54
N10          2.12 [+ or -] 0.02   24.50 [+ or -] 0.44
N11          1.04 [+ or -] 0.02   26.50 [+ or -] 0.67
Ref                 0.97                 47.33
LOD (ng/L)         0.001                 0.007

Number               Sr                    Ca

A01         15.56 [+ or -] 1.92   3572.22 [+ or -] 233.53
A02         15.00 [+ or -] 2.89   4122.22 [+ or -] 296.43
A03         15.56 [+ or -] 1.92   4288.89 [+ or -] 211.04
A04         13.33 [+ or -] 2.89   4027.78 [+ or -] 511.08
A05          9.44 [+ or -] 0.96   3816.67 [+ or -] 464.58
A06         13.33 [+ or -] 3.33   3711.11 [+ or -] 475.61
A07         13.89 [+ or -] 2.54   4027.78 [+ or -] 500.09
A08         12.78 [+ or -] 3.47   3872.22 [+ or -] 453.79
A09         12.22 [+ or -] 3.85   3661.11 [+ or -] 431.51
A10         12.78 [+ or -] 3.47   3466.67 [+ or -] 180.28
A11         12.78 [+ or -] 3.47   3200.00 [+ or -] 251.66
N01         15.00 [+ or -] 2.89   3666.67 [+ or -] 196.50
N02         12.22 [+ or -] 3.85   3605.56 [+ or -] 9.62
N03         15.00 [+ or -] 2.89   3450.00 [+ or -] 196.50
N04         13.89 [+ or -] 2.54   3600.00 [+ or -] 50.00
N05         14.44 [+ or -] 3.85   3205.56 [+ or -] 110.97
N06         15.56 [+ or -] 1.92   2688.89 [+ or -] 25.46
N07         10.00 [+ or -] 5.77   3733.33 [+ or -] 72.65
N08         15.00 [+ or -] 2.89   3866.67 [+ or -] 183.33
N09         16.67 [+ or -] 0      3405.56 [+ or -] 161.88
N10         15.56 [+ or -] 1.92   4283.33 [+ or -] 120.18
N11         15.00 [+ or -] 2.89   3783.33 [+ or -] 44.10
Ref                33.33                 12633.33
LOD (ng/L)         0.003                   0.021

Number                 Mg                        Mn

A01          2244.44 [+ or -] 133.68   1500.00 [+ or -] 317.54
A02          2266.67 [+ or -] 109.29    888.89 [+ or -] 350.13
A03          2305.56 [+ or -] 250.19    350.00 [+ or -] 152.75
A04          2633.33 [+ or -] 217.94   1038.89 [+ or -] 330.96
A05          2333.33 [+ or -] 389.80    966.67 [+ or -] 327.87
A06          2161.11 [+ or -] 110.97    872.22 [+ or -] 227.51
A07          2277.78 [+ or -] 409.72   1088.89 [+ or -] 460.17
A08          2127.78 [+ or -] 437.90   1194.44 [+ or -] 395.23
A09          2072.22 [+ or -] 58.53     905.56 [+ or -] 194.60
A10          2427.78 [+ or -] 122.85   1588.89 [+ or -] 242.86
A11          2327.78 [+ or -] 48.11     694.44 [+ or -] 327.59
N01          2227.78 [+ or -] 25.46     827.78 [+ or -] 82.21
N02          2294.44 [+ or -] 9.62     1394.44 [+ or -] 25.46
N03          2022.22 [+ or -] 75.15     650.00 [+ or -] 33.33
N04          1950.00 [+ or -] 16.67    1161.11 [+ or -] 34.69
N05          1972.22 [+ or -] 85.53     727.78 [+ or -] 19.24
N06          2183.33 [+ or -] 28.87     788.89 [+ or -] 9.62
N07          2822.22 [+ or -] 41.94    1705.56 [+ or -] 25.46
N08          2216.67 [+ or -] 92.80    1483.33 [+ or -] 16.67
N09          2566.67 [+ or -] 57.74     872.22 [+ or -] 19.24
N10          2183.33 [+ or -] 44.10    1244.44 [+ or -] 9.62
N11          2400.00 [+ or -] 100.00    733.33 [+ or -] 16.67
Ref                  2800.00                   1333.33
LOD (ng/L)            0.005                     0.039

Table 6: Predicting results obtained by SVM.

Data type                Sensitivity      Specificity

Isotope data            0.9394 (31/33)   0.9091 (10/11)
Metallic element data   0.8788 (29/33)   0.4545 (5/11)
Coupled data            0.9697 (32/33)     1 (11/11)

Data type                  Accuracy

Isotope data            0.9318 (41/44)
Metallic element data   0.7727 (34/44)
Coupled data            0.9773 (43/44)

Table 7: Rank results of isotope and metallic element data.

       Isotopes

Ranking       Variable

1         [[delta].sup.2]H
2         [[delta].sup.18]O
3         [[delta].sup.15]N
4         [[delta].sup.13]C

 Metallic element

Ranking   Variable

1            Cs
2            Cu
3            Ca
4            Rb
5            Sr
6            Ba
7            Cd
8            Cr
9            Ni
10           Zn
11           Ti
12           Mg
13           Mn
14           Co
COPYRIGHT 2017 Hindawi Limited
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2017 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Lou, Yun-xiao; Fu, Xian-shu; Yu, Xiao-ping; Ye, Zi-hong; Cui, Hai-feng; Zhang, Ya-fen
Publication:Journal of Analytical Methods in Chemistry
Article Type:Report
Date:Jan 1, 2017
Words:7081
Previous Article:Development of a LC-MS/MS Method for the Simultaneous Detection of Tricarboxylic Acid Cycle Intermediates in a Range of Biological Matrices.
Next Article:Electrochemical Evaluation of trans-Resveratrol Levels in Red Wine Based on the Interaction between Resveratrol and Graphene.
Topics:

Terms of use | Privacy policy | Copyright © 2022 Farlex, Inc. | Feedback | For webmasters |