Diversification and Soil Management Effects on the Quality of Organic Apricots

Nutritional and commercial value of the apricot fruits depends on the interaction among several factors, including the management strategies, the environmental conditions, and the genotype. This study aimed to evaluate the effect of orchard management and pedoclimatic conditions on apricot fruit quality traits (soluble solids, titratable acidity, dry matter and soluble solids and titratable acidity ratio) and phytochemical content (total phenolic content and antioxidant activity). Two apricot cultivars (Kioto* and Pieve*) were subjected to three different organic systems: (i) a standard management (BAU), (ii) an innovative diversified system with reduced tillage based on the natural cover conservation (INC) and (iii) an innovative diversified system with introduction of different cover crops (ICC), for two consecutive years (2019–2020). Overall, besides the effect of the climatic conditions, the multivariate analysis revealed the role of management systems, rather than genotype, on the quality of apricot fruits. In particular, INC soil management positively affected soluble solids, titratable acidity and dry matter parameters, especially in environments with low precipitation and high solar radiation (20.7 ± 1.7 g 100 g−1 FW, 19.9 ± 2.5 g MA 100 g−1 FW and 0.23 ± 0.2 g 100 g−1 FW respectively). In addition, the results indicate that ICC management showed the highest antioxidant values (1.45 ± 0.19 μg of Trolox equivalent (TE) 100 mg−1 FW), especially in the environment with high precipitation. This study was a useful tool for farmer to choose an effective organic system management in order to obtain fruits with high quality parameters.


Introduction
Apricot (Prunus armeniaca L.) is a fruit species widely spread in the temperate zones of the world, whose fruits are highly appreciated for their characteristic flavor, sweetness and juiciness. Italy is the fourth world's apricot producer after Turkey, Uzbekistan and Iran, and it is a leading country in Europe [1]. In 2019, the Italian apricot annual production was 279,900 tons (6.7% of the total world production) and the cultivated area covered 17,910 ha (3.1% of the world area), of which about 14.9% were organic [1]. The increasing market demand in production and consumption of fresh and dried apricots is mainly due to the positive effects on human health and disease risk reduction [2][3][4][5][6]. Indeed, apricots are rich in several phytochemicals such as sugars (mainly glucose, fructose and sucrose), organic acids (mainly malic, citric and quinic acids) [7], and bioactive compounds (phenolics, flavonoids, anthocyanins, carotenoids, dietary fiber, minerals and vitamins). All these compounds coupled with the fruit's physical and chemical characteristics are responsible for the apricot quality. In particular, phenols strongly contribute to the antioxidative potential of the fruit, suggesting a key role in the defense system against oxidative stress in human tissues [8].
Furthermore, consumers usually assess fruits for appearance, size, sweetness (sugar composition or sugar/acid ratio) and taste [9]. In addition, farmer's choices in terms of management and cultivars are done to satisfy consumers expectations in terms of quality and, especially for organic products, sustainability of production. Research activities should be then carried out by the active involvement of farmers in the research process, to implement their actual needs/research demand in the activities [10]. The modern apricot industry needs commercial cultivars characterized by highly fruit quality attributes [8], which can be intensified by the interaction of several factors, including the cultivation techniques under different environmental conditions and the genotype [3,11,12].
As for the agronomic management systems, conventional agriculture practices are still widely adopted by farmers, resulting in the highest rate of soil erosion and pollution [13][14][15]. However, in recent years, crop diversification, optimized use of energy, soil and water resources, as well as strategies to counteract adverse future climatic conditions are objective of growing attention, to contribute to long-term sustainability of the agroecosystems, climate change mitigation and adaptation [16]. In accordance, the recently published, and Farm to Fork strategy (FFS) of the European Commission [17] aims to a proportion of 25% organically managed farmland in the EU by 2030 towards more environmentally friendly and less intensive agricultural systems, with a significant reduction in fertilizer consumption and use of pesticides. The strategy refers to biodiversity, which with a separate strategy, "Bringing nature back into our lives" [18] is seen as crucial to increase the resilience of European societies. In this context, the introduction of agroecological practices can be the first step of this transformative process towards a more sustainable agriculture [19]. Sayara et al. [20] suggested to increase the orchard organic management, based on reducing the intensity of agronomic practices (e.g., tillage, fertilization, irrigation, etc.) and recycling of agricultural key compounds (nitrogen, phosporus potassium and/or organic matter) deriving from to the human activity and waste production, which can positively contribute to environmental sustainability. Moreover, the introduction of cover crops rotation, as well as intercropping, before or between principal crops as well as between shrubs or trees of plantation crops, may enhance soil chemical, physical and biological characteristics, contributing to soil health and sustaining cash crop development and yield [21,22].
In light of these considerations, our study aims to evaluate the effects of preharvest factors and orchard management on organic apricot fruit postharvest quality. In an Italian Long-Term Experiment (LTE), conducted throughout the direct engagement of local stakeholders (i.e., farmers and technician), cultivars-by-environment-by-management interactions were investigated to evaluate their potential impact in maximizing apricot quality, in terms of fruits pomological (color and weight), physico-chemical (soluble solids content, titratable acidity, dry matter and soluble solids and titratable acidity ratio) and nutraceutical characteristics (total phenolic content and antioxidant activity). The response of two cultivars (Kioto* and Pieve*) was compared under three different organic management systems based on different level of diversification for two consecutive years: a control system (Business As Usual) managed according to local standardized management practices for organic fruit orchards, an innovative diversified system with reduced tillage based on the natural cover conservation and use of an innovative diversified system with introduction of different cover crops. The study was carried out in 2019 and 2020 in the "maintenance of organic orchards" long-term experiment (MAIOR LTE) of CREA-Research Centre for Olive, Fruit and Citrus Crops, located in Rome, Central Italy (latitude 41.8000 • N, longitude 12.5690 • E, altitude 86 m a.s.l.). The MAIOR LTE experimental design was a split plot with 3 blocks. The main plot was assigned to the organic agronomic management, with three compared systems with a different level of diversification: (i) Business As Usual (BAU-a local common practices for organic fruit orchards management-soil tillage and organic commercial fertilizer), (ii) Innovative diversified system with Natural cover and a green Compost use (INC-no soil tillage at transplanting, natural cover, Municipal Waste Compost MCW), (iii) Innovative diversified system with introduced Cover and Compost use (ICCintroduction of Agroecological Service Crops, tillage and MWC). The main features of the three compared organic systems are reported in Table 1. Distribution once per year in autumn

