Compositional Nutrient Diagnosis Methodology and Its Effectiveness to Identify Nutrient Levels in Yerba Mate ( Ilex paraguariensis )

: Yerba mate is a forest species of both cultural and economic importance growing in the subtropical regions of South America, especially in the south of Brazil. Despite its importance, yerba mate has never received enough attention from researchers, so the nutritional sufficiency ranges and critical levels have not yet been determined. This research aimed to establish these parameters for yerba mate to enable its foliar diagnosis. A total of 167 leaf samples were collected from production fields located in the five yerba mate-growing regions in Rio Grande do Sul, and the leaf nutrients were determined by standard chemical methods. The yield of each production field was accessed, and the cutoff value separating low-and high-yield groups was calculated in 16.75 Mg ha − 1 . The multivariate compositional nutrient diagnosis (CND) standards were determined, and nutrient interactions were estimated by correlation and principal component analyses. There was no positive correlation between any single nutrient and yield, even in the high-yield population, evidencing that a higher yield is the outcome of the balance among all nutrients. Excess of B occurred in one-third of the low-yield samples, while deficiency of Cu and K occurred in one-fourth of these samples. Finally, we established the adequate leaf nutrient levels for yerba mate.


Introduction
Yerba mate (Ilex paraguariensis St. Hil.) is an endemic tree species of subtropical and temperate forest regions of South America, with wild occurrence restricted to northeastern Paraguay, the southern states of Brazil, the northeastern region of Argentina, and Uruguay [1].Its dried leaves are traditionally consumed in the form of an infusion, a tea-like drink that has been part of the local tradition since the pre-Columbian era [2].Later on, in the seventeenth century, yerba mate was already one of the most important items of trade in South America [3].Nowadays, in addition to the usual consumption, there is a growing demand for products derived from yerba mate.Mainly due to its various pharmacological properties, industrialized beverages, cosmetics, and medicines are increasingly being produced, both nationally and internationally [4].
In Rio Grande do Sul, the southernmost state of Brazil, yerba mate is one of the most important agricultural products, resulting in an average monthly revenue of more than USD 5 million [5].It is estimated that around 35,000 hectares are cultivated with yerba mate in Rio Grande do Sul, a large part of which is destined for export to countries such as Germany, Argentina, Uruguay, and the USA.Indeed, Rio Grande do Sul accounts for about 80% of the Brazilian yerba mate exports to the international market [6].
Despite the socioeconomic and cultural relevance of yerba mate in Rio Grande do Sul, there is a lack of technical information for the plant nutrient management.Unlike other important crops, the sufficiency ranges (SRs) and critical levels (CLs) for the interpretation of foliar analysis results for yerba mate have not yet been determined [7].It is important to notice that yerba mate is a perennial crop from which large amounts of plant biomass are removed during the harvesting process, making it necessary to quickly replace the nutrients in a balanced way [8,9].
In this context, the nutritional imbalance of plants and the lack of parameters to measure this imbalance are factors that can be contributing to the low productivity of yerba mate in Rio Grande do Sul.Among the different approaches to evaluating the nutritional balance of plants and defining CLs and SRs, compositional nutrient diagnosis (CND) stands out [8,10].The method is considered one of the best ways to express the balance, deficiency, or excess of nutrients in plant tissue, with advantages over the univariate critical value (CV) or the bivariate diagnosis and recommendation integrated system (DRIS) approaches [11].
Adopting CLs or SRs will contribute to defining when fertilizers must be applied to yerba mate.This approach helps minimize excess nutrient buildup in the soil, which in turn reduces the risk of contaminating the nearby surface and subsurface waters [12].Additionally, it can help prevent nutrient deficiencies in yerba mate plants, caused by the interaction of nutrients in the soil or solution, which can occur with greater intensity in soils with excess nutrients.
Thus, the cultivation of yerba mate requires parameters that help define the amount and form of replacement fertilization to ensure sustainable productivity in the long term.This research aimed to establish the nutritional sufficiency ranges for yerba mate to enable its foliar diagnosis.

