Next Article in Journal
Photodynamic Inactivation of Foodborne Bacteria: Screening of 32 Potential Photosensitizers
Next Article in Special Issue
Chemical Characterization and Temporal Variability of Pasta Condiment By-Products for Sustainable Waste Management
Previous Article in Journal
Food Waste from Campus Dining Hall as a Potential Feedstock for 2,3-Butanediol Production via Non-Sterilized Fermentation
Previous Article in Special Issue
Varietal Authenticity Assessment of QTMJ Tea Using Non-Targeted Metabolomics and Multi-Elemental Analysis with Chemometrics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Potential Use of Near Infrared Spectroscopy (NIRS) to Determine the Heavy Metals and the Percentage of Blends in Tea

by
Isabel Revilla
1,
Miriam Hernández Jiménez
1,
Iván Martínez-Martín
1,
Patricia Valderrama
2,
Marta Rodríguez-Fernández
1 and
Ana M. Vivar-Quintana
1,*
1
Food Technology, Universidad de Salamanca, E.P.S. de Zamora, Avenida Requejo 33, 49022 Zamora, Spain
2
Department of Chemistry, Universidade Tecnológica Federal do Paraná (UTFPR), Via Rosalina Maria dos Santos 1233, Campo Mourão 87301-899, Paraná, Brazil
*
Author to whom correspondence should be addressed.
Foods 2024, 13(3), 450; https://doi.org/10.3390/foods13030450
Submission received: 21 December 2023 / Revised: 26 January 2024 / Accepted: 27 January 2024 / Published: 31 January 2024

Abstract

:
The following study analyzed the potential of Near Infrared Spectroscopy (NIRS) to predict the metal composition (Al, Pb, As, Hg and Cu) of tea and for establishing discriminant models for pure teas (green, red, and black) and their different blends. A total of 322 samples of pure black, red, and green teas and binary blends were analyzed. The results showed that pure red teas had the highest content of As and Pb, green teas were the only ones containing Hg, and black teas showed higher levels of Cu. NIRS allowed to predict the content of Al, Pb, As, Hg, and Cu with ratio performance deviation values > 3 for all of them. Additionally, it was possible to discriminate pure samples from their respective blends with an accuracy of 98.3% in calibration and 92.3% in validation. However, when the samples were discriminated according to the percentage of blending (>95%, 95–85%, 85–75%, or 75–50% of pure tea) 100% of the samples of 10 out of 12 groups were correctly classified in calibration, but only the groups with a level of pure tea of >95% showed 100% of the samples as being correctly classified as to validation.

1. Introduction

After water, tea is one of the most widely consumed non-alcoholic beverages in the world. Tea production worldwide had increased from 5.9 billion kg in 2018 to 6.8 billion kg in 2023 with a projected growth rate of 2.7%. The global tea market was valued at around 120.7 billion U.S. dollars in 2023 and is expected to grow annually by 5.3% until 2028 [1].
Tea is processed from the buds and new stems of Camellia sinensis which is native to southwest China [2]. The chemical composition of tea includes carbohydrates, amino acids, proteins, minerals, polyphenols, and alkaloids; it provides potential health benefits [3,4] and has important physiological properties [5]. Tea is considered to be antioxidant, anti-inflammatory, antiproliferative, antimutagenic, and antibacterial; it confers protection against cardiovascular diseases, hyperglycemia, metabolic disorders, and some cancers and may also reduce the risk of osteoporotic fractures in the elderly [6]. In addition to all of this, it contains caffeine, the consumption of which can provide a refreshing sensation, relieve pain and inflammation, encourage metabolism, and has diuretic effects [7].
Tea can be classified into six main categories based on their different processing techniques, which have been named green tea (unfermented), yellow and white teas (slightly fermented, 10–20%), Oolong tea (semi-fermented, 30–60%), black tea (fully fermented, 80–100%), and dark tea (post-fermented, ≤100%). In addition to these, Pu’er, also called red tea, is a unique microbial fermented tea produced from the sun-dried leaves of large-leaf tea species in the Yunnan province of China; it is attracting increasing interest which may be due to its multiple health benefits [8]. The composition, bioactive compound content, and antioxidant activity of tea are influenced by the different processing technologies applied, which means that different types of tea obtained from the same cultivar may vary greatly [9]. Green tea has thus been shown to provide protection against cancer by scavenging free radicals, dark tea may decrease the risk of cardiovascular disease, and black and oolong teas seem to possess outstanding antioxidant properties [10].
Based on their mineral composition, the leaves of Camellia sinensis are a source of elements such as zinc, manganese, iron, copper, magnesium, titanium, aluminum, strontium, bromine, sodium, potassium, phosphorus, iodine, and fluorine. Some of these minerals such as Fe, Mn, and Zn are components of important enzymes or are involved in a number of physiological processes, which means they are considered essential [11,12]. However, some elements are toxic to human health, such as As, Cr, Cd, Co, Ni, and Pb [13,14].
Moreover, different growing areas show variations in growing conditions, which include climate, rainfall, altitude, soil, fertilizers, microelements, and processing procedures [15], and these are all factors which contribute to the chemical composition of tea leaves [16]. The mineral composition of tea leaves appears to be influenced by the geographical origin, soil type, tea variety, climatic conditions, and other seasonal changes [17]. The presence of toxic heavy metals such as mercury, lead, and cadmium in tea leaves is probably related to soil contamination, atmospheric deposition, chemical and organic fertilizer applications, or irrigation with contaminated water and waste disposal [6]. High levels of heavy metals in the soil could inhibit the growth of tea trees by altering the biomass while reducing tea quality by changing the content of major ingredients [18] and therefore the quality of the tea. For instance, the application of phosphate fertilizers may be a major source of Cd in agricultural soils [19,20]. On the other hand, inorganic amendments were mainly responsible for the increase of Cr and Cu in tea growing soils while organic amendments were the influential factors in the enrichment of other metals [19].
As a result, the determination of mineral composition in tea leaves is relevant to industry and consumers. It is usually carried out by means of atomic absorption spectrometry (AAS) [12,21,22,23], which has low operational costs and provides good analytical performance although only a few elements can be determined simultaneously [24] by inductively coupled plasma mass spectroscopy (ICP-MS) or inductively coupled plasma atomic emission techniques (ICP-AES) which are accurate and expensive [25,26,27]. X-ray fluorescence (XRF) spectroscopy is an emerging screening technique for determining the mineral/elemental content in teas [4,11]. It is fast and accurate, non-destructive, and requires only minimal sample preparation. One of its limitations is the partial sensitivity of the method for high mass elements, which is due to high X-ray scattering (which can result in increased background intensities) [24].
In contrast, near infrared spectroscopy (NIRS) is a rapid, simple, multi-parametrical, and non-destructive technique. Moreover, NIRS shows excellent analytical performance in the prediction of the main contents of tea, including tea polyphenols, catechins, caffeine, theaflavins, chlorophyll, and sensory properties [7,28,29]. However, although NIR spectroscopy has been successfully used to predict the mineral composition in several matrices such as meat, cheese, pulses, quinoa, alfalfa, and propolis [30,31,32,33,34,35], as far as we know it has not been used for this purpose in tea.
Blending of teas with different characteristics is a common practice in the tea industry for several reasons, such as reducing costs by using lower grade teas [36], improving quality by obtaining a better taste [37] and fulfilling the consumer expectations by preparing teas with new characteristics or standardizing the sensory properties of a brand [36]. Moreover, blending different herbal teas of different origins and characteristics with tea to take advantage of the specific sensory properties and some medicinal values of this plants is increasingly consumed as wellness beverages [38,39]. In fact, a wide range of blends has been studied, from very low [40] to high blending percentages, up to 50% in some case. These blends can be made from different types of tea or by mixing tea with other plants for the preparation of other herbal teas [38,41,42]. In this scenario, product quality control and fraud prevention by operators is becoming increasingly important. In fact, fraud control is mandatory for industries certified according to GFSI (Global Food Safety Initiative) standards, such as IFS Food (International Featured Standards) and BRCGS Food (British Retail Consortium Global Standards). For example, mixing with leaves of lower quality species is considered adulteration [43]. Therefore, nowadays clean, low-cost and non-invasive techniques have been tested to assess the mixture, such as hyperspectral imaging [40], electronic nose [41], or visible and near-infrared spectroscopy [42,44].
Based on this background, the main purpose of this work was therefore to study the potential of near infrared spectroscopy for predicting the metal composition (Al, Pb, As, Hg and Cu) of green, red, and black tea samples and their mixes. In addition, the possibility of establishing classification models for the three types of tea (green, red, and black) and their blends using mineral composition and NIR spectra was also considered.

2. Materials and Methods

2.1. Samples

A total of 322 samples were analyzed. Samples of commercial black (7 samples), green (11 samples), and red (8 samples) teas were purchased in local markets (Zamora and Salamanca, Spain). The binary blends of two out three tea types were prepared in variable percentages ranging from 99 to 50% of tea of interest (1 to 50% of second tea) and thoroughly mixed. As the samples were purchased in supermarkets or small shops data on the geographical origin, the processing conditions, and the season they were picked were not available. The 16 pure samples, corresponding to 7 black, 11 green and 8 red tea samples, were analysed for their mineral content. In addition, the 65 samples that were used for the external calibration of the model were also analysed, these samples consisted of randomly selected mixed samples in different proportions. For the remaining samples (241) the metal contents were calculated based on the metal concentrations of each of the samples used to make each mixture
To ensure that the best mixing of the samples was obtained, a mixing protocol was established. The protocol conditions were as follows: the individual amounts of each of the teas used were weighed and placed in a hand-operated drum mixer with 500 g capacity containers for 15 s. Both the mixtures (241) for which the metal content was calculated and the 65 mixtures for which the metal content was established by laboratory analysis were prepared following this protocol.