Irrigation
Drip system (About 1500 m 3 ha −1 from June to September for 2 times a week in 12 h at a time) Drip system (About 1500 m 3 ha −1 from June to September for 2 times a week in 12 h at a time) Drip system (About 1500 m 3 ha −1 from June to September for 2 times a week in 12 h at a time)

Cover management
Natural cover, mowed twice per year (spring, autumn) Natural cover, mowed twice per year (spring, autumn) Natural cover (spring to autumn) mowed in autumn Cover crop mixtures and strips sowed in autumn and terminated in spring 1 BAU: Businness As Usual; 2 INC: Innovative diversified system with Natural cover and Compost; 3 ICC: Innovative diversified system with introduced Cover and Compost use.
The split plot was assigned to the apricot (Prunus armeniaca L.) cultivars. Two selfcompatible apricot cultivars with different vigor were chosen: Kioto* (KM) is a medium ripening (3th week of June in Central Italy) with contained vigor and high productivity, whereas Pieve* (PM) is a medium-late ripening cultivar (4th week of June in Central Italy) characterized by high vigor and good productivity. They were grafted on the same rootstock Myrobalan 29C, a clonal selection of Prunus cerasifera (able to contain the vigor of plants and suitable for all kinds of soils) and planted at spacing of 4.4 by 3 m, and an area pertaining to each tree of 13.2 m 2 . Each sub-plot was composed by 3 × 4 trees (132 m 2 ). The LTE started in winter 2017, at the end of a participatory process involving farmers in the identification of the studied factors to include the local stakeholders needs in the research activity [23].
Apricot samples were harvested at the homogeneous maturation grade (physiological ripening) evaluated by a tristimulus colorimeter (Chroma Meter CR-400, Minolta, Milan, Italy). Three sets of fruits (about 1.5 kg) were randomly sampled in each plot and three replicates in each plot were carried out for each cultivar. Quality analyses were carried out on fresh pitted fruits, whereas determination of phytochemical characteristics was performed on samples stored at −80 • C.

