Diagnostic Methods to Assess the Nutritional Status of the Carrot Crop

Leaf analysis may identify nutrient deficiency or excess in plant tissue. The aim of this study was to determine reference values and generate nutritional diagnosis from the results of leaf analysis of the carrot crop (Daucus carota L.) by the methods of Critical Level, Sufficiency Range, Diagnosis and Recommendation Integrated System, and Compositional Nutrient Diagnosis.


INTRODUCTION
Interpretation of the results of leaf analysis may identify nutrient deficiency or excess in plant tissue, thus allowing adjustments in the crop fertilization program (Malavolta et al., 1997).Despite the importance and application of leaf analysis, adequate contents for carrot (Daucus carota L.) growing found in the literature are out of date as they do not consider high yielding crops (Malavolta et al. 1997;Hanlon and Hochmuth, 2009).In the Alto Paranaíba region, yields are high over 80 Mg ha -1 (Aquino et al., 2015), twice that suggested in official recommendations for the state of Minas Gerais, Brazil (CFSEMG, 1999).
The carrot crop has a short cycle (90-130 days) and continually puts forth leaves during the growing season.There are no phenotypic changes over the market cycle (it does not include the reproductive phase).Thus, establishing a phenological stage to collect leaf samples is complex.The youngest fully expanded or physiologically mature leaf, around 40 days after sowing, as the diagnostic leaf (Malavolta et al., 1997).However, climatic conditions and the cultivar may lead to distinct growth and nutrient uptake characteristics after this period of cultivation.Nevertheless, definition of leaf physiologically maturity (recommended for sampling) is difficult because the period among leaf emergence is short.In order to overcome these problems, sample collection time can be established at harvest time (Hanlon and Hochmuth, 2009), and instead of sampling specific leaves, the whole shoot of the plants can be collected.The standardization of leaf collection time minimizes factors that may interfere in interpretation of analyses, such as nutrient dilution and concentration (Malavolta et al., 1997).
The methods most used for interpretation of leaf analysis have been Critical Level (CL), Sufficiency Range (SR), Diagnosis and Recommendation Integrated System (DRIS), and Compositional Nutrient Diagnosis (CND) (Coelho et al., 2013;Gott et al., 2014).The CL method consists of mathematical models to establish the leaf nutrient contents that allow production of a certain fraction (90, 95, or 99 %) of maximum yield.The SR establishes an optimal range of leaf content in accordance with the average content of nutrients in the leaves of the most productive plants (or fields).Although these are classical methods and represent the main information contained in the literature (Camacho et al., 2012;Santos et al., 2013), they are univariate and analyze each nutrient separately.
The DRIS and CND methods stand out, in comparison to the CL and SR methods, through being bivariate and multivariate, respectively.Diagnosis and Recommendation Integrated System was developed by Beaufils (1973), and it is based on comparison of indices calculated in accordance with the mutual relationship between two nutrients.An advantage of this method is minimizing dilution or concentration effects because it is based on a relationship of equilibrium between nutrients (Beaufils, 1973;Jones, 1981).Furthermore, this method enables identification of nutritional imbalances even when all nutrients are above the critical level (Baldock and Schulte, 1996).
The Compositional Nutrient Diagnosis (CND) method, developed by Parent and Dafir (1992), differs from DRIS because the content of each nutrient in the sample is corrected in accordance with the geometric mean of nutritional composition.Therefore, CND is based on multiple interactions among all the nutrients under diagnosis.
The CND method initially presents advantages compared to DRIS since it considers the interaction of all nutrients simultaneously, and not only dual interactions (Parent and Dafir, 1992).The DRIS and CND methods have been proposed for nutritional diagnosis in crops such as soybean (Urano et al., 2007), cotton (Serra et al., 2010), sugarcane (Santos et al., 2013), and common bean plants (Partelli et al., 2014)