2.2. Chemical Analyses

Element concentrations of Al, As, Cu, Hg, and Pb were determined by ICP-MS in pure teas (7 black, 11 green, and 8 red teas) and in the 65 samples of the external validation test. Approximately 0.2 g of each sample was placed in Teflon vessels with 5 mL of HNO3 and then placed in a Milestone microwave digestion system (Ethos Sel Milestone, Ontario, ON, Canada) and 1000 watts of power was applied for 5 min. The samples were allowed to cool and then 5 mL of HNO3 and 1 mL of H2O2 (30%) were added. The samples were then treated with 1000 watts for 10 min. When the samples were at room temperature, they were made up to 100 mL with distilled water and stored at 4 °C until analysis. The concentrations of the analyzed elements were quantified using an Agilent 7800 ICP mass spectrometer (Agilent, Santa Clara, CA, USA). The working conditions were as follows: reference power 1550 W, plasma airflow 15 L/min, auxiliary airflow 0.9 L/min, nebulizer airflow 0.99 L/min. Certified standard solutions of 1g/L concentration (Panreac, Barcelona, Spain) were used, grouping the elements into two multi-elemental standards. The results were expressed in mg kg−1 of sample for As, Cu, Hg and Pb and in g kg−1 for Al.

2.3. Near Infrared Spectroscopy

A total of 10 g of each of 322 samples (7 black, 11 green, and 8 red teas and the binary mixtures of black and red, black and green, and green and red tea prepared ranging from 50 to 99% of the tea of interest) were analyzed using the Foss NIRSystem 500 (Hillerod, Denmark) with the transport half-cup which has a 11.6 × 3.7 cm window. Spectra were recorded from 1000–2498 nm every 2 nm and 32 spectra were performed at each recording. Three spectra of each sample were recorded and then the average spectrum was calculated. The diffuse reflectance signal of the NIR spectrum was expressed as log (1/R) (R = reflectance). The software used for data manipulation and chemometric analysis was Win ISI 4.10. In order to develop the models 80% of the samples (257) were taken in the calibration group and 20% of the samples (65) in that of external validation, the samples were place in the two groups at random.

2.4. Chemometric Techniques

The metal content was quantified by the modified partial least squares (MPLS) regression method [45]. In MPLS the NIR residuals at each wavelength were obtained after each factor had been calculated and standardized (by dividing the standard deviations of the residuals at each wavelength) before calculating the next factor. To minimize the undesired contributions to the NIR signals, various spectral pretreatments were used together with derivatives and smoothing. The effects of scattering were removed using multiplicative scatter correction (MSC), standard normal variate (SNV), detrend (DT), or SNV-DT and combinations of these [28,33,46]. Moreover, several mathematical treatments were tested in the development of the NIRS calibrations where each of one was coded with four digits (2,4,4,1, for example), in which the first digit is the number of the derivative, the second is the gap over which the derivative is calculated, the third is the number of data points in a running average or smoothing, and the fourth is the second smoothing [33]. The number of variables was reduced by principal component analysis (PCA). A sample with an H statistic and ≥3.0 standardized units from the mean spectrum was defined as a global H outlier and eliminated from the calibration set. The calibration process was implemented with the spectra of the resulting samples and their chemical data. Using the criterion T ≥ 2.5, samples with high residual values when predicted were removed from the set. The statistical parameters of the calibration were obtained for each of the components after removing the samples for spectral (H criterion) or chemical (T criterion) reasons. To avoid overfitting during the development of the equations using MPLS, it is recommended to apply a cross-validation to select the optimal number of factors [47]. Therefore, the calibration set is divided into groups, in this case four due to the number of samples. Then, each group is validated using a calibration developed on the other samples. Finally, validation errors are combined into a standard error of cross-validation (SECV) [48]. The SECV has been reported to be the single best estimate of the predictive ability of the equation. This statistic is similar to the mean standard error of prediction (SEP) of ten randomly chosen prediction sets [49]. The statistics R squared (RSQ, multiple correlation coefficients) and SECV were used to select the best equations.

2.5. Pattern Recognition

Unsupervised pattern recognition by means of Principal Component Analysis (PCA) was performed to highlight the heavy metals which could be markers to discriminate the different teas (green, black, or red). The calculations were performed using the heavy metal results of the pure samples through the MATLAB software R2007B (The MathWorks Inc., Natick, MA, USA) and PLS-toolbox 5.2.

2.6. Discriminant Analysis

The samples were discriminated according to their purity by using different approaches and methods. The first approach was to study the feasibility of discriminating each type of the pure samples (black, red, and green) from each type of the following blends: black blends containing more than 50% of black tea; red blends containing more than 50% of red tea; and green blends containing more than 50% of green tea, which means that six groups were analyzed. The second approach was to study the feasibility of discriminating the percentage of the blend. To do so, four categories of each type of tea were considered: blends containing more than 95%, blends containing between 95% and 85%, blends containing between 85% and 75%, and blends containing between 75% and 50% of the selected tea (black, red, or green). Twelve categories were therefore considered in this approach
The SIMCA method (soft independent modelling of class analogy) was used to classify the samples using the whole NIR spectrum. The data (700 pieces of data corresponding to the spectrum from 1100 to 2498 nm, each measuring 2 nm) were normalized, scaled, and mean-centered for SIMCA modelling. First of all, a PCA was applied, and a new set of variables were calculated by linearly transforming the original variables into a set of new components (PCs) which retain the information of the original set. The number of PCs for classifying purposes was selected using the criteria of having an eigenvalue of >1. The projection of the samples in the space determined by the first PCs allowed the detection of groups present in the samples. Afterwards, the OPLS-DA (orthogonal projection latent structure discriminant analysis) method was carried out on the training set (80% of the data) and on the validation set (20% of the samples) in order to test the robustness of the discriminant model. The software used was SIMCA-P software version 14.1 (Umetrics, Sartorius Stedim Biotech AS, Umeå, Sweden).
The RMS-X residual (residual mean squares) method was also used. In this method the square root of the arithmetical mean of the residual mean squares (RMS or the quadratic mean) is calculated to obtain the RMS-X residual values. Different mathematical treatments, derivatives, and smoothing procedures coded as previously reported (i.e., Detrend (2,4,4,1)) were tested both for the calibration set (80% of the samples) and for the validation set (20% of the samples). The sensitivity (the proportion of samples belonging to one specific group which are correctly identified) and selectivity (the proportion of samples not belonging to that group which are correctly identified) were also calculated using the expressions proposed by [50]. The combination giving the highest percentage of correctly classified samples together with the higher selectivity and specificity was selected as the most suitable. The WinISI 4.10 software was used for this analysis.

3. Results and Discussion

3.1. The Heavy Metal Content of the Samples

The contents of the five heavy metals quantified by ICP-MS in pure samples of green, black, and red tea are shown in Table 1. Statistical analysis showed significant differences between samples for all of the heavy metals. Black tea samples showed significantly lower contents of Al, Pb, and As, while green tea showed the lowest values of cupper and was the only type of sample to show the presence of Hg. On the other hand, red tea showed the highest values of As and Pb.
The mean values observed for aluminum are above the range previously reported for black 0.4–1.2 g/kg [27,51,52,53], green 0.5–1.41 g/kg [53,54,55], and red tea 1.44 g/kg [26]. In accordance with our results, these studies highlight the fact that green and red tea showed higher values of aluminum than black tea [26,53]. As regards Pb, our results are fall within the ranges previously reported in the literature for black tea (0.15 to 2.31 mg/kg), green tea (0.10 to 1.93 mg/kg), and red tea (0.10 to 4.66 mg/kg) [9,21,55,56,57,58,59,60,61]. Moreover, the highest amounts of Pb were observed in red tea followed by green and black tea as reported in the aforementioned studies and others comparing the three types of tea [59]. Less research has been carried out on the arsenic and mercury content. The reported values for As varied between 0.09 and 0.39 mg/kg in black tea [57,62,63,64], between 0.038 and 0.088 mg/kg in green tea [55], and between 0.07 and 0.25 in red tea [9,60]. The levels of this heavy metal shown in Table 1 were within the range described for red tea and in the lower part of the range for black and green tea. As far as Hg is concerned, the reported intervals are 0.04–0.61, 0.032–0.22, and 0.01–0.03 mg/kg for black, green and red tea, respectively, [55,60,62,65]; the results found in the teas commercialized in Spain were therefore lower than previously reported. Finally, the results reported for copper show a high variability ranging from 11.29 to 44.38 mg/kg and from 4.7 to 36.5 mg/kg for black [27,57] and green tea [58], respectively. However, the most frequent values are between 10 and 26.5 mg/kg for both types of tea [21,27,53,55,57,61,62]. On the other hand, the reported values for red tea showed less variability and are within this range (14.34–21 mg/kg) [9,26,60,66,67] as was observed in the results obtained during this study. The current results also agreed with the studies of Gabrowska et al. [65], Konieczynsky et al. [67] Mackenzi et al. [26], and Czernicka et al. [66] who found that black tea showed higher levels of copper than green tea, while red tea tended to show similar values to black tea [26,65].

3.2. Exploratory Analysis of the Samples according to Their Heavy Metal Content

PCA was applied by means of the autoscale preprocess to extract more information from the heavy metal results obtained in the pure samples. Two major components (PCs) reflected 74.29% of the variance explained. The projection of the samples and the variables on the two first PC’s is shown in Figure 1.
PC2 reveals a distinction for the red tea samples on the positive side. As far as the loadings are concerned, samples of green tea are grouped by the Al and Hg results, while red tea is separated based on the Cu, As, and Pb results. In general, green tea presents higher concentrations of Al and Hg compared with red and black tea. On the other hand, the red tea samples showed higher concentrations of Cu, As, and Pb than green and black tea. The black tea samples presented a lower concentration of heavy metals. These results therefore indicated that it was possible to classify the tea types according to their heavy metal content.