Climatic Conditions
Monthly data on mean air temperature at 2 m (AT_mean), average maximum temperature at 2 m (AT_max), average minimum temperature at 2 m (AT_min), soil temperature at 0.2 m (ST), relative humidity at 2 m (RH), solar radiation (SR), and rainfall (RF) were collected from weather station located nearby the experimental site during the 2-year period (from July 2018 to June 2020).
According to Köppen classification, Rome has a Mediterranean climate, characterized by mild winters and warm to hot summers. On average, January is the coldest month (7.5 • C), whilst July and August are the hottest (25.1 and 25.4 • C, respectively) and driest. November and December are the wettest months, with an average rainfall of 110-120 mm/month. During the apricot growing season (from March to June), the longterm mean air temperature recorded "1971-2000" varied from 10 • C to 23 • C, with an average minimum temperature above 5 • C and an average maximum temperature below 32 • C. Precipitation was equally distributed during March-May, with monthly average values equal to 75-80 mm, while June was slightly drier (45 mm).
The climate variables recorded during the two experimental years (from July 2018 to June 2020) are shown in Table 2. The monthly mean air temperature ranged between 6.0 • C (January 2019) and 26.4 • C (July and August 2019). The highest mean temperature was recorded in August 2019 (AT_max = 33.7 • C), whereas the lowest in January 2019 (AT_min = 3.3 • C). In winter, minimum temperatures rarely dropped below 0 • C, whereas in summer maximum temperatures were often over 30 • C, especially in 2019. The mean air temperature during the growing season was 16.1 • C for 2019 and 16.5 • C for 2020. Soil temperature followed the same trend as the air temperature: mean ST never dropped below 6.2 • C (January 2019) or rise above 26.4 • C (July 2018). Mean solar radiation varied from 53 W m −2 (December 2019) to 277 W m −2 (June 2019). The rainfall distribution differed greatly from the long-term trend: during 2019 and 2020 growing season, rainfall was significantly higher and lower than long-term averages (+16% and −51%, respectively). Moreover, May 2019 was an extremely rainy month (212.8 mm), while it never rained in the following month. In the same period of 2020, the total amount of rainfall was very low, especially in May (14.6 mm). At last, the maximum monthly average relative humidity was recorded in November 2019 (88.2%), corresponding to the maximum rainfall (302.6 mm), whereas the minimum RH in July 2019 (63.8%). Legend: AT_min = average minimum temperature; AT_max = average maximum temperature; AT_mean = mean air temperature; RH = relative humidity; SR = solar radiation; ST = soil temperature; RF = rainfall.

Soil Characterizations
Soils from each management system were sampled (0-25 depth) each year in early spring. The samples were dried for 48 h at 70 • C and then analyzed for Total Nitrogen (N) and Total Organic Carbon (TOC) contents by elemental LECO Nitrogen and TOC analyzers (mod. FP-528 and RC-612, respectively; St. Joseph, MI, USA) using Dumas and dry combustion methods, respectively. The main total nitrogen and total organic carbon soil content of the different systems management are presented in Figure 1.
Soils showed significant highest TOC content (p < 0.05) for the INC system in the period 2019-2020 (2.9 ± 0.3 and 3.2 ± 0.2 g 100 g −1 dry weight (dw), respectively). The lowest TOC content was observed for BAU and ICC during 2019 (2.6 ± 0.3 g 100 g −1 dw). As for N content, significant different values were observed between BAU and ICC (6.3 ± 0.4 and 7.0 ± 0.1 ppm dw, respectively), and INC (3.7 ± 0.6 ppm dw) during 2020. In general, ICC and BAU showed similar trends for N and TOC contents, whilst INC strongly differed from the others.

Fruit Quality Traits
Fruits were weighed, and the ground skin color on the external opposite sides (two readings in the equatorial perimeter) of all unpeeled fruits was assessed by the CIELab color space coordinates of a tristimulus colorimeter (Chroma Meter CR-200; Minolta, Milan, Italy), equipped with a D65 illuminant [9]. The instrument was calibrated with a standard calibration plate and the results, L*, a*, b*, were reported as means of the two opposite fruit sides. Similar procedure was applied to determinate the flesh color. An aliquot of fresh fruit was used to determine the dry matter (DM) using an oven (Tecnocalor 2000 model M60-VF Tecnovetro, Monza, Italy) calibrated at 105 • C ± 1, until a constant weight was reached. In order to measure the titratable acidity (TA), fruits (about 400 g per each cultivar, year and system management) were deprived of stones and homogenized using an automatic titration system (785 DPM Titrino, Metrohm Ldt, Herisau, Switzerland), according to Ceccarelli et al. [24]. The TA content was expressed as gram (g) of malic acid (MA) 100 g −1 fresh weight (FW) of edible fruits. Soluble solids content (SSC) was determined on fruit juice with a digital refractometer (Refracto 30PX, Mettler Toledo, Greifensee, Switzerland) and expressed as g 100 g −1 FW (Brix degrees).
Agronomy 2021, 11, x FOR PEER REVIEW 6 of 18 Figure 1. Soil total organic carbon (TOC) and total nitrogen (N) variability among the different system management during the three years. Legend: Values belonging the same parameters without common letters are statistically different according to Kruskal Wallis test (p < 0.05). Upper case letters refer to differences among years for each treatment; lower case letters refer to differences among treatments within the same year.
Soils showed significant highest TOC content (p < 0.05) for the INC system in the period 2019-2020 (2.9 ± 0.3 and 3.2 ± 0.2 g 100 g −1 dry weight (dw), respectively). The lowest TOC content was observed for BAU and ICC during 2019 (2.6 ± 0.3 g 100 g −1 dw). As for N content, significant different values were observed between BAU and ICC (6.3 ± 0.4 and 7.0 ± 0.1 ppm dw, respectively), and INC (3.7 ± 0.6 ppm dw) during 2020. In general, ICC and BAU showed similar trends for N and TOC contents, whilst INC strongly differed from the others.