Materials and Methods
Rio Grande do Sul state has five yerba mate-growing regions, known as Polos Ervateiros (Figure 1).The climate of these regions is humid subtropical, predominantly classified as Cfa and as Cfb in the higher altitudes of Alto Taquari, Nordeste Gaúcho, and Alto Uruguai, according to Köppen's classification [13].Although occupying a smaller territorial area, Cfb is the preferable climate for yerba mate [14].The altitude is between 400 and 1000 m [14,15].Soils are preponderantly Leptosols, Regosols, Cambisols, Chernosols, Alisols, and Ferralsols, being derived from basalts, in the lower altitudes, and dacites, in altitudes above 600 m [16].
as Germany, Argentina, Uruguay, and the USA.Indeed, Rio Grande do Sul accounts for about 80% of the Brazilian yerba mate exports to the international market [6].
Despite the socioeconomic and cultural relevance of yerba mate in Rio Grande do Sul, there is a lack of technical information for the plant nutrient management.Unlike other important crops, the sufficiency ranges (SRs) and critical levels (CLs) for the interpretation of foliar analysis results for yerba mate have not yet been determined [7].It is important to notice that yerba mate is a perennial crop from which large amounts of plant biomass are removed during the harvesting process, making it necessary to quickly replace the nutrients in a balanced way [8,9].
In this context, the nutritional imbalance of plants and the lack of parameters to measure this imbalance are factors that can be contributing to the low productivity of yerba mate in Rio Grande do Sul.Among the different approaches to evaluating the nutritional balance of plants and defining CLs and SRs, compositional nutrient diagnosis (CND) stands out [8,10].The method is considered one of the best ways to express the balance, deficiency, or excess of nutrients in plant tissue, with advantages over the univariate critical value (CV) or the bivariate diagnosis and recommendation integrated system (DRIS) approaches [11].
Adopting CLs or SRs will contribute to defining when fertilizers must be applied to yerba mate.This approach helps minimize excess nutrient buildup in the soil, which in turn reduces the risk of contaminating the nearby surface and subsurface waters [12].Additionally, it can help prevent nutrient deficiencies in yerba mate plants, caused by the interaction of nutrients in the soil or solution, which can occur with greater intensity in soils with excess nutrients.
Thus, the cultivation of yerba mate requires parameters that help define the amount and form of replacement fertilization to ensure sustainable productivity in the long term.This research aimed to establish the nutritional sufficiency ranges for yerba mate to enable its foliar diagnosis.

Materials and Methods
Rio Grande do Sul state has five yerba mate-growing regions, known as Polos Ervateiros (Figure 1).The climate of these regions is humid subtropical, predominantly classified as Cfa and as Cfb in the higher altitudes of Alto Taquari, Nordeste Gaúcho, and Alto Uruguai, according to Köppen's classification [13].Although occupying a smaller territorial area, Cfb is the preferable climate for yerba mate [14].The altitude is between 400 and 1000 m [14,15].Soils are preponderantly Leptosols, Regosols, Cambisols, Chernosols, Alisols, and Ferralsols, being derived from basalts, in the lower altitudes, and dacites, in altitudes above 600 m [16].A total of 167 leaf samples were collected from different production fields, located in the five Polos Ervateiros, from 2019 to 2022.On average, the soils where yerba mate was collected had a medium organic matter content, low pH, and high Mn and Al contents (Table 1).Leaf samples were dried in an air-forced oven at 60 • C until constant weight.Then samples were ground and sieved through a 2 mm mesh.Plant tissue chemical analyses were performed by standard methods [17].Acid digestion was performed with H 2 O 2 /H 2 SO 4 for the analysis of N, P, K, Ca, and Mg; and with HNO 3 /HClO 4 for the analysis of Zn, Cu, Mn, Fe, and S, in addition to burning in a muffle furnace for determination of B. The concentrations of Ca, Mg, Cu, Fe, Zn, and Mn were determined by atomic absorption spectrometry; P, B, and S were determined by visible-light spectrometry; N via NH 3 distillation; and K by flame spectrometry.Yerba mate yield, in Mg ha −1 , was informed by farmers, based on the previous harvest.