3.3. Determination of the Heavy Metal Content Using near Infrared Spectroscopy

3.3.1. Spectra of the Samples

In order to perform the calibration equations, binary mixtures of black and red, black and green, and green and red tea were prepared ranging from 1% to 50% of each pair. The spectra of the pure teas (7 black, 11 green, and 8 red teas) and the mixtures were measured using the transport device of Foss NIRSystem 5000 (DK-3400, Hillerød, Denmark), equipped with the transport half-cup which has a 11.6 × 3.7 window. Figure 2a shows the spectra of the 322 samples directly obtained from the equipment, together with the average spectrum of the pure samples of each tea type, and Figure 2b the spectra after the pre-treatment Detrend (2,4,4,1) which gave the best calibration and validation results as was subsequently observed.
The main absorption peaks were found at 1440 and 1950 nm corresponding to O-H stretching and O-H deformation; information was therefore provided on the water content of the samples. According to the spectra, these peaks were larger for black tea followed by red and green tea. The peak founds at 1220 and 1670 nm corresponded to the C-H stretch second overtone and the C-H first overtone, while the peak found at 1740 corresponded to the -CH2 overtone. These bands have been linked to the presence of ingredients such as catechins, tea polyphenol, and alkaloids [68,69,70,71]. The region between 2210 and 2500 nm included a combination of C-H stretching and C-H deformation and a 2nd overtone of C-H deformation. This spectral region was linked to the presence of polyphenols, catechins, and teaflavins [72] but also to that of cellulose and caffeine [28,73]. Other peaks observed in the spectra were found at 2150 and 2170 nm and associated with NH2 vibrations in amide I groups + NH bending in amide III groups, which possibly corresponded to free amino acids and proteins [74]. These spectral features were even more clearly observed when a second derivative was applied (Figure 2b). In addition to the region between 2210 and 2500 nm or the region between 2100 and 2200 nm which showed a large number of peaks related to the phenolic and protein composition of the samples, respectively, already described, another spectral interval which showed a high number of peaks was the region between 1600 and 1800 nm corresponding to the CH-first overtone and the S-H first overtone of both total catechins and caffeine [60]. The peaks found between 1870 and 1900 corresponding to the C=O stretch overtone of TF2A or to the O-H combination and the C-H first overtone of TF2B are noteworthy [75]. Other peaks which can be observed more clearly in the second derivative are those at 1330 and 1400 nm which could be attributed to the first overtone of the C-H combinations of free amino acids [72].
The mean spectra of pure samples were similar to those previously described for black [46,69] green [74,76], and red teas [60] (Figure 2a). Black tea showed higher absorbance in all of the regions of the spectrum, while although the red and green tea spectrum showed similar absorbance, red tea always presented higher absorbance, especially in the band located at 1940 nm. This result suggests a higher water content in black tea followed by red and green tea. According to previous results the water content of different types of tea is highly variable and while some reports found a higher water content in black than in green teas [66], other studies reported the opposite [65] or similar values of water for black and green [77] or for black and red teas [65]. On the other hand, in general black teas showed lower levels of catechins and tannins and a higher caffeine content [65,67,77] because as the fermentation progresses the catechin content decreases and the caffeine increases [78]. As far as the protein content is concerned, it has been reported that it is very low in Purer tea [78] while the values for black and green tea showed a strong variability. In some cases, the values were similar for black and green tea and in others they were higher in green than in black tea [77,78,79].
Based on the mineral composition, a relationship between a particular absorption band and a particular mineral content cannot be established. This is due to the fact that there is no direct relationship between the mineral and the spectrum. Therefore, its determination is possible thanks to the associations of the minerals with other compounds present (this factor being different for each food).

3.3.2. Calibration Equations

The calibration was carried out by using the modified partial least squares regression (MPLS) method. To do so, the first step was to perform the principal component analysis (PCA). The results of the PCA showed that for all of the minerals using ten principal components the spectral variability explained was higher than 97%. In order to remove the samples which did not fulfil the spectral requirements, the Mahalanobis distance (H distance) was applied so that those samples with H values equal to or higher than 3 (1 sample) were removed from the sample set. Subsequently a calibration set including 80% of the samples (257) and a validation set with 20% of the samples (65) was built in which the mean, maximum values, and minimum values were similar for both sets (Table 2).
Table 3 shows the statistical descriptors of the best calibration equation obtained for each of the mineral elements (Al, As, Cu, Hg, and Pb). The best equation was determined according to RSQ, SEC, SECV, and RPD values. The final number of samples used for calibrations was the result of applying the T criterion. The samples with a T value exceeding 2.5 therefore had to be removed from the set because they were different from the population for chemical reasons. This number of samples depended on the component calibrated and ranged from 12 to 21. The calibration margins (the estimated minimum, the estimated maximum) were similar to or even broader than the range of values determined in the samples.
It can be observed that according to the results obtained with an RSQ of between 0.89 and 0.97 it is possible to predict all of the heavy metals in teas of varying compositions without any previous treatment or handling.

3.3.3. Validation

The internal validation of the model was carried out by means of the cross-validation method. In order to do so, the calibration set was divided into four cross-validation groups, of which three were taken for the calibration set and one for the prediction set. This process was repeated four times in order to allow all of the groups to pass through a calibration and a prediction set. With this step the prediction capacity of the model was checked. Figure 3 shows the correlation of the reference values (obtained by ICP-MS) and the predicted values using NIR spectra with RSQ values ranging from 0.820 for aluminum to 0.971 for copper.
The prediction capacity of the model was evaluated with the RPD parameter; a value greater than 2.5 indicated that the model can be considered adequate. All of the RPD values were above 2.5, ranging from 3.11 for aluminum to 5.78 for copper, which reveals that the prediction capacity of the model was suitable for predicting heavy metal contents in unknown samples of tea.
To check the robustness of the equations we applied them to the NIR spectra of the samples included in the external validation set (n = 65), the heavy metal composition of which was representative of the whole sample set (Table 2). Figure 4 shows the comparison between the values obtained after applying the equations and the reference values obtained by ICP-MS. The RSQ values ranged from 0.635 for lead to 0.929 for copper while the SEP and SEP(C) values were very similar for all of the heavy metals analyzed. The values obtained by both methodologies were compared using Student’s test-values for paired values. The p values (Table 4) were higher than 0.05 for all of the components; it should be pointed out that there were no differences between reference and predicted values. It can therefore be concluded that both analytical methods give similar results. Table 4 also shows the RMSE (root mean standard error) values which were low for all of the heavy metals.
These results are in accordance with those of previous studies which indicated that it was possible to predict the mineral content of different vegetable matrices with good RSQ reaching values of 0.98 for K and Na for example in alfalfa and flour, respectively, [30,80]. However, for the microelements which are found in smaller quantities in the samples, the RSQ values attained lower values ranging from 0.41 to 0.64 for Mg, 0.64 to 0.95 for Fe, and 0.41 to 0.84 for Mn in products such as alfalfa, rocket (Eruca versicaria), artichoke, quinoa, Curcubita pepo, and Brachiaria [31,34,81,82,83,84]. Data on the calibration of heavy metals with NIRS are scarce; it has been reported that copper can be predicted with RSQ values ranging from 0.52 to 0.76 [81,85], while RSQ values of between 0.70 and 0.79 have been reported for aluminum and lead [86]. NIR spectroscopy can detect macronutrients such as N, P, and S because they are major constituents of NIR-sensitive organic compounds, whereas micronutrients and macronutrients which exist mostly in inorganic forms such as Ca, Mg, and K are detected by means of association with organic compounds and indirect correlation with organic compounds [87]. Therefore, there is not a direct absorption peak of metal components in the NIR region [88].
However, according to the general lineal model, the predicted concentration of a mineral (Y) can be explained by the next equation Y = β0 + β1xλ1 + β2xλ2 + β3xλ3 +……+ βnxλn, in which β are the coefficients measuring the contribution of each wavelength and xλ1, xλ2, xλ3, ………. xλn, are the wavelengths. Therefore, the β coefficients can help to correlate the metal content with specific wavelengths which in turn are correlated with the bonds of different functional groups because the highest the β value the highest the contribution of this wavelength.
Figure 5 shows the β coefficients corresponding to the calibration equations of the five heavy metals calibrated. It can be observed that Al had two large peaks at 1282 and 1290 nm, the two highest peaks for Pb were found at 1164 and 1984 nm, As showed the highest values for β coefficients at 1264 and 1580 nm, Hg at 1244 and 1252 and finally Cu at 1310 and 1614 nm. Besides this, Figure 5 allows to observe that all of the peaks with the highest values, both negatives and positives, for each of the calibrated metals were in the spectral regions between 1110 and 1300 nm, in the region between 1500 and 1700 nm pointing out that these regions had a large influence on the developed models. The first one (1110–1300 nm), which corresponds with the C-H second overtone, can be attributed to the presence of the catechins, tea polyphenols and the TF2B theaflavin, while the region between 1500 and 1700 correspond to the O-H and N-H bond and can be related with the presence of catechins and caffeine. Other spectral regions that showed high values of β coefficients were sited at 1800 which could be due to C=O bond stretch second overtone of the TF2A theaflavin and O-H combination and C-H first overtone of the TF2B [72]. There was, therefore, a strong correlation between phenolic compounds and metal content. This was because phenolic compounds, especially flavonoids, tend to form stable complexes with metals such as iron, chromium, nickel, copper, or lead [89,90].
Previous work has shown high correlations between trace minerals and the spectral region 1100–1300 in forages and legumes as previously reported while [33] found that trace mineral concentration in lentils showed a strong correlation with absorbance in the range 1510–1550 nm. A significant correlation between minerals and the absorption bands near 1400 nm and the spectral region between 1600 and 1700 has also been described for the calibration of minerals in grape derivatives [85,88] as observed in tea. Similar results, although with differences in the specific wavelength where the absorption maximum appears, have been described for other vegetables such as rocket, pumpkin fruits or flours [80,81,82]. This is due to the fact that mineral trace elements are bound to complexes and chelates that appear to be different among the plant matrix.
Regarding the metals calibrated in this work, Al content in propolis has been correlated with absorption peaks at 1336 and 1556 nm while Cu showed correlation with β coefficient peaks at 1244 and 1356 nm [32], as were observed in this work for tea, which were attributed to phenolic content. In addition, [87] also found absorption peaks in the 2200–2400 range that were attributed to the affinity of copper for proteins and chloroplasts. With regard to Pb, this heavy metal showed two important peaks in tea, the previously mentioned at 1984 and another at 2000 nm in agreement with the values reported in propolis which were 1820 and 1968 [32] although another study showed that the bands at 548 and 718 nm were most likely related to the lead content in tea [91]. Arsenic levels also appear to be correlated with bands in the visible part of the spectrum (674 nm) and with the bands at 1916, 1956 and 1966 nm observed in prostrate amaranth [92]. These authors considered that as absorption affected cellular absorbers such as chlorophyll and other macromolecules such as proteins and lipids. These results are in agreement with the peak observed in Figure 5 at 1966 nm for As. The data related to mercury are very sparse and indicate that the wavelengths at 1011, 1173 and 1295 nm, which are in the same range reported in this work for tea, were particularly important for the estimation of Hg in soils [93]. Finally, it is noteworthy that Al, Pb, As and Cu showed a peak at 1290 nm with a high and positive value for Al and Cu but with a negative value for Pb and As. This opposite behavior for the same wavelength depending on the mineral was also found by Cozzolino et al. [88].