Fruit Quality Traits
Fruits were weighed, and the ground skin color on the external opposite sides (two readings in the equatorial perimeter) of all unpeeled fruits was assessed by the CIELab color space coordinates of a tristimulus colorimeter (Chroma Meter CR-200; Minolta, Milan, Italy), equipped with a D65 illuminant [9]. The instrument was calibrated with a standard calibration plate and the results, L*, a*, b*, were reported as means of the two opposite fruit sides. Similar procedure was applied to determinate the flesh color. An aliquot of fresh fruit was used to determine the dry matter (DM) using an oven (Tecnocalor 2000 model M60-VF Tecnovetro, Monza, Italy) calibrated at 105 °C ± 1, until a constant weight was reached. In order to measure the titratable acidity (TA), fruits  Figure 1. Soil total organic carbon (TOC) and total nitrogen (N) variability among the different system management during the three years. Legend: Values belonging the same parameters without common letters are statistically different according to Kruskal Wallis test (p < 0.05). Upper case letters refer to differences among years for each treatment; lower case letters refer to differences among treatments within the same year.

Extraction and Total Phenolic Content and Antioxidant Activity Determinations
About 5 g of fresh whole fruit were homogenized with a blender (Ultra-Turrax T25, IKA Labortechnik, Staufen, Germany) in 25 mL of acidified methanol water solution (methanol/water 70/30 v/v with 5 mM HCl) to determine both total phenolic content (TPC) and antioxidant activity (AA), according to Ceccarelli et al. [24]. The extraction was carried out under shaking in a thermostatic bath at 37 • C for 2 h and centrifugation at 5 • C and 8000× g for 15 min (centrifuge mod. 4239R, ALC International, Milan, Italy). Finally, the supernatant was recovered. The TPC was determined by the Folin-Ciocalteu method, as reported by Ceccarelli et al. [24]. Results were calculated and expressed as mg of gallic acid equivalent (GAE) 100 g −1 FW. Determination of GAE was performed by GA standard curve (0.025-0.5 mg mL −1 ). Briefly, 0.4 mL of extract were added to 16.0 mL water, 2.0 mL of Folin-Ciocalteu phenol reagent, 6.0 mL of 1M sodium carbonate and the final volume was adjusted to 25 mL with the same solution used for the extraction. The TPC of the extracts was determined spectrophotometrically at 760 nm after 2 h using an Evolution 300 UV-Vis Spectrophotometer (Thermo Electron Scientific Instruments, Madison, WI, USA). In the end, AA of the extracts was determined by DPPH (1,1-diphenyl-2-picrylhydrazyl) method, as described by Ceccarelli et al. [24]. 1.5 mL of DPPH solution was added to 1.5 mL of fruit extract and, after 15 min, absorbance at 515 nm was determined using the UV-Vis Spectrophotometer. Trolox (0.5-10 µg mL −1 ) was used as reference compound and AA was expressed as µg of Trolox equivalent (TE) 100 mg −1 FW.

Statistical Analysis
Descriptive statistics, including mean, standard deviation (SD), minimum (Min) and maximum (Max) value, 25th percentile (Q1), median, and 75th percentile (Q3) were calculated for SSC, TA, DM, TPC, and AA. Generalized linear model (GzLM) was applied to assess the influence of three factors (cultivar, year and management system) and all their two-way interactions on SSC, TA, DM, TPC, and AA. In GzLM, the response variable Y has a probability distribution from the exponential family with mean E[Y] = µ. GzLM estimates the relationship between the mean value of the response variable and a linear combination of the explanatory variables X, which are connected by a link function g as follows: where g(µ) = [g(µ 1 ), . . . , g(µ n )] is the n × 1 vector of link functions, X is the n × (p + 1) predictors matrix, and β is the (p + 1) × 1 vector of unknown parameters.
In our study, we assumed that each of the response variables Y (SSC, TA, DM, TPC or AA) followed a gamma distribution and the logarithm of its expected value E(Y) could be modelled by a linear combination of parameters. The gamma distribution is commonly used for continuous, non-negative and skewed response variables. In our study, gamma regression models GzLM with the logarithm as the canonical link function and the Gamma distribution function were applied. In detail the statistical model was: where Y = SSC or TA or DM or TPC or AA; β 0 = average effect common to all observations; τ i = cultivar (CV), i = 1, 2 (1 = KM, 2 = PM); δ l = year, l = 1, 2 (1 = 2019, 2 = 2020); ζ k = management system (M), k = 1, 2, 3 (1 = BAU, 2 = ICC, 3 = INC); ε = error term.
Parameters in GzLM are estimated using maximum likelihood. Hypothesis tests are constructed using Wald-type statistics based on the maximum likelihood estimators. The overall fit of GzLM model could be measured by the deviance statistic, i.e., twice the difference between two log-likelihoods, which can be approximated by a chi-square distribution for sufficiently large sample sizes. Therefore, a Likelihood Ratio Chi-Square was used to test the goodness-of-fit. The deviance generalizes the sum of squared errors in normal linear models.
Finally, Principal Component Analysis (PCA) was performed on the quality attributes and nutraceutical parameters (SSC, TA, DM, TPC or AA) in order to highlight similarities among apricot samples and clustering tendencies for cultivars and management systems.
Non-parametric Kruskal Wallis test was also applied to evaluate significant differences in TOC and N values among management systems during the 2018-2020 period.