Parameter Mean Range
--------g kg −1 --------Clay 373 27 to 77 OM * 33 1.4 to 6.9 -------mg dm The step-by-step calculation and rationale of CND are fully detailed in Parent and Dafir [10], Traspadini et al. [18], Rozane et al. [19], Rodrigues et al. [20], Tadayon et al. [11], and González-Vences et al. [21].In short, the leaf tissue was considered a closed system formed by the sum of the 11 analyzed nutrients (N, P, K, . ..), in the same unit of measurement, and non-analyzed plant components (other undetermined elements, carbohydrates, . ..), gathered in a term called R. This forms a d-dimensional nutrient arrangement; that is, a simplex (S d ) (Equation ( 1)) arising from the proportions of d + 1 nutrients that include the d elements and a residual value (R d ) (Equation ( 2)): where N, P, K . . .are the proportions of nutrients determined in the dry matter, and Rd is calculated by difference, as follows: In our case, the nutrients were expressed in mg kg −1 , and the sum of all plant components totalled 1,000,000.The geometric mean was obtained by calculating the 12th root of the product of the 12 components of the leaf sample (11 analyzed nutrients + R).The multivariate relationship between plant components was obtained by the natural logarithm of the ratio of each component to the geometric mean of the sample.Nutrient proportions become scale invariant after they are divided by the geometric mean (G) (Equation ( 3)) of d + 1 components, including R d [22]: To express each component of the simplex in relation to all the others (interaction study), it is sufficient to define new variables (V) (Equation ( 4)), which undergo centered logarithmic transformation (natural or Neperian logarithm); that is, in relation to the geometric mean of the observed values, and are expressed as follows: and, by definition, V N + V P + V K + . . .+ V Rd = 0 Next, the database was divided into low-and high-yield groups, using the procedure described by Khiari et al. [23].CND norms consisted of the mean and standard deviation of the multivariate ratios of the nutrients in the high-yield group, which correspond to the relations of the centered logarithmic transformation of VX of d nutrients for high-yielding plants, i.e., V N *, V P *, V K * . . .V R * . . .and SD N *, SD P *, SD K *, . . .SD R *, respectively.
CND indices were calculated by subtracting the row-centered log of each plant component from the respective norm mean and then dividing the result by the norm standard deviation (Equation ( 5)).
where V X * and SD X * are, respectively, the mean and standard deviation of element X in the high-yield subpopulation, and I X is the CND index of element X.
Independence among the data is ensured by the centered logarithmic transformation [22].The CND indices are normed, and the variables made linear as dimensions of a circle (d + 1 = 2), a sphere (d + 1 = 3), or a large sphere (d + 1 > 3), in a dimensional space of d + 1.The nutritional imbalance index r 2 is distributed as a variable x d 2 if the CND indices are independent reduced variables [24], and it is calculated by the equation (Equation ( 6)): At last, the sufficiency range of each nutrient was determined by linear regression analysis between the nutrient content and the respective CND index [11,19].The critical level was obtained by equating the CND index to zero in the equation, and the sufficiency range was delimited by a 2/3 standard deviation around this value.
Nutrient interactions were estimated by correlation [19] and principal component analyses (PCAs) [11].Before performing the PCAs, our dataset underwent evaluation through the Kaiser-Meyer-Olkin (KMO) method and Bartlett's sphericity test utilizing the EFAtools package within R software (version 4.3.1).The KMO value was 0.519, and Bartlett's sphericity test was significant (p < 0.001; χ2 = 349.31),indicating that the data were suitable for factor analysis.All CND calculations were made using Microsoft Excel software (version Office 365).The correlation and t-test analyses were performed using SigmaPlot software (version 11).For the PCA analysis, we used PAST software (version 4.05), adopting the selection criterion (SC) proposed by Collins and Ovalles [25], by which eigenvectors for each PC were considered significant when having a value larger than that calculated using the formula SC = 0.5/(PC eigenvalue) 0.5 .