3.4. Discrimination of the Percentage of the Blends Using near Infrared Spectroscopy

The spectra of the pure tea samples and the different binary mixtures were analyzed in order to study the feasibility of discriminating each type of the pure samples (black, red and green teas) from their blends and the feasibility of discriminating the percentage of the blend for each type of tea. Two methods were assayed (SIMCA and RMS-X residuals) and for each of them the total of the samples was divided into a training set with 80% of the samples while the remaining 20% of the samples built the validation set.
The analysis of the whole spectra with the SIMCA method for discriminating the pure samples vs. their blends involved a PCA analysis of the samples which gave a result of three PCs with an eigenvalue of >1 which explained the 99.6% of the total variance. The representation of the samples in the space defined by the three PCs did not show a good separation between the six groups. Subsequently a classification by the OPLS-DA method was carried out and the results are shown in Table 5. It can be seen that 100% of the pure black samples were correctly classified for both the calibration and validation sets and 98% and 93% (calibration and validation) of the pure red samples were correctly classified. However, only 45% of the pure green samples were assigned to the correct class in the calibration step and this value fell to 10% for the validation group while green blends were better classified (77% and 75% in the calibration and validation groups). On the other hand, only 25% of the classifications were correct for the red blends in both sets while 0% of the black blends were correctly classified.
The OPLS-DA score plot of the six groups (Figure 6) showed that pure red, green, and black samples were well separated, but the blends (mainly black and red blends and green and red blends) appeared distributed along the axes connecting the pure samples.
The results of the OPLS-DA for the groups which differed as to the tea type (black, red, and green) and in the percentage of blending (>95%, 95–85%, 85–75% or 75–50%) are shown in Table 5. Once again, the discrimination results for pure or very high percentages of each tea sample were good. Therefore, 100% of the samples containing >95% of the black or green teas were correctly classified for both the calibration and validation sets, while for red teas the values were 91% and 88%, respectively. However, the model was not able to classify correctly any of the samples of the remainder of the blends.
As far as RMS-X residuals are concerned, this method implies a previous handling of the whole spectra with different combinations of mathematical treatments (MSC, SNV, DT or SNV-DT), derivatives, and smoothing procedures. After optimization the highest percentage of correctly classified samples was obtained by applying the Detrend (2,4,4,1) combination which was also the pre-treatment selected for heavy metal calibration. The results are shown in Table 5; a better performance was obtained by using RMS-X residuals than SIMCA for both analytical approaches, pure tea vs. blends, and the discrimination of the percentage of blending for each tea type.
As far as the classification of the teas according to the purity of the samples is concerned, 100% of the pure samples were correctly classified in calibration and only for pure red tea the classification was successful for 91% of the samples in validation. To consider the different blends, the percentage of correctly classified samples was between 93% for black blends and 98% for red blends in calibration and between 71% and 97% for black and green blends respectively. The high sensitivity and in particular the high specificity of the RMS-X residuals method is noteworthy in both the calibration and the validation sets, which means that only the samples belonging to a certain group are included in this group.
Finally, the classification according to the percentage of blending for each type of tea showed very good results for the classification stage with values ranging from 70% for blends containing between 75% and 50% of black tea to 100%. However, although the percentage of samples correctly classified in the validation set was also 100% for the blends containing >95% of each tea, the values were much lower for the remainder of the blends and even reached 0% for the 75–85% blends. One of the reasons is that there were 12 categories and the number of samples in some categories were extremely low. Despite this, the specificity of the method was again very high (0.87–1).
These results agree with those of previous studies which showed the feasibility of discriminating teas by using NIR technology. The treatment of NIR spectra with different chemometric tools has therefore allowed the successful discrimination of teas according to their quality [29,59], variety [68,69], degree of fermentation [94], geographical origin [15], and blending [41]. On the other hand, RMS-X residuals showed a satisfactory level of discrimination which is in agreement with previous studies which pointed out the feasibility of using this method to distinguish the samples according to their origin [95,96] and even according to the percentage of the mixture and variety [80]. Indeed, this method has shown better performance than other classification methods such as SIMCA as has previously been reported by Revilla et al. [97].

4. Conclusions

The results allowed us to conclude that it was possible to discriminate pure tea samples according to their content in heavy metals (As, Pb, Cu, Al, Hg). Red teas were therefore characterized by their higher content of As and Pb, green teas were the only ones containing Hg, and black teas showed the highest levels of copper.
Moreover, NIRS technology has been shown to be a suitable technique for determining the heavy metal composition of tea leaves. It was possible to predict all of the heavy metals analyzed with high RSQ (>0.89) and RPD (>3) values. The range of applicability of the equations is comparable to the reference method, and no significant differences between the reference and the predicted values were observed. Finally, NIR spectroscopy together with the RMS-X residual method allowed the discrimination of pure tea samples (black, red and green) and their blends with an accuracy of 98.3% in calibration while 92.3% of the total samples were correctly classified in validation. When the samples were discriminated according to the percentage of blending, 100% of the samples of 10 out of 12 groups were correctly classified in calibration although only the groups with a level of pure tea of >95% showed 100% of the samples to be correctly classified at the validation stage. NIR spectroscopy has therefore been shown to be a suitable technique for detecting the presence of other tea types in pure samples and for quantifying the heavy metal content in a multiparametric, fast, and simple way without prior treatment of the sample.

Author Contributions

Resources and Project administration A.M.V.-Q.; Methodology, A.M.V.-Q.; Investigation, M.R.-F. and P.V., Formal analysis, I.M.-M., M.H.J. and P.V. Writing-original draft, I.R. Writing-reviewing and editing, I.R. and A.M.V.-Q.; visualization, M.R.-F.; supervision, I.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research has been carried out as part of the “Shuttle” framework program (LANZADERA_TCUE_18-20P_002), administered by the Fundación General de la Universidad de Salamanca as part of the TCUE Plan 2018–2020 and has been co-funded by the European Regional Development Fund (ERDF) and the Junta de Castilla y León.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Acknowledgments