Fruit Apricot General Aspects
The three tested factors (cultivar, year and organic system management) showed no significant effect on apricot fruits weight and color. Mean values of fresh weight ranged from 23.6 g to 49.9 g for Pieve* and from 23.5 g to 48.3 g for Kioto*, highlighting a high variability inside each cultivar (Table S1). The variability found in each cultivar was clearly higher than that due to the harvest years and system management. Regarding skin color (Table S1), coordinates showed a slight variability among the organic system managements, harvest years and cultivars (overall mean values L* 51.32 ± 3.31, a* 19.84 ± 3.24, b* 32.24 ± 3.64), confirming similar apricot ripening state. Analogous behaviors were also observed for fruit flesh color coordinates (overall mean values L* 54.82 ± 2.72, a* 12.01 ± 3.04, b* 40.23 ± 3.14).

Cultivars Variability
As previously reported, the CIELab coordinates were measured in field on fruits skin to assess the ripening stage and successively on the flesh for both harvest experimental  Table S1.
Regarding other qualitative parameters, the cultivars considerably differed in their SSC and DM values ( Table 3). The soluble solids content was higher for Pieve* than Kioto* (median 18.5 g 100 g −1 and 16.7 g 100 g −1 , respectively), as well as dry matter (median 0.21 g 100 g −1 FW and 0.19 g 100 g −1 FW, respectively). TA amount was similar for the two cultivars, although Kioto* was slightly higher (median 20.3 g MA 100 g −1 FW for KM and 19.6 g MA 100 g −1 FW for PM). High SSC/TA ratio was calculated for the two cultivars: Pieve* ranged from 0.89 to 1.09 and Kioto* between 0.78 and 0.98 (Table 3). A high variability in the total phenol content and the antioxidant activity was observed for both apricot cultivars (Table 3). Indeed, TPC ranged from 211 to 528 mg GAE 100 g −1 FW for Kioto* and from 273 to 565 mg GAE 100 g −1 FW for Pieve*, whereas AA varied between 0.80 and 1.83 µg TE 100 mg −1 FW for Kioto* and between 0.81 and 1.70 µg TE 100 mg −1 FW for Pieve*.

Management Systems Variability
The organic management systems had no significant effect on color coordinates CIELab, SSC and TA. The lowest SSC (median 17.4 g 100 g −1 ) and the highest TA (median 21.5 g MA 100 g −1 FW) values were observed when ICC was applied, contrary to what was recorded for BAU (Table 3). Consequently, the lowest SSC/TA ratio was observed in fruits harvested when ICC management was used (0.8), whereas the highest one was found in BAU (1.0) (Table S2). Likewise, DM was lower for the ICC system than for BAU and INC. Significant differences and high variability were observed for the total phenolic content and the antioxidant activity. TPC for ICC system was greatly higher than BAU and INC (median 462 mg GAE 100 g −1 FW (ICC PM), 449 mg GAE 100 g −1 FW (ICC KM), 411 mg GAE 100 g −1 FW (BAU PM), 230 mg GAE 100 g −1 FW (BAU KM), and 337 mg GAE 100 g −1 FW (INC PM), 306 mg GAE 100 g −1 FW (INC KM)). Finally, ICC and BAU systems showed higher AA values than INC (median 1.14, 1.14 and 0.99 µg TE 100 mg −1 FW, respectively). Table 3. Descriptive statistics, including mean, standard deviation (SD), minimum (Min) and maximum (Max) values, 25th percentile (Q1), median, and 75th percentile (Q3), for quality parameters, as a function of cultivar (Kioto* KM and Pieve* PM) and agronomic treatment (Business As Usual-BAU; Innovative diversified system with Natural cover and Compost use-ICC; Innovative diversified system with introduced Cover and Compost use-INC). Legend: SSC = soluble solids content (g 100 g −1 ), TA = titratable acidity (g MA 100 g −1 FW), DM = dry matter (g 100 g −1 FW), TPC = total phenolic content (mg GAE 100 g −1 FW) and AA = antioxidant activity (µg TE 100 mg −1 FW).