Results
The yield of the 167 yerba mate production fields from where the leaf samples were collected varied from 3 to 30, with an average of 13.07 Mg ha −1 .The means and ranges of the plant nutrient contents are shown in Table 2.The correlation analyses between the isolated nutrient content and the plant yield resulted in only one significant result: B was negatively correlated with it (Table 3).It was also possible to verify the significant influence of some single-nutrient leaf concentration on the concentration of other nutrients.The cutoff between the low-and high-yield groups was determined by establishing cubic functions between yerba mate yield and the 12 cumulative variance ratios and equating the second derivative of each equation to zero to find its inflection point, as described by Khiari et al. [23] (Table 4).The cutoff was established at 16.75 Mg ha −1 , the highest value within the database yield range.
Once the low-and high-yield groups were delimited, correlation analyses between isolated nutrients were performed again, now within each yield group.In the low-yield group, the yerba mate yield was negatively correlated with B, as verified for the whole population (Table 5).In addition, Fe and Mn were also negatively correlated with yield.In turn, in the high-yield group, the yield was negatively correlated with K and S (Table 6).Comparing the average leaf nutrient content between both groups, the only statistically significant difference was the higher level of B in the low-yield group (Table 7).This result reinforces the negative influence of B on the yerba mate yield.Nutrient interactions were also accessed via PCA.In the low-yield group, the first four PCs explained 68.65% of the total variation (Table 8).The A group, formed by K, P, and Cu, was positively loaded on PC1, while PC2 separated Mg from Zn and S. The third PC correlated negatively with P and Zn.Finally, the fourth PC showed a positive correlation with N. The first four PCs explained 74.1% of the total variance in the high-yield group (Table 9).The first PC clustered P, Ca, Mg, and B in the positive quadrant, while locating Zn in the negative quadrant.P, K, and Cu formed a significant cluster located in the positive quadrant of PC2, and located Mn in the negative quadrant.The third PC was positively correlated with S and B, and negatively with N, while the fourth PC was positively correlated with Fe and negatively with Zn.The CND norms as means and standard deviations of V N , V P , V K , V Ca , V Mg , V S , V B , V Cu , V Fe , V Mn , V Zn , and V R for the high-yield group, as well as the critical levels and sufficiency ranges for each nutrient are presented in Tables 10 and 11.On average, the frequencies of nutritional deficiency, balance, and excess were practically the same between low-(14.5, 67.5, and 15.0%, respectively) and high-yield (16.3, 67.9, and 15.8%, respectively) groups (Figure 2).However, it is interesting to notice that the excess of B occurred in 32.3% of the low-yield samples.The deficiency of K and Cu was observed in 27.1 and 31.6% of the low-yield samples, respectively.Furthermore, when we consider the average of the CND indices, we can see that B, K, and Cu were the nutrients with the highest imbalance, with indices above 0.25 in modulus (Figure 3).

Discussion
Establishing the sufficiency ranges and the optimal nutrient levels is crucial to evaluating the nutritional status of a crop.Among the different approaches, CND stands out [11].Many authors compared CND with other methods, such as DRIS and Mathematical Chance (MCh), and verified its efficiency and precision.Caires et al. [26] assessed the nutritional diagnosis of soybean leaves and seeds using DRIS, MCh, and CND.The authors observed that the optimal nutrient levels estimated by DRIS and CND closely matched the average content in the reference population, whereas the MCh method yielded slightly different values.Santos et al. [27] confirmed the same trend across these three methods in their study on sugar cane.Similar results between CND and DRIS were also observed for corn [28].CND and DRIS usually produce closely aligned results and are highly correlated [26].However, CND can produce narrower SRs and establish multivariate relationships among nutrients [28], allowing an accurate association between plant nutrition and yield [19].These results supported our decision to use the CND method for estimating the SRs and CLs of yerba mate.
As previously mentioned, the yield average of the 167 analyzed yerba mate production fields was 13.07 Mg ha −1 .This value is higher than the national and the Rio Grande do Sul state average yields in 2019, which were 7.7 and 10.0 Mg ha −1 , respectively [29].Although orchard yields above 12 Mg ha −1 can be considered high [30], yields higher than 20 Mg ha −1 are not unusual in orchards with high technological levels [31].Indeed, 68% of the evaluated yerba mate production fields were above the state average yield, while only 20% were considered in the high-yield group delimited by the cutoff established at 16.75 Mg ha −1 .These observations reinforce the idea that the Brazilian yield is much lower than the potential yield of yerba mate and has been decreasing in the last decades [29], showing the need for research and technological innovation for this productive chain.