MATERIALS AND METHODS
To determine reference values for diagnosis of nutritional status based on leaf nutrient content, a database with information from commercial carrot fields in the Alto Paranaíba region, MG, was created.Two hundred and ten commercial planting fields were sampled (155 in the winter and 55 in the summer) in the 2012 and 2013 crop seasons in areas ranging from 2 to 15 ha.The average duration of the crop growth cycle for summer and winter cultivars was 105 and 125 days, respectively.The average application rates of N, P 2 O 5 , and K 2 O were as follows: 118, 650, and 398 kg ha -1 , respectively.The main summer cultivars were Juliana and Poliana, and the winter cultivars were Baltimore, Belgrado, Maestro, Músico, Nancy, Nandrim, and Soprano.Winter is the preferred carrot growing season since it provides better weather conditions, which explains the larger area and number of cultivars.
Average altitude of the growing areas was around 1,100 m.The regional climate is classified as Cwa according to the Köppen-Geiger system, which is characterized by a dry season and a well-defined rainy season between October and March.The temperature and rainfall data during the sampling times for data collection are shown in the figure 1.
The soils were classifed as Latossolos Amarelos, Vermelhos and Vermelho-Amarelo (Oxisol), with clayey or very clayey texture.The vast majority of the fields under study have high fertility soils according to CFSEMG (1999).The physical and chemical properties of the soils are shown in table 1.
The database was composed of root yield and leaf nutrient content.Root yield in each field was obtained by harvesting four samples from 4 m of double rows, and their average weight was converted into tons per hectare.The time of sampling was that stipulated by farmers of the region based on growth, size standardization and firmness of the roots (an average of 105 days in the summer and 125 days in the winter).
Leaf collection consisted of the whole shoot at the time of harvest.Collection of the whole shoot was adopted to minimize sampling errors because the juxtaposed emergence of leaves hampers the individualization of each leaf position on the plant.Futhermore, the effects of nutrient dilution and concentration are minimized by collection at harvest time due to greater leaf dry matter (Hanlon and Hochmuth, 2009).
The leaves were washed and rinsed in deionized water, followed by drying in an air circulation laboratory oven at 70 °C for 72 h.They were then ground in a Willey mill equipped with a 1.27 mm sieve.Nutrient content was determined according to methods described by Malavolta et al. (1997).Contents of N, P, K, Ca, Mg, S, B, Cu, Fe, Mn, and Zn and root yields were critically analyzed (outlier analysis) and processed to generate reference values for leaf diagnoses by the CL, SR, DRIS, and CND methods.
To determine the CL, the yield data were organized in ascending order and grouped into 14 classes (number and classes = √n; in which "n" is the number of samples).
The maximum root yield obtained among the 14 classes was 111.7 Mg ha -1 and this was used to fit regressions of CL s .Thus, the CL was established for the yield of 100.5 Mg ha -1 , which is 90 % of the maximum yield.The interval composing the classes was obtained by dividing the total amplitude of the nutrient by the number of classes.Subsequently, linear mathematical models were fitted to study the relationship between the leaf content of each nutrient and the yield of the 14 classes, as performed by Gott et al. (2014).
Each nutrient was evaluated separately from others.
The significance of the parameters of the mathematical models was evaluated by the t test at 5 % probability.When the slope of the linear regression was significant by the t test and positive, the CL was estimated as the content corresponding to 90 % of the maximum yield (Gott et al., 2014).When the slope of the linear regression was negative and/or not significant, the CL was considered as the average content of the nutrient in the database.To establish standard scores or norms by the RS, DRIS, and CND methods, the population was divided into two subpopulations based on root yield.The threshold of these classes was defined as the average root yield plus 2/₃ of the standard deviation of this variable, which corresponded to 87.8 Mg ha -1 .After the data normality test, the subpopulation with the greatest yield was considered nutritionally balanced or a reference for establishment of norms/standard scores for interpretation of the methods.The subpopulation with high yield was represented by 64 fields (30.5 % of the fields) and low yield by 146 fields (69.5 % of the fields).
The average (x) and standard deviation (s x ) of the leaf nutrient contents were estimated in the high yielding population.The threshold levels of the SR were defined by the average content plus or minus the standard deviation of this variable (FS = x̄ ± S x̄) , as performed by Gott et al. (2014).As for the CL, each nutrient was evaluated separately from the others.
The mean (x) and variance (s²) of the dual relationships of the leaf nutrient contents were determined by the DRIS method, which composed the DRIS norms.The relationships were obtained from the ratio of macro-and micronutrient contents in g kg -1 and mg kg -1 , respectively.The functions of the dual relationships were obtained as described by Jones (1981) and chosen through the F method (Letzch, 1985): in which X = nutrient for which the index is to be calculated; Y 1 , Y 2 , ..., Y n = nutrients that are in the denominator of relationships with nutrient X; Z 1 , Z 2 , ..., Z m = nutrients that are in the numerator of relationships with nutrient X; Z/X = relationship between the contents of the nutrients Z and X of the sample to be subjected to DRIS; X/Y = relationship between the contents of the nutrients X and Y of the sample to be subjected to DRIS; z/x = average relationship between the contents of the nutrients Z and X provided by DRIS norms; x/y = average relationship between the contents of the nutrients X and Y provided by DRIS norms; and s = standard deviation of the dual relation in the reference population.
The DRIS index for each nutrient was calculated in accordance with the equation proposed by Beaufils (1973): in which X/Y and Y/X = functions of the dual relationships chosen for calculation of the DRIS index; m = number of functions in which the nutrient X is in the denominator; and n = number of functions in which the nutrient X is in the numerator.
For establishment of the Compositional Nutrient Diagnosis (CND), the proposal of Khiari et al. (2001) was adopted, which includes adaptations of the original model as described by Parent and Dafir (1992): in which V X = multinutrient variable for the nutrient X; X i = leaf concentration of the nutrient to be analyzed (mg kg -1 ); and G = geometric mean of dry matter constituents.
To calculate the geometric mean (G), the following equation was used: Rev Bras Cienc Solo 2016;v40:e0140813 in which G = geometric mean of dry matter constituents; X 1 , X 2 , ..., Xn = concentration of the nutrients in the dry matter (%); R = residue or concentration of the unmeasurable nutrients in the dry matter (%); and n = number of nutrients evaluated.
The CND index for each nutrient was calculated by the following equation: in which CND Xi = CND index for the nutrient; VX i = multinutrient variable of the field evaluated; VX´ = average of the multinutrient variable of the reference population; and S X´ = standard deviation of the multinutrient variable of the reference population.
For fields that composed the reference population, optimum leaf content was obtained by the DRIS and CND methods through fitting a regression equation between the nutrient index (independent variable) and its respective leaf content (dependent variable).Subsequently, nutrient content was obtained in the leaf that provides a null or balanced index.
The optimum range for each nutrient was established by leaf contents corresponding to ± ²/₃ indexes, the normal ranges of Beaufils (Urano et al., 2007;Santos et al., 2013).
The order of nutritional limitation was generated in two forms: in accordance with the optimum range of leaf content or in accordance with fertilization response potential.
From the optimum range, the order of nutritional limitation was estimated for the SR, DRIS, and CND methods, as well as from existing literature for purposes of comparison.
From the optimum content ranges of each method and from leaf nutrient concentration in the samples, the percentage of fields in imbalance (deficient or excessive) for each nutrient and diagnostic method was established.
Nutritional restriction was considered as deficiency or excess when the leaf nutrient content was below or above the optimum range established by diagnostic methods, respectively.The order of nutritional limitation by deficiency was generated in descending order of nutrients in regard to the frequency of fields under deficiency (DRIS and CND).
The order of nutritional limitation was also generated as a function of potential response to fertilization for DRIS and CND.In this case, the fields were first placed in five groups in accordance with fertilization response potential (FRP), namely: positive, positive or zero, zero, negative or zero, and negative (Wadt et al., 1998).
Positive FRP occurs when the DRIS or CND index of a nutrient is negative and greater in magnitude than the medium DRIS or CND.Positive or zero FRP occurs when the DRIS or CND index is negative and its module is larger than the mean DRIS or CND; however, it is not the nutrient with the lowest index.Zero FRP occurs when the DRIS or CND index is negative, but the module is lower than the mean DRIS or CND.Negative or zero FRP occurs when the nutrient index is positive and lower than the mean DRIS or CND.Negative FRP occurs when the nutrient index is positive and simultaneously the greatest of the indexes and greater than the mean DRIS or CND.Subsequently, the frequency of fields grouped into classes with positive and positive or zero response for each nutrient was calculated.In these classes of fertilization response, the nutrient can be considered as limiting for cultivation.
Field rankings in regard to nutritional status were performed by grouping the fields in three classes for each nutrient, namely: limiting by deficiency (LD), limiting by excess (LE) and not limiting (NL).The ranking was also generated by two different methods: as a function of the optimum content range and of the fertilization response potential.
In regard to the first method, those fields that presented leaf contents below or above the range limits were classified as LD or LE, respectively.When the nutrient content was in the range of the interval, it was considered NL for the field.The ranking of fertilization response potential was conducted for DRIS and CND.The fields with positive and positive or zero response were classified as LD for the respective nutrient.The fields with zero response were classified as NL, and the other classes of fertilization response (negative and negative or zero) were grouped in the LE class (Urano et al., 2006).
For both ranking methods (optimum range or fertilization response potential), the frequency of fields and average yield of each class (LD, LE, and NL) were presented for all nutrients, and leaf diagnosis methods (SR, DRIS, and CND) and data presented in the literature (Hanlon and Hochmuth, 2009) were evaluated.The t test was carried out for comparing the mean yields of the LD and LE classes with the yield of the NL class.
The calculations for determining reference values (content and optimum ranges) by the CL, SR, DRIS, and CND methods were carried out using Microsoft Excel® software.
The significance of the parameters of the fitted equations was analyzed by the t test at 5 % probability by SAS software 8.2.