Year Variability
Fruit flesh and skin did not show significant differences in lightness, blue-yellow and green-red values between the two harvest years (Figure 2).  Apricot mean weight showed no significant differences between the two years (from 23.6 to 50.0 g in 2019 and from 24.1 to 49.5 g in 2020), while fruits harvested in 2020 (year characterized by the lower rainfall and relative humidity) exhibited the highest SSC and DM values (Figure 2). SSC ranged from 17.4 to 24.1 g 100 g −1 in 2020, and from 13.6 to 20.8 g 100 g −1 in 2019. The DM varied between 18.1 and 28.3 g 100 g −1 in 2020 and between 13.6 and 24.0 g 100 g −1 in 2019. Quite similar TA values were observed for the two years ( Figure 2). Furthermore, average SSC/TA ratio for 2019 showed the highest values (1.0) respect to 2020 (0.8).

Generalized Linear Models
As a first step, normality of the data distributions for all parameters was checked by Shapiro-Wilk test. Because the test responses indicated that data were not normal distributed and considering that all distributions were continuous, non-negative and skewed, generalized linear models based on gamma distribution for each response variable were applied, in order to assess the overall influence of three factors (cultivar, year and management system) on all variables. No significant differences were observed for fruit fresh weight and CIELab color coordinates of fruit skin and flesh. Therefore, these parameters were not considered for GzLM. Table 4 showed the results of the analysis concerning the effect of cultivar, organic management system and year on five parameters. GzLM model for SSC showed that two single factors (cultivar and year) were statistically significant (p = 0.0002 and p < 0.0001, respectively), as was their two-way interaction (p = 0.0002). Legend: CV = cultivar M = management system. Numbers in bold were significant (p < 0.05).
Titratable acidity was not affected by any of the three factors. Regarding DM, all factors resulted highly significant (p < 0.0001), whether as single factors or in the two-way interactions, with the exception of the cultivar × treatment interaction (p = 0.6113). Agronomic treatment was the most significant factor for TPC (χ 2 = 11.76, p = 0.0006). The year, although statistically significant (χ 2 = 4.78, p = 0.0288), appeared less important than treatment, while the treatment × year interaction showed a higher weight (χ 2 = 11.75, p = 0.0006). The genotype did not affect AA, whereas agronomic treatment was the main contributing factor to AA (χ 2 = 12.87, p = 0.0003), followed by year (χ 2 = 11.70, p = 0.0003) and their two-way interaction (χ 2 = 12.88, p = 0.0003).

Clustering Tendencies
Principal Component Analysis was used to visualize simultaneously the cultivars, the management systems and harvest years on the base of their chemical properties and antiradical potential similarity and to highlight those metabolites contributing mostly to their grouping. The first two principal components (PC1 and PC2) accounted for 49.9% and 23.4% of the total observed variance, respectively. PC1 was positively correlated with SSC, AA and DM content, whereas PC2 was positively correlated with TPC and TA content (Figure 3). PCA emphasized the different response of the two cultivars and the three organic management systems over the two years. Indeed, the biplot showed two clusters closely related to the environmental conditions (i.e., years). The first group was on the right side and depicted the behavior of the apricot samples in 2020, characterized by high values of SSC, DM and AA. The second was on the left side and represented the behavior of the samples in 2019, with low values of SSC, DM and AA.