Discussion
Establishing the sufficiency ranges and the optimal nutrient levels is crucial to evaluating the nutritional status of a crop.Among the different approaches, CND stands out [11].Many authors compared CND with other methods, such as DRIS and Mathematical Chance (MCh), and verified its efficiency and precision.Caires et al. [26] assessed the nutritional diagnosis of soybean leaves and seeds using DRIS, MCh, and CND.The authors observed that the optimal nutrient levels estimated by DRIS and CND closely matched the average content in the reference population, whereas the MCh method yielded slightly different values.Santos et al. [27] confirmed the same trend across these three methods in their study on sugar cane.Similar results between CND and DRIS were also observed for corn [28].CND and DRIS usually produce closely aligned results and are highly correlated [26].However, CND can produce narrower SRs and establish multivariate relationships among nutrients [28], allowing an accurate association between plant nutrition and yield [19].These results supported our decision to use the CND method for estimating the SRs and CLs of yerba mate.
As previously mentioned, the yield average of the 167 analyzed yerba mate production fields was 13.07 Mg ha −1 .This value is higher than the national and the Rio Grande do Sul state average yields in 2019, which were 7.7 and 10.0 Mg ha −1 , respectively [29].Although orchard yields above 12 Mg ha −1 can be considered high [30], yields higher than 20 Mg ha −1 are not unusual in orchards with high technological levels [31].Indeed, 68% of the evaluated yerba mate production fields were above the state average yield, while only 20% were considered in the high-yield group delimited by the cutoff established at 16.75 Mg ha −1 .These observations reinforce the idea that the Brazilian yield is much lower than the potential yield of yerba mate and has been decreasing in the last decades [29], showing the need for research and technological innovation for this productive chain.
Corroborating the findings of Rozane et al. [19] and Krug et al. [9], who applied CND to grapevines and citrus, respectively, there was no positive correlation between any single nutrient and yield, even in the high-yield population.This result evidences that a higher yield is the outcome of the appropriate balance among all nutrients.On the other hand, correlation analyses revealed a negative relationship between leaf B content and yerba mate yield.Moreover, in comparison to the high-yield group, the average B content was significantly higher in the low-yield group.Finally, when estimating the percentage of nutritional deficiency, balance, and excess by CND, it could be noticed that about onethird of the low-yield samples showed an excess of B, and it was the most imbalanced nutrient.It is known that B is a nutrient with a very narrow window between toxicity and deficiency [32].Although B toxicity is usually associated with saline soils in regions of low precipitation, it can also occur in highly acidic soils, especially if there is an excess of K or a deficiency of Mg [33].In general, adequate tissue B concentrations for the optimal growth of dicot plants range from 20 to 70 mg kg −1 [34].We estimated the normal range of B for yerba mate to be between 47 to 66 mg kg −1 .Thus, although yerba mate is a plant well-adapted to soil acidity, capable of tolerating high levels of Al and Mn, and also capable of accumulating high levels of B without showing visual symptoms of toxicity [35], elevated levels of B can be one of the factors that reduces its productivity.Even without causing direct toxicity, the high levels of B may be interfering with the absorption of other nutrients.For instance, Pawlowski et al. [36] found that excess B decreased the levels of Mg and Ca in soybean leaves.
More than 25% of the low-yield samples were classified as deficient in Cu and K. PCA, a multivariate technique that reduces the data dimensionality while preserving the main contributors to variance [37,38], also revealed a similar trend in the levels of K and Cu in the yerba mate leaves.As in the case of B, the availability of Cu is linked to soil characteristics, which mainly depend on its parent material [39].Soils originating from basic rocks are naturally richer in metals, including Cu.Thus, considering the yerba mate production areas in Rio Grande do Sul, we can state that soils originating from basalts are richer in Cu than soils originating from dacites, which predominate above 600 m in altitude [16,40].Motta et al. [41], evaluating yerba mate under low input systems in the Brazilian southern states, observed Cu contents in the soil ranging from 0.6 to 10.7 mg dm −3 and in the leaf from 7.3 to 14.0 mg kg −1 .In our research, these ranges were even broader (Tables 1 and 2).Those authors emphasized the importance of the soil parent material in the chemical composition of yerba mate.Also, it is important to notice that soil may contain high levels of Cu, but part of the Cu can be retained in the root system of plants, reducing the amount transported to the aerial part [42].
Regarding K, the normal range was estimated by CND to be between 11.9 and 16.7 g kg −1 (Table 11).As a result, 27% of the low-yield group samples had insufficient K levels.The first study aiming to define the sufficiency range for K in yerba mate was carried out by Reissmann et al. [43], who estimated that the adequate foliar K contents would be between 14 and 18 mg kg −1 .In the same sense, in a study that evaluated the effect of potassium fertilization on yerba mate, Santin et al. [44] identified that the leaf content of 16 mg kg −1 was sufficient to obtain high yields.These values are close to those estimated in the present study.The demand for K by yerba mate can be calculated based on the analysis of nutrient exports at harvest.Considering the K export potential in the yerba mate harvest as 15.05 g kg −1 , determined by Souza et al. [45], and the average yerba mate yield of 13.07 Mg ha −1 obtained in the present work, it is possible to estimate that about 197 kg K ha −1 , equivalent to 237 kg K 2 O ha −1 , are exported in a harvest.This value is proportionally higher as the yield increases, demonstrating that yerba mate is a K-demanding crop.Thus, monitoring K availability in the soil, as well as the nutritional condition of the plants, is fundamental for obtaining and maintaining high yield levels.
Significant correlations were found between the CND indices and the respective nutrient concentrations in plant tissue, but not between productivity and the CNDr 2 index (CND-r 2 = −0.1344yield + 13.856; R 2 = 0.0154; p = 0.111).A similar result was obtained by Rozane et al. [19] and Krug et al. [9].Those authors concluded that the nutritional balance would explain only part of the productivity of the vineyards and citrus orchards they evaluated, which were of different cultivars.In our case, in addition to the genetic variability of yerba mate, the evaluated production areas have diverse management systems, with distinct forms of fertilization, plant densities, and pruning frequency, directly interfering with the yield of the orchards.Regardless, this research brings the first definition of CND norms, critical levels, and sufficiency ranges for the leaf nutrient diagnosis of yerba mate.It will contribute to the better nutrition management of a crop with high cultural importance, in regional terms, and of growing economic importance worldwide.