RESULTS AND DISCUSSION
Analysis of chemical properties of the soil (Table 1) showed that acidity was corrected (pH in water 6.3) in addition to adequate P contents, available S, and exchangeable cations (average concentration of P and K in the very high fertility class, and medium contents of Ca, Mg, and S in the high fertility class), according to CFSEMG (1999).
In contrast, the soils showed imbalances in micronutrients since the average Mn content was classified as low, B as medium, Fe as adequate, and Cu and Zn as high, in accordance with the classification proposed by CFSEMG (1999).Moreover, the soils had medium contents of organic C (2 dag kg -1 ) and low remaining P -P-rem (10.6 mg L -1 ).The low P-rem value shows that the soils of the region are highly buffered for this nutrient, i.e., it is necessary to add large amounts to the soil to increase the P fraction available to the plants (Bedin et al., 2003;Broggi et al., 2011).
Regression equations with positive and significant slope for K, B, and Mn were fitted; thus, CL was established for these nutrients in the crop shoots.The CL for K was 49.7 g kg -1 , 56.2 mg kg -1 for B, and 79.9 mg kg -1 for Mn (Table 2).It was not possible to fit a regression equation with significant parameters for the CL of Ca, Mg, and Cu.Thus, the CLs of these nutrients were established as their average contents in the database.The higher proportion of roots/leaves in high yielding cultivars along with higher Ca partition to leaves may explain the lack of fit of a positive regression for Ca content.
High yielding plants may have lower Ca content in the leaves due to lower transpiration (smaller number of leaves) and higher allocation of Ca in the roots.Magnesium is allocated primarily in the roots and its content in leaves may decline with increasing yield (Cecílio Filho and Peixoto, 2013).The equations fitted to describe the leaf content of N, P, S, Fe, and Zn as a function of root yield were decreasing, i.e. with a slope lower than zero.Thus, the CLs of these nutrients were also defined as the average of their contents in the database.The fitting of decreasing functions means that there were dilution effects of the contents of these nutrients in the leaves as yield increased.
The root yield of the reference class for the SR, DRIS, and CND methods ranged from 87.8 to 120.5 Mg ha -1 , with an average of 98.2 Mg ha -1 .Averages and variances of the dual relationships (DRIS norms) of the high yielding population (>87.8Mg ha -1 ) are shown in table 3 and the CND norms in table 4.
Among the mathematical models to describe the relationship between indices (DRIS and CND) and nutrient contents in leaf (Tables 4 and 5), quadratic models obtained the best fit for cationic micronutrients (except for Zn in the CND index), whereas for the other nutrients, linear models were fitted.As for S, it was not possible to fit a significant mathematical model by the CND method (Table 4).
The contents and optimum ranges of the nutrients in the carrot shoots at harvest (Table 6) had little influence on the diagnostic method.In regard to optimum content, small variations were observed between adequate levels obtained by CL compared to DRIS and CND for cationic micronutrients.In this case, the CL established optimum contents with concentrations greater than the other methods.This is due to CL being more sensitive than CND to uncontrolled effects of environment, such as nutrient concentration or dilution in plant tissue (Wadt, 2008).Included in the calculation method of CND is the effect of variations in dry matter on interpretation of nutrient contents, correcting them and thus reducing the effects of dilution or concentration (Parent and Dafir, 1992).
An additional explanation for the difference in contents of cationic micronutrients between the methods is the high variability (coefficient of variation >46 % -data not shown) obtained for the content of these micronutrients in the plant tissue.Since the DRIS and CND methods are not univariate, i.e., they do not analyze each nutrient individually, they are not highly influenced by the variability of leaf nutrient concentrations.For the other nutrients (macronutrients and B), there were practically no differences between the optimum contents established by the CL, DRIS, and CND methods.
The sufficiency range of leaf nutrients, intervals and amplitudes were similar for the SR and DRIS methods, except for Fe (Table 6).The CND method showed narrower ranges that were included more in the optimum range established by SR and DRIS.For Fe, the SR established optimum range with an amplitude greater than DRIS, and these methods, for their part, established wider ranges than the CND.The smaller amplitudes of optimum ranges generated by CND favored more accurate diagnosis of the nutritional status of the carrot plants (Gott et al., 2014).However, small ranges generated for some nutrients may lead to wrong diagnosis due to error that arises in the analysis of determination of leaf content (laboratory errors).
Optimum ranges of leaf contents established by Hanlon and Hochmuth (2009) for carrot growing in the United States have higher upper limits for N and P compared to the limits of the present study.However, for K, Ca, and B, the optimum range stated in this study have higher limits (upper and lower) than those indicated by Hanlon and Hochmuth (2009).In regard to Mg, the optimum range proposed by these authors is wider and encompasses the range generated by SR, DRIS, and CND in this study.For Mn and Zn, the ranges obtained by SR, DRIS, and CND in this study were wider and encompassed the optimum ranges obtained in the United States.However, there was high discrepancy between the optimum ranges obtained by Hanlon and Hochmuth (2009) and those generated by the methods of SR, DRIS, and CND for Cu and Fe.For these two nutrients, the contents indicated by the authors are lower than those obtained in this study.This disagreement demonstrates the importance of having regionalized standard scores for leaf diagnosis, or generated for each growing condition (Serra et al., 2010;Gott et al., 2014).Similar conclusion regarding the importance of regionalization of optimum leaf nutrient levels for the Pera orange and sugarcane crops, respectively, was came by Camacho et al. (2012) and Santos et al. (2013).
The content and the optimum ranges generated by different methods (CL, SR, DRIS, and CND) were consistent for nearly all the nutrients in this study.Similar results were obtained by René et al. (2013), who proposed leaf diagnostic rates by the CL, DRIS, and CND  3. Average (AVER) and variance (VAR) of the dual relationships of leaf content in carrot shoots (DRIS) for the high yield (<87.8 Mg ha -1 ) population methods for poplar hybrid (Populus maximowiczii) and concluded that the methods generate consistent indices.Among the reasons that led to consistency in the ranges proposed by different diagnostic methods, the high yields of carrot obtained in the Alto Paranaíba region, MG, and the large number of fields sampled stand out.
The yield potential of the carrot crop is from 100 to 120 Mg ha -1 (Cecílio Filho and Peixoto, 2013).These high yields may be related to high soil fertility and the nutrient balance conditions of the crop.It is possible that most of these fields were nutritional balanced because to establish standard scores for leaf analysis by SR, DRIS, and CND, the reference population covered areas with yields higher than 87.8 Mg ha -1 (and reached 120.5 Mg ha -1 ).For common bean (Phaseolus vulgaris L.), Partelli et al. (2014) found that methods based on nutrient balance (DRIS and CND) showed high agreement in leaf diagnoses, though different from those obtained by SR.These authors evaluated 55 commercial fields and used those with yields higher than 2.7 Mg ha -1 as a reference population.The potential yield of the bean crop in Brazil estimated by the Blackman method is 6 Mg ha -1 (Oliveira et al., 2011).The low number of fields evaluated (compared to this study) and the limit yield of the reference class far below the maximum potential estimated for the crop, may have influenced the results obtained by Partelli et al. (2014).
There is consensus among the methods (SR, DRIS, and CND) that Mn and Mg are the most limiting nutrients for carrot growing in the general and low yield populations (Table 7).After these nutrients, K and Ca are considered limiting for production.However, these elements were not highlighted by all the diagnostic methods in sequence to Mn and Mg.Therefore, one can induce the need to increase K, Ca, and Mg supply for carrot cultivation.However, because of the limitation caused by Mn, the supply of K, Ca, and Mg should be made from sources that do not alter the pH of the soil.Increasing pH in this case would worsen the deficiency limitation posed by Mn.
The Soil Fertility Commission of Minas Gerais suggests at least 8 mmol dm 3 of Mg 2+ in the soil for carrot cultivation (CFSEMG, 1999).In the soils sampled, the Mg 2+ content was above this critical level (Table 1), and it was still the macronutrient most limiting to cultivation.In addition to content in the soil, another factor which must be taken into account is cation ratios.The average ratio of Ca:Mg and Mg:K was 3.2:1 and 3.5:1, respectively (Table 1).The limitation from Mg deficiency was greater than that of Ca and K for the carrot crop (Table 7).Future studies should investigate whether lower Ca:Mg and greater Mg:K ratios may provide nutritional balance conditions for Mg and increased yield.Among the factors that may have contributed to limitation of carrot growing by Mn deficiency, the low content of this nutrient in the soil may be highlighted (Table 1).
In addition, lime application resulted in an increase in pH, which, combined with high rates of P at sowing (average of 640 kg ha -1 of P 2 O 5 ), reduced Mn availability.High pH and high P application rates are factors that lead to reduced Mn availability to plants (Moreira et al., 2006;Gonçalves et al., 2011).
The importance of nutritional limitations caused by Fe and Zn was not significant, and the importance of B and Cu was moderate (Table 7).The lesser importance of these micronutrients compared to Mn can be explained by their levels in the soil (Table 1).
In accordance with the classification proposed by CFSEMG (1999), Zn and Cu contents are considered high in the soil, while those of B and Fe are moderate and good, respectively.
Analysis of the order of nutritional limitation (deficiency or excess) in the high yield population (>87.8Mg ha -1 ) does not designate the nutrients most limiting for carrot cultivation.The order of limitation varied as a result of the method of diagnosis.There was no agreement among the optimum ranges for most of the nutrients established in this study and those proposed by Hanlon and Hochmuth (2009).Thus, there was also no agreement among the orders of limitation obtained from the diagnostic indices generated in this study compared to the contents considered adequate by Hanlon and Hochmuth (2009).
The SR and DRIS methods showed higher frequency of fields as NL (non-limiting) compared to the frequency shown in the CND (Table 8).In contrast to this study, Serra et al. (2010) found that CND classified a higher number of fields as NL compared to DRIS for cotton plants (Gossypium hirsutum var.latifolium Hutch).In regard to the manners of interpretation (fertilization response potential [FRP] and optimum range) for DRIS and CND, the FRP classified a small number of fields as limiting (LD, limiting by deficiency or LE, limiting by excess) compared to the classification given by the optimum range of leaf contents.