M. Hernández Jimenez acknowledges the Own Program III: Grants for Pre-doctoral Contracts of the University of Salamanca co-funded by Banco Santander.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Vuo, H. Tea: Market Data & Analysis. Available online: https://www.statista.com/study/48824/tea-market-data-and-analysis/16-09-2023 (accessed on 20 December 2023).
  2. Zhang, L.; Ho, C.-T.; Zhou, J.; Santos, J.S.; Armstrong, L.; Granato, D. Chemistry and Biological Activities of Processed Camellia sinensis Teas: A Comprehensive Review. Compr. Rev. Food Sci. Food Saf. 2019, 18, 1474–1495. [Google Scholar] [CrossRef] [PubMed]
  3. Karori, S.M.; Wachira, F.N.; Wanyoko, J.K.; Ngure, R.M. Antioxidant Capacity of Different Types of Tea Products. Afr. J. Biotechnol. 2007, 6, 2287–2296. [Google Scholar] [CrossRef]
  4. Lim, C.M.; Carey, M.; Williams, P.N.; Koidis, A. Rapid Classification of Commercial Teas According to Their Origin and Type Using Elemental Content with X-Ray Fluorescence (XRF) Spectroscopy. Curr. Res. Food Sci. 2021, 4, 45–52. [Google Scholar] [CrossRef] [PubMed]
  5. Chang, C.-W.; Wang, S.-H.; Jan, M.-Y.; Wang, W.-K. Effect of Black Tea Consumption on Radial Blood Pulse Spectrum and Cognitive Health. Complement. Ther. Med. 2017, 31, 1–7. [Google Scholar] [CrossRef] [PubMed]
  6. Abd El-Aty, A.M.; Choi, J.-H.; Rahman, M.M.; Kim, S.-W.; Tosun, A.; Shim, J.-H. Residues and Contaminants in Tea and Tea Infusions: A Review. Food Addit. Contam. Part A 2014, 31, 1794–1804. [Google Scholar] [CrossRef]
  7. Wang, Y.-J.; Li, T.-H.; Li, L.-Q.; Ning, J.-M.; Zhang, Z.-Z. Micro-NIR Spectrometer for Quality Assessment of Tea: Comparison of Local and Global Models. Spectrochim. Acta A Mol. Biomol. Spectrosc. 2020, 237, 118403. [Google Scholar] [CrossRef]
  8. Lv, H.; Zhang, Y.; Lin, Z.; Liang, Y. Processing and Chemical Constituents of Pu-Erh Tea: A Review. Food Res. Int. 2013, 53, 608–618. [Google Scholar] [CrossRef]
  9. Lv, H.; Yang, T.; Ma, C.; Wang, C.; Shi, J.; Zhang, Y.; Peng, Q.; Tan, J.; Guo, L.; Lin, Z. Analysis of Naturally Occurring 3″-Methyl-Epigallocatechin Gallate in 71 Major Tea Cultivars Grown in China and Its Processing Characteristics. J. Funct. Foods 2014, 7, 727–736. [Google Scholar] [CrossRef]
  10. Xu, J.; Wang, M.; Zhao, J.; Wang, Y.-H.; Tang, Q.; Khan, I.A. Yellow Tea (Camellia sinensis L.), a Promising Chinese Tea: Processing, Chemical Constituents and Health Benefits. Food Res. Int. 2018, 107, 567–577. [Google Scholar] [CrossRef]
  11. Desideri, D.; Meli, M.A.; Roselli, C.; Feduzi, L. Polarized X Ray Fluorescence Spectrometer (EDPXRF) for the Determination of Essential and Non Essential Elements in Tea. Microchem. J. 2011, 98, 186–189. [Google Scholar] [CrossRef]
  12. Polechońska, L.; Dambiec, M.; Klink, A.; Rudecki, A. Concentrations and Solubility of Selected Trace Metals in Leaf and Bagged Black Teas Commercialized in Poland. J. Food Drug Anal. 2015, 23, 486–492. [Google Scholar] [CrossRef]
  13. Karak, T.; Bhagat, R.M. Trace Elements in Tea Leaves, Made Tea and Tea Infusion: A Review. Food Res. Int. 2010, 43, 2234–2252. [Google Scholar] [CrossRef]
  14. Başgel, S.; Erdemoğlu, S.B. Determination of Mineral and Trace Elements in Some Medicinal Herbs and Their Infusions Consumed in Turkey. Sci. Total Environ. 2006, 359, 82–89. [Google Scholar] [CrossRef]
  15. Yan, S.-M.; Liu, J.-P.; Xu, L.; Fu, X.-S.; Cui, H.-F.; Yun, Z.-Y.; Yu, X.-P.; Ye, Z.-H. Rapid Discrimination of the Geographical Origins of an Oolong Tea (Anxi-Tieguanyin) by Near-Infrared Spectroscopy and Partial Least Squares Discriminant Analysis. J. Anal. Methods Chem. 2014, 2014, e704971. [Google Scholar] [CrossRef]
  16. Heaney, S.; Koidis, A.; Morin, J.M. Tea and Flavoured Tea. In Handbook of Food Authenticity: A Guide to Food Authenticity Issues and Analytical Solutions; Morin, J.-F., Lees, M., Eds.; Elsevier: Amsterdam, The Netherlands, 2018; pp. 315–331. [Google Scholar]
  17. Zhao, H.; Yu, C.; Li, M. Effects of Geographical Origin, Variety, Season and Their Interactions on Minerals in Tea for Traceability. J. Food Compos. Anal. 2017, 63, 15–20. [Google Scholar] [CrossRef]
  18. Chen, Y.; Xu, J.; Yu, M.; Chen, X.; Shi, J. Lead Contamination in Different Varieties of Tea Plant (Camellia sinensis L.) and Factors Affecting Lead Bioavailability. J. Sci. Food Agric. 2010, 90, 1501–1507. [Google Scholar] [CrossRef]
  19. Karak, T.; Bora, K.; Paul, R.K.; Das, S.; Khare, P.; Dutta, A.K.; Boruah, R.K. Paradigm Shift of Contamination Risk of Six Heavy Metals in Tea (Camellia sinensis L.) Growing Soil: A New Approach Influenced by Inorganic and Organic Amendments. J. Hazard. Mater. 2017, 338, 250–264. [Google Scholar] [CrossRef]
  20. Rashid, M.H.; Fardous, Z.; Chowdhury, M.A.Z.; Alam, M.K.; Bari, M.L.; Moniruzzaman, M.; Gan, S.H. Determination of Heavy Metals in the Soils of Tea Plantations and in Fresh and Processed Tea Leaves: An Evaluation of Six Digestion Methods. Chem. Cent. J. 2016, 10, 7. [Google Scholar] [CrossRef] [PubMed]
  21. Gallaher, R.N.; Gallaher, K.; Marshall, A.J.; Marshall, A.C. Mineral Analysis of Ten Types of Commercially Available Tea. J. Food Compos. Anal. 2006, 19, S53–S57. [Google Scholar] [CrossRef]
  22. Nkansah, M.A.; Opoku, F.; Ackumey, A.A. Risk Assessment of Mineral and Heavy Metal Content of Selected Tea Products from the Ghanaian Market. Environ. Monit. Assess. 2016, 188, 332. [Google Scholar] [CrossRef]
  23. Martín-Domingo, M.C.; Pla, A.; Hernández, A.F.; Olmedo, P.; Navas-Acien, A.; Lozano-Paniagua, D.; Gil, F. Determination of Metalloid, Metallic and Mineral Elements in Herbal Teas. Risk Assessment for the Consumers. J. Food Compos. Anal. 2017, 60, 81–89. [Google Scholar] [CrossRef]
  24. Drivelos, S.A.; Georgiou, C.A. Multi-Element and Multi-Isotope-Ratio Analysis to Determine the Geographical Origin of Foods in the European Union. TrAC Trends Anal. Chem. 2012, 40, 38–51. [Google Scholar] [CrossRef]
  25. Atasoy, A.D.; Yesilnacar, M.I.; Yildirim, A.; Atasoy, A.F. Nutritional Minerals and Heavy Metals in Tea Infusions and Daily Intake of Human Body. Turk. J. Agric.-Food Sci. Technol. 2019, 7, 234–239. [Google Scholar] [CrossRef]
  26. McKenzie, J.S.; Jurado, J.M.; De Pablos, F. Characterisation of Tea Leaves According to Their Total Mineral Content by Means of Probabilistic Neural Networks. Food Chem. 2010, 123, 859–864. [Google Scholar] [CrossRef]
  27. Salahinejad, M.; Aflaki, F. Toxic and Essential Mineral Elements Content of Black Tea Leaves and Their Tea Infusions Consumed in Iran. Biol. Trace Elem. Res. 2010, 134, 109–117. [Google Scholar] [CrossRef] [PubMed]
  28. Lee, M.-S.; Hwang, Y.-S.; Lee, J.; Choung, M.-G. The Characterization of Caffeine and Nine Individual Catechins in the Leaves of Green Tea (Camellia sinensis L.) by near-Infrared Reflectance Spectroscopy. Food Chem. 2014, 158, 351–357. [Google Scholar] [CrossRef] [PubMed]
  29. Dong, C.; Li, J.; Wang, J.; Liang, G.; Jiang, Y.; Yuan, H.; Yang, Y.; Meng, H. Rapid Determination by near Infrared Spectroscopy of Theaflavins-to-Thearubigins Ratio during Congou Black Tea Fermentation Process. Spectrochim. Acta A Mol. Biomol. Spectrosc. 2018, 205, 227–234. [Google Scholar] [CrossRef]
  30. González-Martín, I.; Hernández-Hierro, J.M.; González-Cabrera, J.M. Use of NIRS Technology with a Remote Reflectance Fibre-Optic Probe for Predicting Mineral Composition (Ca, K, P, Fe, Mn, Na, Zn), Protein and Moisture in Alfalfa. Anal. Bioanal. Chem. 2007, 387, 2199–2205. [Google Scholar] [CrossRef]
  31. González Martín, M.I.; Wells Moncada, G.; Fischer, S.; Escuredo, O. Chemical Characteristics and Mineral Composition of Quinoa by Near-infrared Spectroscopy. J. Sci. Food Agric. 2014, 94, 876–881. [Google Scholar] [CrossRef]
  32. González-Martín, M.I.; Escuredo, O.; Revilla, I.; Vivar-Quintana, A.M.; Coello, M.C.; Riocerezo, C.P.; Moncada, G.W. Determination of the Mineral Composition and Toxic Element Contents of Propolis by Near Infrared Spectroscopy. Sensors 2015, 15, 27854–27868. [Google Scholar] [CrossRef]
  33. Lastras, C.; Revilla, I.; González-Martín, M.I.; Vivar-Quintana, A.M. Prediction of Fatty Acid and Mineral Composition of Lentils Using near Infrared Spectroscopy. J. Food Compos. Anal. 2021, 102, 104023. [Google Scholar] [CrossRef]
  34. González-Martín, M.I.; Hernández-Hierro, J.M.; Revilla, I.; Vivar-Quintana, A.; González-Pérez, C.; García, L.G.; Riocerezo, C.P.; Ortega, I.A.L. Differentiation of Organic and Non-Organic Ewe’s Cheeses Using Main Mineral Composition or near Infrared Spectroscopy Coupled to Chemometric Tools: A Comparative Study. Talanta 2011, 85, 1915–1919. [Google Scholar] [CrossRef]
  35. González-Martín, I.; Hernández-Hierro, J.M.; Bustamante-Rangel, M.; Barros-Ferreiro, N. Near-Infrared Spectroscopy (NIRS) Reflectance Technology for the Determination of Tocopherols in Alfalfa. Anal. Bioanal. Chem. 2006, 386, 1553–1558. [Google Scholar] [CrossRef]
  36. Turgut, S.S.; Küçüköner, E.; Karacabey, E. TeaPot: A Chemometric Tool for Tea Blend Recipe Estimation. Appl. Food Res. 2021, 1, 100006. [Google Scholar] [CrossRef]
  37. Gogoi, R.C. Blending of Tea—The Development. Two Bud 2014, 61, 53–56. [Google Scholar]
  38. Rehder, A.P.; Silva, P.B.; Xavier, A.M.F.; Barrozo, M.A.S. Optimization of Microwave-Assisted Extraction of Bioactive Compounds from a Tea Blend. J. Food Meas. Charact. 2021, 15, 1588–1598. [Google Scholar] [CrossRef]
  39. Builders, P.F.; Mohammed, B.B.; Sule, Y.Z. Preparation and Characterization of a Poly-Herbal Tea with Effective Antioxidant Properties. Sci. World J. 2020, 15, 29–34. [Google Scholar]
  40. Sandasi, M.; Chen, W.; Vermaak, I.; Viljoen, A. Non-Destructive Quality Assessment of Herbal Tea Blends Using Hyperspectral Imaging. Phytochem. Lett. 2018, 24, 94–101. [Google Scholar] [CrossRef]
  41. Xu, S.; Sun, X.; Lu, H.; Zhang, Q. Detection of Type, Blended Ratio, and Mixed Ratio of Pu’er Tea by Using Electronic Nose and Visible/Near Infrared Spectrometer. Sensors 2019, 19, 2359. [Google Scholar] [CrossRef] [PubMed]
  42. Higa, F.; Koppel, K.; Chambers, E. Effect of Additional Information on Consumer Acceptance: An Example with Pomegranate Juice and Green Tea Blends. Beverages 2017, 3, 30. [Google Scholar] [CrossRef]
  43. Vilà, M.; Bedmar, À.; Saurina, J.; Núñez, O.; Sentellas, S. High-Throughput Flow Injection Analysis–Mass Spectrometry (FIA-MS) Fingerprinting for the Authentication of Tea Application to the Detection of Teas Adulterated with Chicory. Foods 2022, 11, 2153. [Google Scholar] [CrossRef]
  44. Cardoso, V.G.K.; Poppi, R.J. Non-Invasive Identification of Commercial Green Tea Blends Using NIR Spectroscopy and Support Vector Machine. Microchem. J. 2021, 164, 106052. [Google Scholar] [CrossRef]
  45. Martens, H.; Naes, T. Multivariate Calibration; Wiley: Hoboken, NJ, USA, 2001; Volume 116. [Google Scholar]
  46. Turgut, S.S.; Entrenas, J.A.; Taşkın, E.; Garrido-Varo, A.; Pérez-Marín, D. Estimation of the Sensory Properties of Black Tea Samples Using Non-Destructive near-Infrared Spectroscopy Sensors. Food Control 2022, 142, 109260. [Google Scholar] [CrossRef]
  47. Shenk, J.; Westerhaus, M. Routine Operation, Calibration, Development and Network System Management Manual; NIRSystems Inc.: Silver Spring, MD, USA, 1995; Volume 3. [Google Scholar]
  48. Williams, P.C.; Sobering, D. How do we do it: A brief summary of the methods we use in developing near infrared calibrations. In Near Infrared Spectroscopy: The Future Waves; Davies, A.M.C., Williams, P.C., Eds.; NIR Publications: Chichester, UK, 1996; pp. 185–188. [Google Scholar]
  49. Williams, P.; Norris, K. Near-Infrared Technology in the Agricultural and Food Industries; Amercian Association of Cereal Chemists, Inc.: St. Paul, MN, USA, 1987. [Google Scholar]
  50. Horcada, A.; Valera, M.; Juárez, M.; Fernández-Cabanás, V.M. Authentication of Iberian Pork Official Quality Categories Using a Portable near Infrared Spectroscopy (NIRS) Instrument. Food Chem. 2020, 318, 126471. [Google Scholar] [CrossRef]
  51. Dambiec, M.; Polechońska, L.; Klink, A. Levels of Essential and Non-Essential Elements in Black Teas Commercialized in Poland and Their Transfer to Tea Infusion. J. Food Compos. Anal. 2013, 31, 62–66. [Google Scholar] [CrossRef]
  52. Wróbel, K.; Wróbel, K.; Urbina, E.M.C. Determination of Total Aluminum, Chromium, Copper, Iron, Manganese, and Nickel and Their Fractions Leached to the Infusions of Black Tea, Green Tea, Hibiscus Sabdariffa, and Ilex Paraguariensis (Mate) by ETA-AAS. Biol. Trace Elem. Res. 2000, 78, 271–280. [Google Scholar] [CrossRef] [PubMed]
  53. Karimi, G.; Hasanzadeh, M.; Nili, A.; Khashayarmanesh, Z.; Samiei, Z.; Nazari, F.; Teimuri, M. Concentrations and Health Risk of Heavy Metals in Tea Samples Marketed in Iran. Pharmacology 2008, 3, 164–174. [Google Scholar]
  54. Ye, Y.; Yan, W.; Peng, L.; He, J.; Zhang, N.; Zhou, J.; Cheng, S.; Cai, J. Minerals and Bioactive Components Profiling in Se-Enriched Green Tea and the Pearson Correlation with Se. LWT 2023, 175, 114470. [Google Scholar] [CrossRef]
  55. Islam, M.A.; Ebihara, M. Elemental Characterization of Japanese Green Tea Leaves and Tea Infusion Residue by Neutron-Induced Prompt and Delayed Gamma-Ray Analysis. Arab. J. Chem. 2017, 10, S677–S682. [Google Scholar] [CrossRef]
  56. Pourramezani, F.; Akrami Mohajeri, F.; Salmani, M.H.; Dehghani Tafti, A.; Khalili Sadrabad, E. Evaluation of Heavy Metal Concentration in Imported Black Tea in Iran and Consumer Risk Assessments. Food Sci. Nutr. 2019, 7, 4021–4026. [Google Scholar] [CrossRef]
  57. Tsushida, T.; Takeo, T. Zinc, Copper, Lead and Cadmium Contents in Green Tea. J. Sci. Food Agric. 1977, 28, 255–258. [Google Scholar] [CrossRef] [PubMed]
  58. Gonzalez-Weller, D.; Rubio, C.; Gutiérrez, A.J.; Pérez, B.; Hernández-Sánchez, C.; Caballero, J.M.; Revert, C.; Hardisson, A. Dietary Content and Evaluation of Metals in Four Types of Tea (White, Black, Red and Green) Consumed by the Population of the Canary Islands. Pharm. Anal. Acta 2015, 6, 10. [Google Scholar] [CrossRef]
  59. Li, F.; Lu, Q.; Li, M.; Yang, X.; Xiong, C.; Yang, B. Comparison and Risk Assessment for Trace Heavy Metals in Raw Pu-Erh Tea with Different Storage Years. Biol. Trace Elem. Res. 2020, 195, 696–706. [Google Scholar] [CrossRef] [PubMed]
  60. Liu, Z.; Xie, H.; Chen, L.; Huang, J. An Improved Weighted Partial Least Squares Method Coupled with Near Infrared Spectroscopy for Rapid Determination of Multiple Components and Anti-Oxidant Activity of Pu-Erh Tea. Molecules 2018, 23, 1058. [Google Scholar] [CrossRef] [PubMed]
  61. Karimi, E.; Oskoueian, E.; Hendra, R.; Oskoueian, A.; Jaafar, H.Z.E. Phenolic Compounds Characterization and Biological Activities of Citrus Aurantium Bloom. Molecules 2012, 17, 1203–1218. [Google Scholar] [CrossRef] [PubMed]
  62. Barman, T.; Barooah, A.K.; Goswami, B.C.; Sharma, N.; Panja, S.; Khare, P.; Karak, T. Contents of Chromium and Arsenic in Tea (Camellia sinensis L.): Extent of Transfer into Tea Infusion and Health Consequence. Biol. Trace Elem. Res. 2020, 196, 318–329. [Google Scholar] [CrossRef] [PubMed]
  63. Portugal, Ó.B.; Flores-Quispe, M. Concentrations of Arsenic, Cadmium, and Lead in Herbal Infusion Tea Bags Marketed in Tacna, Peru. Environ. Monit. Assess. 2022, 194, 534. [Google Scholar] [CrossRef] [PubMed]
  64. Falahi, E.; Hedaiati, R. Heavy Metal Content of Black Teas Consumed in Iran. Food Addit. Contam. Part B 2013, 6, 123–126. [Google Scholar] [CrossRef]
  65. Garbowska, B.; Wieczorek, J.; Garbowska, B.; Wieczorek, Z.; Polak-Śliwińska, M. The Content of Minerals, Bioactive Compounds and Anti-Nutritional Factors in Tea Infusions. J. Elem. 2017, 23, 369–380. [Google Scholar] [CrossRef]
  66. Czernicka, M.; Zagula, G.; Bajcar, M.; Saletnik, B.; Puchalski, C. Study of Nutritional Value of Dried Tea Leaves and Infusions of Black, Green and White Teas from Chinese Plantations. Rocz. Państw. Zakładu Hig. 2017, 68, 237–245. [Google Scholar]
  67. Konieczynski, P.; Viapiana, A.; Wesolowski, M. Comparison of Infusions from Black and Green Teas (Camellia sinensis L. Kuntze) and Erva-Mate (Ilex Paraguariensis A. St.-Hil.) Based on the Content of Essential Elements, Secondary Metabolites, and Antioxidant Activity. Food Anal. Methods 2017, 10, 3063–3070. [Google Scholar] [CrossRef]
  68. Wang, X.; Huang, J.; Fan, W.; Lu, H. Identification of Green Tea Varieties and Fast Quantification of Total Polyphenols by Near-Infrared Spectroscopy and Ultraviolet-Visible Spectroscopy with Chemometric Algorithms. Anal. Methods 2015, 7, 787–792. [Google Scholar] [CrossRef]
  69. Chen, Q.; Zhao, J.; Zhang, H.; Wang, X. Feasibility Study on Qualitative and Quantitative Analysis in Tea by near Infrared Spectroscopy with Multivariate Calibration. Anal. Chim. Acta 2006, 572, 77–84. [Google Scholar] [CrossRef] [PubMed]
  70. Bian, M.; Skidmore, A.K.; Schlerf, M.; Fei, T.; Liu, Y.; Wang, T. Reflectance Spectroscopy of Biochemical Components as Indicators of Tea (Camellia sinensis) Quality. Photogramm. Eng. Remote Sens. 2010, 76, 1385–1392. [Google Scholar] [CrossRef]
  71. Li, F.; Jiang, X.; Zhao, J.; Zhang, S. Graphene Oxide: A Promising Nanomaterial for Energy and Environmental Applications. Nano Energy 2015, 16, 488–515. [Google Scholar] [CrossRef]
  72. Chen, Q.; Chen, M.; Liu, Y.; Wu, J.; Wang, X.; Ouyang, Q.; Chen, X. Application of FT-NIR Spectroscopy for Simultaneous Estimation of Taste Quality and Taste-Related Compounds Content of Black Tea. J. Food Sci. Technol. 2018, 55, 4363–4368. [Google Scholar] [CrossRef] [PubMed]
  73. Panigrahi, N.; Bhol, C.S.; Das, B.S. Rapid Assessment of Black Tea Quality Using Diffuse Reflectance Spectroscopy. J. Food Eng. 2016, 190, 101–108. [Google Scholar] [CrossRef]
  74. Li, C.; Guo, H.; Zong, B.; He, P.; Fan, F.; Gong, S. Rapid and Non-Destructive Discrimination of Special-Grade Flat Green Tea Using Near-Infrared Spectroscopy. Spectrochim. Acta A Mol. Biomol. Spectrosc. 2019, 206, 254–262. [Google Scholar] [CrossRef]
  75. Cai, J.; Wang, Y.; Xi, X.; Li, H.; Wei, X. Using FTIR Spectra and Pattern Recognition for Discrimination of Tea Varieties. Int. J. Biol. Macromol. 2015, 78, 439–446. [Google Scholar] [CrossRef]
  76. Chen, Q.; Zhao, J.; Liu, M.; Cai, J.; Liu, J. Determination of Total Polyphenols Content in Green Tea Using FT-NIR Spectroscopy and Different PLS Algorithms. J. Pharm. Biomed. Anal. 2008, 46, 568–573. [Google Scholar] [CrossRef]
  77. Adnan, M.; Ahmad, A.; Ahmed, A.; Khalid, N.; Hayat, I.; Ahmed, I. Chemical Composition and Sensory Evaluation of Tea (Camellia sinensis) Commercialized in Pakistan. Pak. J. Bot. 2013, 45, 901–907. [Google Scholar]
  78. Lin, J.-K.; Lin, C.-L.; Liang, Y.-C.; Lin-Shiau, S.-Y.; Juan, I.-M. Survey of Catechins, Gallic Acid, and Methylxanthines in Green, Oolong, Pu-Erh, and Black Teas. J. Agric. Food Chem. 1998, 46, 3635–3642. [Google Scholar] [CrossRef]
  79. Horanni, R.; Engelhardt, U.H. Determination of Amino Acids in White, Green, Black, Oolong, Pu-Erh Teas and Tea Products. J. Food Compos. Anal. 2013, 31, 94–100. [Google Scholar] [CrossRef]
  80. Martínez-Martín, I.; Hernández-Jiménez, M.; Revilla, I.; Vivar-Quintana, A.M. Prediction of Mineral Composition in Wheat Flours Fortified with Lentil Flour Using NIR Technology. Sensors 2023, 23, 1491. [Google Scholar] [CrossRef] [PubMed]
  81. Villatoro-Pulido, M.; Moreno Rojas, R.; Muñoz-Serrano, A.; Cardeñosa, V.; Amaro López, M.Á.; Font, R.; Del Río-Celestino, M. Characterization and Prediction by Near-Infrared Reflectance of Mineral Composition of Rocket (Eruca Vesicaria Subsp. Sativa and Eruca Vesicaria Subsp. Vesicaria): Characterization and NIRS Prediction of Mineral Content of Rocket Germplasm. J. Sci. Food Agric. 2012, 92, 1331–1340. [Google Scholar] [CrossRef]
  82. Martínez-Valdivieso, D.; Font, R.; Gómez, P.; Blanco-Díaz, T.; Del Rïo-Celestino, M. Determining the Mineral Composition in Cucurbita Pepo Fruit Using near Infrared Reflectance Spectroscopy. J. Sci. Food Agric. 2014, 94, 3171–3180. [Google Scholar] [CrossRef]
  83. De Oliveira, D.M.; Fontes, L.M.; Pasquini, C. Comparing Laser Induced Breakdown Spectroscopy, near Infrared Spectroscopy, and Their Integration for Simultaneous Multi-Elemental Determination of Micro- and Macronutrients in Vegetable Samples. Anal. Chim. Acta 2019, 1062, 28–36. [Google Scholar] [CrossRef] [PubMed]
  84. Mir-Marqués, A.; Martínez-García, M.; Garrigues, S.; Cervera, M.L.; de la Guardia, M. Green Direct Determination of Mineral Elements in Artichokes by Infrared Spectroscopy and X-Ray Fluorescence. Food Chem. 2016, 196, 1023–1030. [Google Scholar] [CrossRef] [PubMed]
  85. Cozzolino, D.; Moron, A. Exploring the Use of near Infrared Reflectance Spectroscopy (NIRS) to Predict Trace Minerals in Legumes. Anim. Feed Sci. Technol. 2004, 111, 161–173. [Google Scholar] [CrossRef]
  86. Patel, N.; Toledo-Alvarado, H.; Cecchinato, A.; Bittante, G. Predicting the Content of 20 Minerals in Beef by Different Portable Near-Infrared (NIR) Spectrometers. Foods 2020, 9, 1389. [Google Scholar] [CrossRef]
  87. Prananto, J.A.; Minasny, B.; Weaver, T. Chapter One—Near Infrared (NIR) Spectroscopy as a Rapid and Cost-Effective Method for Nutrient Analysis of Plant Leaf Tissues. Adv. Agron. 2020, 164, 1–49. [Google Scholar] [CrossRef]
  88. Cozzolino, D.; Cynkar, W.; Shah, N.; Smith, P. Quantitative Analysis of Minerals and Electric Conductivity of Red Grape Homogenates by near Infrared Reflectance Spectroscopy. Comput. Electron. Agric. 2011, 77, 81–85. [Google Scholar] [CrossRef]
  89. Sharma, A.; Shahzad, B.; Rehman, A.; Bhardwaj, R.; Landi, M.; Zheng, B. Response of Phenylpropanoid Pathway and the Role of Polyphenols in Plants under Abiotic Stress. Molecules 2019, 24, 2452. [Google Scholar] [CrossRef] [PubMed]
  90. Havsteen, B.H. The Biochemistry and Medical Significance of the Flavonoids. Pharmacol. Ther. 2002, 96, 67–202. [Google Scholar] [CrossRef] [PubMed]
  91. Sanaeifar, A.; Zhang, W.; Chen, H.; Zhang, D.; Li, X.; He, Y. Study on Effects of Airborne Pb Pollution on Quality Indicators and Accumulation in Tea Plants Using Vis-NIR Spectroscopy Coupled with Radial Basis Function Neural Network. Ecotoxicol. Environ. Saf. 2022, 229, 113056. [Google Scholar] [CrossRef] [PubMed]
  92. Font, R.; Del Río, M.; Vélez, D.; Montoro, R.; De Haro, A. Use of Near-Infrared Spectroscopy for Determining the Total Arsenic Content in Prostrate Amaranth. Sci. Total Environ. 2004, 327, 93–104. [Google Scholar] [CrossRef]
  93. Xu, X.; Chen, S.; Ren, L.; Han, C.; Lv, D.; Zhang, Y.; Ai, F. Estimation of Heavy Metals in Agricultural Soils Using Vis-NIR Spectroscopy with Fractional-Order Derivative and Generalized Regression Neural Network. Remote Sens. 2021, 13, 2718. [Google Scholar] [CrossRef]
  94. Ning, J.; Sun, J.; Li, S.; Sheng, M.; Zhang, Z. Classification of Five Chinese Tea Categories with Different Fermentation Degrees Using Visible and Near-Infrared Hyperspectral Imaging. Int. J. Food Prop. 2017, 20, 1515–1522. [Google Scholar] [CrossRef]
  95. Hernández-Jiménez, M.; Hernández-Ramos, P.; Martínez-Martín, I.; Vivar-Quintana, A.M.; González-Martín, I.; Revilla, I. Comparison of Artificial Neural Networks and Multiple Regression Tools Applied to near Infrared Spectroscopy for Predicting Sensory Properties of Products from Quality Labels. Microchem. J. 2020, 159, 105459. [Google Scholar] [CrossRef]
  96. Hernández-Jiménez, M.; González-Martín, M.I.; Martínez-Martín, I.; Revilla, I.; Vivar-Quintana, A.M. Carbon Stable Isotopes, Fatty Acids and the Use of NIRS to Differentiate IBERIAN Pigs. Meat Sci. 2021, 182, 108619. [Google Scholar] [CrossRef]
  97. Revilla, I.; Vivar-Quintana, A.M.; González-Martín, M.I.; Hernández-Jiménez, M.; Martínez-Martín, I.; Hernández-Ramos, P. NIR Spectroscopy for Discriminating and Predicting the Sensory Profile of Dry-Cured Beef “Cecina”. Sensors 2020, 20, 6892. [Google Scholar] [CrossRef] [PubMed]