Conclusions
The CNDr 2 indices were effective in establishing the critical levels, sufficiency ranges, and the nutritional status of yerba mate concerning the leaf concentration of the nutrients at deficiency, balance, and excess concentrations.

Figure 1 .
Figure 1.Polos Ervateiros, the yerba mate-growing regions in the state of Rio Grande do Sul, Brazil.Figure 1. Polos Ervateiros, the yerba mate-growing regions in the state of Rio Grande do Sul, Brazil.

Figure 1 .
Figure 1.Polos Ervateiros, the yerba mate-growing regions in the state of Rio Grande do Sul, Brazil.Figure 1. Polos Ervateiros, the yerba mate-growing regions in the state of Rio Grande do Sul, Brazil.

− 1 1 * 13 Figure 2 .
Figure 2. Percentage of nutritional deficiency, balance, and excess in high-and low-yield groups of yerba mate, estimated by CND.Figure 2. Percentage of nutritional deficiency, balance, and excess in high-and low-yield groups of yerba mate, estimated by CND.

Figure 2 .
Figure 2. Percentage of nutritional deficiency, balance, and excess in high-and low-yield groups of yerba mate, estimated by CND.Figure 2. Percentage of nutritional deficiency, balance, and excess in high-and low-yield groups of yerba mate, estimated by CND.

Figure 2 .
Figure 2. Percentage of nutritional deficiency, balance, and excess in high-and low-yield groups of yerba mate, estimated by CND.

Figure 3 .
Figure 3. Mean values of nutrient indices in the studied yerba mate population (n = 167).

Figure 3 .
Figure 3. Mean values of nutrient indices in the studied yerba mate population (n = 167).

Table 1 .
Average physicochemical parameters of the sampled soils.

Table 2 .
Leaf nutrient content of yerba mate database.
* Pearson correlation coefficient is significant at p < 0.05.** Pearson correlation coefficient is significant at p < 0.01.

Table 4 .
Estimation of the yield inflection points of cumulative variance functions for row-centered log-ratios in the yerba mate database population.

Table 5 .
Correlation between leaf nutrient content and yerba mate low-yield population (n = 133).
* Pearson correlation coefficient is significant at p < 0.05.** Pearson correlation coefficient is significant at p < 0.01.

Table 6 .
Correlation between leaf nutrient content and yerba mate high-yield population (n = 34).
* Pearson correlation coefficient is significant at p < 0.05.** Pearson correlation coefficient is significant at p < 0.01.

Table 7 .
Average leaf nutrient content in high and low-yield populations of yerba mate.

Table 8 .
Principal component analysis (PCA) loadings performed on log-transformed data from low-yield group (n = 133).

Table 9 .
Principal component analysis (PCA) loadings performed on log-transformed data from high-yield group (n = 34).

Table 10 .
Compositional nutrient norms based on the high-yield group of yerba mate.

Table 11 .
Linear regression between nutrient concentrations in yerba mate leaves from high-yield group and their respective CND indices, coefficient of determination, critical levels (CLs), and sufficiency ranges (SRs) for each nutrient.