Discussion
Apricot is one of the most important fruit species grown in the world for its satisfying energetic, dietetic, nutritional and nutraceutical values, linked to different chemical compounds which give a good appearance and a high consumer acceptability [25]. In this study, the characterization of two cultivars (Kioto* and Pieve*), aiming to evaluate the effects of different organic system managements, was carried out through chemical and statistical analyses of nutraceutical and pomological traits. In addition, the two-year study, coupled with chemometric processing, allowed evaluating the effects of weather variability on qualitative traits. Our results confirmed the strong relationship between fruit quality parameters and climate and genotypes, as reported by several authors [5,[25][26][27]. In accordance with Leccese et al. [5], our study highlighted no significant differences of flesh and skin CIELab coordinates respect to monitored years and management systems. Indeed, in order to reduce the effects of different fruit maturation on qualitative parameters, fruit were selected at similar ripening stage based on homogeneous CIELab coordinates, as suggested by Bureau et al. [28] Contrary to what reported by Leccese et al. [5] and Bartolini et al. [29], the restricted variation recorded between the two tested cultivars for qualitative parameters, such as fruit fresh weight, titratable acidity TPC and AA, could be due to the discrete genetic proximity of Kioto* and Pieve*, as reported by Salazar et al. [30]. At harvest, both cultivars showed a SSC content higher than 11.2 g 100 g −1 FW, corresponding to the minimum degree of ripeness at packing, according to the mean value below the threshold indicated by the European Commission about apricots marketing (Commission Regulation (EC) n. 112/2001) [31]. For what concerns the SSC, on both cultivars we found higher values than the means reported in literature [32][33][34], reasonably due to the specific genotypes and climatic/management conditions of the trials. Moreover, the

Discussion
Apricot is one of the most important fruit species grown in the world for its satisfying energetic, dietetic, nutritional and nutraceutical values, linked to different chemical compounds which give a good appearance and a high consumer acceptability [25]. In this study, the characterization of two cultivars (Kioto* and Pieve*), aiming to evaluate the effects of different organic system managements, was carried out through chemical and statistical analyses of nutraceutical and pomological traits. In addition, the two-year study, coupled with chemometric processing, allowed evaluating the effects of weather variability on qualitative traits. Our results confirmed the strong relationship between fruit quality parameters and climate and genotypes, as reported by several authors [5,[25][26][27]. In accordance with Leccese et al. [5], our study highlighted no significant differences of flesh and skin CIELab coordinates respect to monitored years and management systems. Indeed, in order to reduce the effects of different fruit maturation on qualitative parameters, fruit were selected at similar ripening stage based on homogeneous CIELab coordinates, as suggested by Bureau et al. [28] Contrary to what reported by Leccese et al. [5] and Bartolini et al. [29], the restricted variation recorded between the two tested cultivars for qualitative parameters, such as fruit fresh weight, titratable acidity TPC and AA, could be due to the discrete genetic proximity of Kioto* and Pieve*, as reported by Salazar et al. [30]. At harvest, both cultivars showed a SSC content higher than 11.2 g 100 g −1 FW, corresponding to the minimum degree of ripeness at packing, according to the mean value below the threshold indicated by the European Commission about apricots marketing (Commission Regulation (EC) n. 112/2001) [31]. For what concerns the SSC, on both cultivars we found higher values than the means reported in literature [32][33][34], reasonably due to the specific genotypes and climatic/management conditions of the trials. Moreover, the SSC/TA ratio was also evaluated, because it is a good indicator of ripening and contributes to the overall organoleptic quality: high values represent a high perception of sweetness for consumer [35]. On average, higher SSC/TA values than the maximum value reported by Amoriello et al. [36] for seven widely spread cultivars were observed. The high recorded dry matter content of both tested cultivars suggested their suitability for processing, especially for drying [26].
According to Bartolini et al. [37] and Tarantino et al. [38] reporting high dependence of SSC, DM, TPC and AA accumulations in apricot fruits on the climatic conditions (air temperature and rainfall rate, in particular), our result confirmed the strong effect of the years and, thus, weather variability on qualitative parameters, with the exception weight, colors and TA. Actually, it is well known how the climate and weather variability holds an important role on the qualitative and nutritional parameters of vegetables and fruits [37][38][39][40]. Bartolini et al. [37] reported a key role of water availability on fruit quality traits, highlighting a lowest SSC and DM values and a higher TA content on fruits grown and ripened during warm-wet spring season. Similar results were observed in our study, confirming low SSC and high TA values during 2019-2020 (characterized by less abundant rainfall). Consequently, also the SSC/TA ratio was strongly influenced by water availability. Furthermore, beside water availability, light intensity and air temperature are related to the antioxidant activity in different fruit species [41]. Contrary to the results of Bartolini et al. [36], in our organic trial apricot fruits had the lowest antioxidant levels and the highest TPC values if high amounts of rainfall occurred during the growing period. This apparent discrepancy could be due to the measure methodology used for the antioxidant activity evaluation (Ferric Antioxidant Power FRAP, 2,2 -Azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) diammonium salt Radical Scavenging Assay ABTS, DPPH etc.), that strongly influenced the evaluation and could represent a limit [42]. Conversely, our results were in accordance with Sahamishiirazi et al. [43], reporting the highest TPC level when the highest humidity occurred during fruit growing period. In detail, high humidity could influence the transpiration and photosynthetic activities and consequently results in high TPC synthesis [43]. Sahamishirazi et al. [43] and Solovchenko and Schmitz-Eiberger [44] also showed that the antioxidant biosynthesis depends on temperature and spectral properties of light. High solar radiation and temperatures could contribute to reach high AA values, as it happened during the harvest year 2019-2020.
The management system adopted on fruits trees could strongly influence the quality parameters [45]. In particular, the different organic management options for orchards might impact the biogeochemical processes and soil properties (e.g., pH, soil biotic activity, organic matter mineralization) dealing with the function of supplying nutrient to the plant [45]. Different intake of nutrients, such as phosphorus and nitrogen, could significantly affect the biosynthesis of fruits SSC and TPC, consequently causing a strong variation of AA [46,47]. Many studies showed an increase of TPC and SSC content when lower nitrogen availability due to the slow-release property of organic nutrient input occurs [46,47]. We observed similar trend in our study: the TPC values for INC system increased during 2020 when soil total nitrogen decreased, as showed in Figure 1. Beside the role of tillage, this effect showed not being related to the fertilizer, with similar trend for BAU and ICC (Figure 1). This result could be related to the low C:N ratio (Table 1) of the used municipal waste compost, leading to a fast N release into the soil, similarly to the commercial fertilizer. However, the different compounds biosynthesis could be due to upon activation of defense mechanisms which could favor the accumulation of these compounds [47]. Moreover, there is evidence that differences in fertilization regimens (different N input) are responsible for higher TPC content in organic harvests [48]. The reduction of plant nitrogen availability could be due to also other factors such as soil temperature and rainfall amount, which affect the transformation of nitrogen into forms that can be assimilated by the plant or the washout. In our study, the higher rainfall associated with low soil temperature that occurred during 2019 could have led to a loss in TOC as reported in Figure 1, in agreement with Nie et al. [49]. Consequently, apricot fruit qualitative parameters differed among the organic management systems in the two harvest years. Nowadays, little is known about the possible impact of conservation agriculture on the content of fruits phytochemicals in relation to the different organic management system. However, our results highlighted a synergic effects of weather conditions and soil management. No tillage management coupled with municipal waste compost fertilization and natural cover in environments with low precipitation and high solar radiation, had particular effect on SSC, TA and DM parameters. Unlike our results, Massantini et al. [50] reported no effect of tillage on tomato SSC and TA values and significant differences for chemical traits when different cover crop was used. Furthermore, high TPC and antioxidant content in apple fruits were observed by Ceccanti et al. [51] when soil was no tilled. On the contrary, we found the highest values when soil tillage was used coupled with municipal waste compost fertilization and cover crop mixtures. According to several studies, high levels of health promoting phytochemicals in fruit productions could be linked to several environmental conditions, such as greater plant stress, peculiar rhizosphere microbial communities, and/or lower available nitrogen [3,16,21,46,51,52]. Actually, in our study the introduction of cover crops in the ICC system could have resulted in a short-term shortage of nutrients' availability in the soil for the fruit crop, since competition for resources during the cover crop growing period (November-April) could have occurred [53,54]. The utilization of Italian ryegrass or oat (highly competitive cereals [55]), in the cover crop mixture, beside the provision of ecosystems services such as the control of the spontaneous flora development, organic matter inputs, soil protection and nutrient leaching prevention during the rainiest period, could have interfered with the apricot plants, stimulating the defense response of the crop, finally benefiting the fruit phenolic content. This result is supported by the mean coverage of the cover crop evaluated at termination (70.8%, standard deviation 20.6) and confirms the close linkage between the system management and fruit product quality and its nutritional value. It also suggests the need for careful management of complexity (e.g., cover crop species choice in the mixtures composition) to drive system towards mitigation of potential disservices (i.e., potential production loss due to competition with cover crops) in favor of services (i.e., increase of TPC content of the fruits).

Conclusions
This paper provides information regarding the effects of different organic management systems during two-year trials on fruit quality traits (SSC, TA, DM, and SSC/TA) and phytochemical content (TPC and AA) of two widely spread apricot cultivars. Multivariate analysis revealed a greater influence of harvest year on SSC and DM, whereas TPC and AA appeared to be mostly affected by the system management.
Nowadays, the strategy to increase soil fertility, the related fruit quality while preserving the environment is based on the use of renewable resources (e.g., compost), recycle of key compounds (nitrogen, phosphorus potassium or organic matter) and reduced tillage that may minimize the environmental hazards and maximize the ecological benefits. In this context, INC management turned out to be really interesting in regard to SSC, TA and DM parameters, especially for environments with low precipitation and high solar radiation. As well-known, apricots with high SSC and DM are commonly used for processing, such as juice, dried fruit or jam. As regards antioxidant activity, ICC management showed the highest values, especially for environment with high precipitation. In this context, PCA was a useful tool to evaluate simultaneously the pomological and nutraceutical parameters respect to the cultivars, harvest years and organic management systems. This information might be useful for farmers to choose the best organic system management in relation to the long-term meteorological data, in order to obtain fruits with high pomological and nutraceutical potential.

Supplementary Materials:
The following are available online at https://www.mdpi.com/article/10 .3390/agronomy11091791/s1, Table S1: Means data, standard deviation (SD), minimum (Min) and maximum (Max) value of fruits weight and color coordinates (CIElab) divided by cultivar (Kioto* KM and Pieve* PM), agronomic treatment (Business As Usual -BAU; Innovative diversified system with Natural cover and Compost use-ICC; Innovative diversified system with introduced Cover and Compost use-INC) and years. Table S2: Means data, standard deviation (SD), minimum (Min) and maximum (Max) value of SSC:TA divided by cultivar (Kioto* KM and Pieve* PM), agronomic treatment (Business As Usual-BAU; Innovative diversified system with Natural cover and Compost use-ICC; Innovative diversified system with introduced Cover and Compost use-INC) and years.