Figure 1 .
Figure 1.Average temperature and monthly rainfall in the Alto Paranaíba region during collection for database formation.

Table 2 .
Leaf content corresponding to 90% of maximum yield (CL), mathematical model fitted to describe leaf content as a function of yield, and coefficient of determination of the models * and **: significant by the t test at 5 and 1 %, respectively.Rev Bras Cienc Solo 2016;v40:e0140813

Table 4 .
Possibly the use of bi-or multivariate methods compared to univariate methods in generating reference values in fields with high yield is not justified.High yields require such nutrient balance where effects of concentration, dilution, and relationship among nutrients are less expressive.Therefore, simpler methods, such as the CL and SR, can be used for accurate diagnosis.Moreover, the number of fields sampled Compositional Nutrient Diagnosis (CND) norms (average -AVE, and standard deviation -SD) of the multinutrient variables (Vi) and geometric mean of nutritional composition (G), mathematical model, amplitude of the indices, and coefficient of determination (R 2 ) of the regressions fitted to describe the leaf content of carrot shoot as a function of the CND index for the population of reference (<87.8Mg ha -1 ) ns : not significant by the t test; * and **: significant by the t test at 5 and 1 %, respectively

Table 5 .
Mathematical model, range of indices and coefficient of determination of regressions (R 2 ) fitted to describe the leaf contents in carrot shoots as a function of the Diagnosis and Recommendation Integrated System (DRIS) index for the high yield (<87.8 Mg ha -1 ) population **: significant by the t test at 1 %.Rev Bras Cienc Solo 2016;v40:e0140813 in this study (210 plots) may have contributed to obtaining similar results among the diagnostic methods evaluated.More representative sampling of carrot growing conditions (management practices and climate) may have reduced the influence of various effects that justify the use of bi-or multivariate methods.

Table 7 .
Order of nutritional limitation by deficiency for the carrot crop in the Alto Paranaíba region, MG, generated by the Compositional Nutrient Diagnosis (CND), Diagnosis and Recommendation Integrated System (DRIS), and Sufficiency Range (SR) methods, compared to the literature