Figure 1. PCA results. (■) black () red. () green tea pure samples.
Figure 1. PCA results. (■) black () red. () green tea pure samples.
Foods 13 00450 g001
Figure 2. (a) Raw spectra of the 322 samples and mean NIR spectra of the black (-), red (-) and green (-) tea pure samples and (b) spectra after second derivative Detrend 2,4,4,1.
Figure 2. (a) Raw spectra of the 322 samples and mean NIR spectra of the black (-), red (-) and green (-) tea pure samples and (b) spectra after second derivative Detrend 2,4,4,1.
Foods 13 00450 g002
Figure 3. Comparison between the reference values determined by ICP-MS and the predicted values for the internal validation set (n = 257) for (a) aluminium, (b) lead, (c) arsenic, (d) mercury, (e) copper. SEP: standard error of prediction; SEP(C): standard error corrected by the bias; RSQ: multiple correlation coefficient.
Figure 3. Comparison between the reference values determined by ICP-MS and the predicted values for the internal validation set (n = 257) for (a) aluminium, (b) lead, (c) arsenic, (d) mercury, (e) copper. SEP: standard error of prediction; SEP(C): standard error corrected by the bias; RSQ: multiple correlation coefficient.
Foods 13 00450 g003
Figure 4. Comparison between the reference values determined by ICP-MS and the predicted values for the external validation set (n = 65) for (a) aluminium, (b) lead, (c) arsenic, (d) mercury, and (e) copper. SEP: standard error of prediction; SEP(C): standard error corrected by the bias; RSQ: multiple correlation coefficient.
Figure 4. Comparison between the reference values determined by ICP-MS and the predicted values for the external validation set (n = 65) for (a) aluminium, (b) lead, (c) arsenic, (d) mercury, and (e) copper. SEP: standard error of prediction; SEP(C): standard error corrected by the bias; RSQ: multiple correlation coefficient.
Foods 13 00450 g004
Figure 5. Values of the β coefficients for the calibration (a) aluminium, (b) lead, (c) arsenic, (d) mercury, and (e) copper.
Figure 5. Values of the β coefficients for the calibration (a) aluminium, (b) lead, (c) arsenic, (d) mercury, and (e) copper.
Foods 13 00450 g005aFoods 13 00450 g005b
Figure 6. The OPLS-DA score plot of the black, red and green pure teas and the black, red and green blends (blends with more than 50% of each tea type) (two predictive and two orthogonal components, R2X = 0.569, R2 = 0.314, Q2 = 0.318).
Figure 6. The OPLS-DA score plot of the black, red and green pure teas and the black, red and green blends (blends with more than 50% of each tea type) (two predictive and two orthogonal components, R2X = 0.569, R2 = 0.314, Q2 = 0.318).
Foods 13 00450 g006
Table 1. Total contents of heavy metal in pure black, green and red tea. Mean value and standard deviation (SD).
Table 1. Total contents of heavy metal in pure black, green and red tea. Mean value and standard deviation (SD).
NAl (g/kg)Pb (mg/kg)As (mg/kg)Hg (mg/kg)Cu (mg/kg)
Black tea72.04 ± 0.23 a0.56 ± 0.28 a0.04 ± 0.01 a0.00 ± 0.00 a17.31 ± 2.60 b
Red tea82.36 ± 0.25 b1.21 ± 0.36 b0.19 ± 0.09 b0.00 ± 0.00 a19.87 ± 2.78 b
Green tea112.52 ± 0.29 b0.98 ± 0.40 b0.09 ± 0.02 a0.01 ± 0.01 b12.19 ± 2.40 a
a,b Different letters in the same column indicate statistically significant differences; analysis of variance (ANOVA); post hoc: Tukey HSD test; p < 0.05.
Table 2. Mineral composition of the calibration and validation sets.
Table 2. Mineral composition of the calibration and validation sets.
Calibration Set (n = 257)External Validation Set (n = 65)
MineralMinimunMaximunMeanSDMinimunMaximunMeanSD
Al (g/kg)1.84903.04702.37050.2701.95402.95242.31530.223
Pb (mg/kg)0.28101.71000.89430.3330.31701.49520.90030.293
As (mg/kg)0.02910.39100.11010.0490.03050.22630.10990.048
Hg (mg/kg)0.00060.01890.00960.0040.00220.01850.00960.004
Cu (mg/kg)8.440023.914015.99843.6418.943023.720015.41393.728
Table 3. Statistical descriptors for the best calibration equations obtained from NIR spectra and mineral composition data obtained by ICP-MS on tea samples.
Table 3. Statistical descriptors for the best calibration equations obtained from NIR spectra and mineral composition data obtained by ICP-MS on tea samples.
MineralNMeanSDMin Est.Max Est.SECSECVRSQRPD
Al (g/kg)2392.35370.25651.5843.1230.0840.1090.8933.11
Pb (mg/kg)2430.87840.32850.0001.8640.0700.0910.9554.71
As (mg/kg)2350.10510.04440.0000.2380.0120.0180.9233.62
Hg (mg/kg)2440.00660.00570.0000.0240.0010.0020.9665.17
Cu (mg/kg)24316.09643.61195.26126.9320.6250.8500.9705.78
N: number of samples after removing the outliers; SD: standard deviation; SEC: standard error of calibration; SECV: standard error of cross validation; RSQ: multiple correlation coefficient; RPD: ratio performance deviation.
Table 4. External validation for the NIR determination of minerals in different types of tea.
Table 4. External validation for the NIR determination of minerals in different types of tea.
Mineralp (Level of Significance)RMSE
Al0.8900.112
Pb0.7921.003
As0.7780.182
Hg0.6840.027
Cu0.7960.021
Table 5. Results of the classification of the pure teas vs. tea blends and of the classification by percentage of blending by SIMCA and by RMSX-residuals method for the best mathematical pre-treatment applied to the NIR spectra.
Table 5. Results of the classification of the pure teas vs. tea blends and of the classification by percentage of blending by SIMCA and by RMSX-residuals method for the best mathematical pre-treatment applied to the NIR spectra.
SIMCA MethodRMSX-Residuals Detrend (2,4,4,1)
% Samples Correctly Classified% Samples Correctly ClassifiedSensitivitySpecificity
ComponentCalibrationValidationCalibrationValidationCalibrationValidationCalibrationValidation
Pure tea vs. blends
Black pure100.00100.00100.00100.001.001.001.001.00
Red pure98.4693.75100.0090.911.001.001.001.00
Green pure45.5610.00100.00100.001.000.911.001.00
Black blends0.000.0093.3371.430.930.710.991.00
Red blends24.0525.0098.7594.740.990.950.990.97
Green blends77.4475.0097.7996.970.980.971.000.91
Percentage of blending
>95% black100.00100.00100.00100.001.001.001.000.99
95–85% black0.000.00100.0050.001.000.501.001.00
85–75% black0.000.00100.000.001.000.001.001.00
75–50% black0.000.0070.0033.330.700.331.000.99
>95% red91.6788.24100.00100.001.001.001.000.95
95–85% red0.000.00100.0020.001.000.201.001.00
85–75% red0.000.00100.000.001.000.001.001.00
75–50% red0.000.00100.0066.671.000.970.990.92
>95% green100.00100.00100.00100.001.001.001.001.00
95–85% green0.000.00100.0062.501.000.971.001.00
85–75% green0.000.00100.000.001.000.961.000.99
75–50% green0.000.0098.5588.241.000.971.000.87
Black blends: blends with more than 50% of black tea with red or green tea. Red blends: blends with more than 50% of red tea with black or green tea. Green blends: blends with more than 50% of green tea with black or red tea.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Revilla, I.; Hernández Jiménez, M.; Martínez-Martín, I.; Valderrama, P.; Rodríguez-Fernández, M.; Vivar-Quintana, A.M. The Potential Use of Near Infrared Spectroscopy (NIRS) to Determine the Heavy Metals and the Percentage of Blends in Tea. Foods 2024, 13, 450. https://doi.org/10.3390/foods13030450

AMA Style

Revilla I, Hernández Jiménez M, Martínez-Martín I, Valderrama P, Rodríguez-Fernández M, Vivar-Quintana AM. The Potential Use of Near Infrared Spectroscopy (NIRS) to Determine the Heavy Metals and the Percentage of Blends in Tea. Foods. 2024; 13(3):450. https://doi.org/10.3390/foods13030450

Chicago/Turabian Style

Revilla, Isabel, Miriam Hernández Jiménez, Iván Martínez-Martín, Patricia Valderrama, Marta Rodríguez-Fernández, and Ana M. Vivar-Quintana. 2024. "The Potential Use of Near Infrared Spectroscopy (NIRS) to Determine the Heavy Metals and the Percentage of Blends in Tea" Foods 13, no. 3: 450. https://doi.org/10.3390/foods13030450

APA Style

Revilla, I., Hernández Jiménez, M., Martínez-Martín, I., Valderrama, P., Rodríguez-Fernández, M., & Vivar-Quintana, A. M. (2024). The Potential Use of Near Infrared Spectroscopy (NIRS) to Determine the Heavy Metals and the Percentage of Blends in Tea. Foods, 13(3), 450. https://doi.org/10.3390/foods13030450

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop